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Abstract 


In  this  thesis,  material  anisotropy  and  the  phenomena  of  void  formation  and  growth  (cavitation) 
in  nonlinearly  elastic  incompressible  solids  are  considered.  We  first  discuss  materials  which  are 
transversely  isotropic  in  the  contexts  of  linear  and  nonlinear  hyperelasticity,  and  derive  important 
constitutive  restrictions  on  the  nonlinear  stored-energy  function,  W,  by  relating  W  to  the 
associated  infinitesimal  elastic  moduli.  We  then  propose  a  class  of  stored-energy  functions 
satisfying  these  conditions  to  model  incompressible  transversely  isotropic  nonlinearly  elastic 
materials.  The  effect  of  material  anisotropy  on  void  nucleation  and  growth  in  incompressible 
nonlinearly  elastic  solids  is  then  examined.  A  bifurcation  problem  is  considered  for  a  solid 
sphere  composed  of  an  incompressible  homogeneous  nonlinearly  elastic  material  that  is 
transversely  isotropic  about  the  radial  direction.  Under  a  uniform  radial  tensile  dead-load,  a 
branch  of  radially  symmetric  configurations  involving  a  traction-free  internal  cavity  bifurcates 
from  the  undeformed  configuration  at  sufficiently  large  loads.  Closed  form  analytic  solutions  are 
obtained  for  a  specific  anisotropic  material  model.  In  contrast  to  the  situation  for  an  isotropic 
sphere,  bifurcation  here  may  occur  locally  either  to  the  right  (supercritical)  or  to  the  left 
(subcritical),  depending  on  the  degree  of  anisotropy.  In  the  latter  case,  the  cavity  has  finite  radius 
on  first  appearance.  Such  dramatic  cavitational  instabilities  were  previously  encountered  by 
Antman  and  Negron-Marrero  [4]  for  anisotropic  compressible  solids  and  by  Horgan  and  Pence 
[20]  for  composite  incompressible  spheres.  Finally,  the  effects  of 
material  anisotropy  and  inhomogeneity  on  void  nucleation  and  growth  in 
incompressible  anisotropic  nonlinearly  elastic  solids  are  examined  via  a  similar  bifurcation 
problem  for  a  composite  sphere  whose  phases  are  each  composed  of  an  arbitrary  homogeneous 
incompressible  transversely  isotropic  material.  Several  types  of  bifurcation  are  found  to  occur. 
Explicit  conditions  determining  the  type  of  bifurcation  are  established  for  the  general  transversely 
isotropic  composite  sphere.  In  particular,  if  each  phase  is  described  by  an  explicit  anisotropic 
material  model  (a  generalization  of  the  classic  neo-Hookean  model  to  anisotropic  materials), 
phenomena  which  were  not  observed  for  the  homogeneous  anisotropic  sphere  nor  for  the 
composite  neo-Hookean  sphere  may  occur.  The  stress  distribution  is  also  examined. 
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Chapter  1 

INTRODUCTION 
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Void  nucleation  and  growth  in  solids  play  a  fundamental  role  in  fracture  and  other 
failure  mechanisms.  (See  e.g.  Tvergaard  [42]  for  a  recent  review  of  void  growth  in 
metals.)  Sudden  void  formation  ("cavitation")  in  vulcanized  rubber  has  also  been 
observed  experimentally  by  Gent  and  Lindley  [14];  see  also  Williams  and  Schapery  [45]. 
A  recent  review  on  cavitation  in  rubber  is  that  of  Gent  [13].  Nonlinear  theories  of  solid 
mechanics  have  been  extensively  used  recently  to  model  such  phenomena.  The  impetus 
for  many  of  these  developments  has  been  supplied  by  the  theoretical  work  of  Ball  [5]  on 
a  class  of  bifurcation  problems  for  the  equations  of  nonlinear  elasticity  that  model  the 
appearance  of  a  cavity  in  the  interior  of  a  solid  homogeneous  isotropic  elastic  sphere  or 
cylinder  once  a  critical  external  load  is  attained.  An  alternative  interpretation  for  such 
problems  in  terms  of  the  sudden  rapid  growth  of  a  pre-existing  microvoid  has  been  given 
by  Horgan  and  Abeyaratne  [19];  see  also  Sivaloganathan  [32].  As  pointed  out  for 
example,  in  [19],  cavitation  is  an  inherently  nonlinear  phenomenon  and  cannot  be 
modeled  using  linearized  solid  mechanics  theories. 

In  the  comprehensive  work  of  Ball  [5]  on  radially  symmetric  solutions,  bifurcation 
and  stability  analyses  are  carried  out  for  displacement  and  traction  boundary-value 
problems  in  n-dimensions  for  both  incompressible  and  compressible  isotropic  materials. 
For  incompressible  materials  the  results  are  extensive,  while  those  in  the  more  difficult 
compressible  case  are  comparatively  limited  and  require  several  constitutive 
assumptions.  Further  studies  in  the  compressible  case  were  carried  out  in  [39, 31,  19,  32, 
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33,  10,  28,  18,  40].  Anisotropic  compressible  materials  were  considered  by  Antman  and 
Negron-Marrero  [4].  Other  contexts  in  which  cavitation  for  compressible  materials  was 
investigated  include  consideration  of  non-radially  symmetric  solutions  [26], 
elastodynamics  [30],  and  elastic  membrane  theory  [17,  37];  (see  also  [16]  for  a 
membrane  made  of  an  incompressible  material).  For  incompressible  materials,  finite 
strain  plasticity  models  were  investigated  in  [9]  while  the  effects  of  rate  dependence  were 
examined  in  [1].  Further  studies  for  incompressible  materials  were  carried  out  in  [6,  24] 
for  elastostatics  and  in  [7]  for  elastodynamics.  The  effects  of  material  inhomogeneity  on 
cavitation  in  incompressible  materials  were  investigated  by  Horgan  and  Pence  [20-22]. 
Void  collapse  for  both  incompressible  and  compressible  materials  has  been  examined  in 
[2],  Further  work  in  plasticity  was  carried  out  in  [23, 25, 43]. 

The  purpose  of  this  dissertation  is  to  examine  the  effects  of  material  anisotropy  on 
cavitation  for  incompressible  nonlinearly  elastic  solids.  Two  previous  investigations  [4, 
20]  provided  particular  motivation  for  this  work.  For  classes  of  compressible 
transversely  isotropic  nonlinearly  elastic  cylinders  and  spheres,  Antman  and  Negron- 
Marrero  [4]  found  dramatic  material  instabilities.  In  particular  they  showed  that 
solutions  describing  cavitation  differ  significantly  from  those  obtained  for  isotropic 
materials.  For  example,  for  radially  reinforced  spheres  (i.e.  spheres  transversely 
isotropic  about  the  radial  direction),  it  was  shown  in  [4]  that  a  jumping  phenomenon  can 
occur  which  is  reminiscent  of  snap-through  buckling  in  structural  mechanics.  A 
discontinuous  change  in  stable  equilibrium  configurations  was  encountered  which 
corresponds  to  the  sudden  appearance  of  a  cavity  of  finite  radius  at  a  critical  load.  The 
same  "snap  cavitation"  phenomenon  was  observed  by  Horgan  and  Pence  [20,  21]  in  an 
investigation  of  cavitation  for  composite  incompressible  nonlinearly  elastic  materials 
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with  isotropic  phases. 

The  problems  investigated  in  [4]  and  [20]  required  completely  different  analyses. 
For  the  classes  of  compressible  materials  considered  in  [4],  giving  rise  to  singular 
ordinary  differential  equations,  phase-plane  methods  were  employed  to  yield  qualitative 
information  on  solution  behavior.  On  the  other  hand,  for  incompressible  materials,  the 
zero  volume  change  constraint  allows  for  immediate  determination  of  the  functional  form 
of  the  radial  deformation  field  and  the  hydrostatic  pressure  is  readily  obtained  by  direct 
integration  of  the  equilibrium  equations  (see  e.g.  [20]).  Thus  results  on  cavitation  for 
incompressible  materials  are  generally  obtained  in  much  more  explicit  form  than  is 
possible  for  compressible  materials.  In  fact,  for  the  illustrative  example  of  a  composite 
sphere  composed  of  two  neo-Hookean  materials  perfectly  bonded  across  a  spherical 
interface,  closed  form  analytic  results  were  obtained  in  [20],  [21]. 

The  present  study  was  motivated  by  [4],  [20],  and  [21]  in  several  ways.  First,  it  is 
clearly  of  interest  to  take  advantage  of  the  constraint  of  incompressibility  to  obtain  more 
explicit  information  on  the  effect  of  anisotropy  on  cavitation.  Secondly,  as  was 
remarked  at  the  end  of  [20],  the  "homogenization"  approach  to  composite  material 
mechanics  suggests  it  worthwhile  to  examine  an  "equivalent"  anisotropic  sphere  as  a 
model  for  a  composite  sphere  isotropic  phases.  Finally,  since  consideration  of  a 
composite  isotropic  sphere  [20,  21]  provided  interesting  results  on  cavitation  and 
material  instabilities,  it  is  worth  investigating  the  combined  effects  of  anisotropy  and 
material  inhomogeneity. 

It  is  worth  noting  here,  as  was  also  observed  in  [4],  that  the  radial  anisotropy 
considered  in  [4]  and  in  this  dissertation  does  indeed  arise  in  technological  applications. 
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A  striking  example  occurs  in  casting  of  metals  where  the  temperature  gradient  in  the 
freezing  process  induces  molecular  structure  resulting  in  transverse  isotropy  about  the 
radial  direction  (see  e.g.  Fig.  p.  321  of  [44]). 

Chapter  2  is  concerned  with  material  properties,  and  focuses  mainly  on  the  type  of 
anisotropy  of  concern  in  this  dissertation.  We  consider  the  stored-energy  functions 
corresponding  to  transversely  isotropic  materials  in  the  contexts  of  both  the  linear  and 
nonlinear  theories  of  elasticity,  and  relate  the  infinitesimal  elastic  moduli  to  the  nonlinear 
stored-energy  function.  We  then  obtain  restrictions  (inequalities)  on  the  nonlinear 
stored-energy  consistent  with  those  known  for  the  linearly  elastic  moduli.  Finally,  we 
obtain  further  results  involving  the  nonlinearly  elastic  transversely  isotropic 
incompressible  stored-energy  function  and  then  apply  the  considerations  of  this  Chapter 
to  a  special  class  of  such  materials.  Chapter  3  is  concerned  with  the  cavitation  problem 
for  a  homogeneous  sphere  composed  of  material  which  is  incompressible  and 
transversely  isotropic  about  the  radial  direction,  while  Chapter  4  considers  this  problem 
for  such  a  sphere  which  is  inhomogeneous.  For  simplicity  and  illustration,  the  stability 
analysis  in  Chapter  3  is  carried  out  for  a  specific  material  model.  In  Chapter  4,  cavitation 
is  analyzed  for  a  nonlinear  elastic  composite  sphere  whose  phases  are  composed  of  two 
(different)  completely  general  transversely  isotropic  materials,  and  a  subsequent 
illustrative  example  is  then  considered.  The  general  results  of  Chapter  4  reduce  to  the 
corresponding  general  results  for  the  homogeneous  sphere  with  the  proper  choice  of 
parameters.  1  Finally,  in  Chapter  5  we  summarize  some  of  the  important  and  interesting 
results  of  Chapters  2-4,  discuss  some  additional  aspects  of  the  cavitation  problem,  and 

1  Let  the  Chapter  4  parameters  a,  f  of  (4A.2\,  (4.4.9),  respectively,  be  such  that  f  =  a  =  1,  and/or  set  W2  of  (4.4.1)  identically 
zero  to  recover  the  general  results  for  the  homogeneous  sphere  from  those  of  Chapter  4. 


indicate  some  areas  for  future  research. 


Chapter  2 

ANISOTROPY,  TRANSVERSE  ISOTROPY,  AND  ELASTIC  PARAMETERS 
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2.1.  Introduction 

Material  properties  are  a  basic  aspect  of  continuum  mechanical  modeling. 
Characterization  of  the  types  of  materials  being  considered  is  necessary  to  properly 
formulate  the  mathematical  problem  as  well  as  describe  its  applicability  to  real  materials. 
The  purpose  of  this  Chapter  is  to  give  a  mathematical  as  well  as  physical  interpretation 
for  the  types  of  elastic  materials  of  concern  in  this  dissertation.  In  Section  2.2  we 
describe  what  is  meant  by  isotropic,  anisotropic,  and  transversely  isotropic  materials, 
within  the  context  of  nonlinear  hyperelasticity.  These  concepts  have,  of  course,  been 
formulated  in  the  literature,  and  are  presented  here  as  an  aid  to  the  reader.  In  Section  2.3 
we  consider  the  elastic  parameters  associated  with  infinitesimal  deformations  of  a 
compressible  transversely  isotropic  material,  and  relate  the  five  independent  physical 
moduli  to  parameters  appearing  in  two  different  but  equivalent  representations  for  the 
stored-energy  function.  In  Section  2.4  we  use  the  results  of  Section  2.3  to  recover  the 
elastic  moduli  for  an  incompressible  transversely  isotropic  linearly  elastic  material,  and 
then  relate  these  to  the  nonlinearly  elastic  stored-energy  function  in  Section  2.5.  Finally 
in  that  Section,  we  examine  the  implications  of  these  results  for  a  specific  material  model 
proposed  in  Chapter  3,  Section  5,  and  motivated  by  the  considerations  given  in  that 
Chapter. 
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2.2.  Isotropy,  anisotropy,  and  transverse  isotropy 

We  refer  to  Figure  2.1  and  consider  a  body  B  composed  of  elastic  material  which  in 
its  undeformed  or  reference  configuration  has  associated  with  any  point  of  B,  the 
coordinate  vector  X  =  (X! ,  X2,  X3 ).  A  deformation  of  this  body  is  then  defined  as 

x  =  x(X)  (2.2.1) 

such  that 

J  =  det  Vx  >  0 ,  (2.2.2) 

where  x  =  (x!,x2,  X3)  refers  to  the  position  vector  of  a  point  in  the  deformed 
configuration,  and  the  quantity  J  defined  in  (2.2.2)  is  physically  significant  in  that  it 
represents  (locally)  the  deformed  to  reference  volume  ratio  of  a  material  volume  element. 
The  restriction  J  >  0  in  (2.2.2)  arises  due  to  considerations  of  physically  reasonable 
deformations,  i.e.,  that  material  cannot  be  annihilated  (J  *  0),  and  that  volume  elements 
are  positive.  Associated  with  any  hyperelastic  material  is  a  strain  energy  density  (per 
unit  undeformed  volume)  or  stored-energy  function,  W,  which  is  a  measure  of  the  elastic 
potential  energy  stored  in  the  body  due  to  deformation.  This  scalar  function  W  depends 
on  the  deformation  gradient  tensor,  F,  defined  as 

F=§  (2.2.3) 

with 

J  =  detF  >  0  (2.2.4) 

in  view  of  (2.2.2),  (2.2.3).  As  depicted  in  Fig.  2.1,  if  we  now  rotate  the  deformed  body 
via  Q,  an  arbitrary  proper  orthogonal  tensor,  i.e., 


8 


QQt  =  QtQ  =  I  ,  detQ  =  1 ,  (2.2.5) 

where  I  is  the  3x3  identity  tensor,  then  the  resulting  body  whose  particles  have  the 
position  vector  Qx  has  associated  with  it  a  stored-energy  W(QF).  Then  the  principle  of 
material  frame-indifference,  which  is  also  known  by  several  other  names,  including  the 
principle  of  material  objectivity,  (see  e.g.,  [41])  states  that 

W(F)  =  W(QF) ,  (2.2.6) 

and  thus  a  change  of  observer  does  not  affect  the  elastic  stored  energy.  Equation  (2.2.6) 
is  assumed  to  hold  regardless  of  the  type  of  material  being  considered,  that  is,  regardless 
of  whether  the  elastic  material  is  isotropic  or  anisotropic  (we  define  these  two  terms  more 
precisely  in  what  follows). 


Figure  2.1.  Principle  of  Material  Frame  Indifference:  W(F)  =  W(QF  ) 


Consider  now  the  situation  depicted  in  Figure  2.2: 


undeformed 


rotated 


W(FQ  ) 


deformed 


Figure  2.2.  Isotropy:  W(F)  =  W(FQ  ) 

If  the  undeformed  body  is  first  rotated  via  Q  satisfying  (2.2.5)  so  that  particles  of  the 
rotated  body  have  the  position  vector  QX,  and  is  then  deformed,  the  stored-energy 
associated  with  the  final  deformed  configuration  of  the  body  is  W(FQ).  The  material  is 
isotropic  if  it  further  holds  that 

W(F)  =  W(FQ) ,  (2.2.7) 

so  that  the  orientation  of  an  isotropic  material  cannot  be  detected  from  experiments.  The 
response  of  an  isotropic  material  thus  exhibits  no  preferred  direction,  and  material 
properties  are  the  same  in  all  directions.  If  (2.2.7)  does  not  hold  for  all  Q  satisfying 
(2.2.5),  the  material  is  said  to  be  anisotropic,  and  the  comments  following  (2.2.7)  do  not 
hold  for  such  a  material.  For  anisotropic  materials,  equation  (2.2.7)  is  replaced  by 

W(F)  =  W(FQ),  (2.2.8) 

A  A 

where  Q  satisfies  (2.2.5),  and  the  set  of  all  rotation  tensors  Q  which  satisfy  (2.2.8)  is  a 
proper  subset  of  the  set  of  all  proper  orthogonal  tensors  Q. 

The  type  of  anisotropy  of  concern  in  this  dissertation  is  transverse  isotropy.  A 
transversely  isotropic  material  is  one  which  exhibits  a  single  preferred  direction  at  every 
point,  so  that  material  properties  in  the  direction  of  transverse  isotropy  are  different  than 
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those  in  other  directions.  For  a  material  which  is  transversely  isotropic  about  an  axis,  the 

A 

rotation  tensors  Q  satisfying  (2.2.8)  describe  all  rotations  about  that  axis. 

2.3.  Elastic  parameters:  linear,  compressible,  transversely  isotropic  materials 

The  stored-energy  function,  w,  for  a  compressible  linearly  elastic  material  which  is 
transversely  isotropic  about  an  axis  in  a  Cartesian  coordinate  system  is  given  in  [12].  We 
rewrite  this  quadratic  form  for  a  linearly  elastic  material  that  is  transversely  isotropic 
about  the  1  direction  in  an  orthonormal  basis  with  associated  indices  (1,2,3)  as 

2w  =  a  (E22  +  E33)  +  a  Eij  +  2(a  -  2|i)E22E33  +  2b(E22  +  E33)En 

+  4|l(Ei2  +  E13)  +  4|lE23  7 

where  E  is  the  infinitesimal  strain  tensor  defined  as 

E  =  1/2(F  +  Ft)-I,  (2.3.2) 

the  Ey  (i,j  =  1,2,3)  are  the  components  of  E,  and  a,  a,  b,  |i,  and  jl  are  elastic  parameters 
where  |i  is  the  shear  modulus  associated  with  the  2-3  direction,  and  jl  the  shear  modulus 
associated  with  the  1-2  and  1-3  directions.1  In  [12]  the  relationship  between  the 
parameters  appearing  in  (2.3.1)  and  the  elastic  constants  of  physical  significance  is  also 
given  to  be 

^  _  4p.(aa  -  b2  -  a|i)  -  aa  -  b2  -  a[i 

E- - —  -r -  ,  E=  - 

a  a  -  b  a  -  Li 

(2.3.3) 

aa-b2-aji  -  b 

v  = - t—  V  = - 

A  -  t  2  *  ^rA  .  7 

aa-b2  2(a  -|i) 

_  1  Here  and  in  equation  (2.3.3)  we  depart  slightly  from  the  notation  employed  in  [12]  so  that  the  elastic  moduli  p,  p  and 
E,  E,  v,  v  appearing  in  (2.3.1)  and  (2.3.3),  respectively,  are  such  that  those  with  a  superposed  bar  correspond  to  elastic  moduli 
associated  with  material  properties  in  the  1,  1-2,  or  1-3  directions,  while  those  without  the  superposed  bar  correspond  to  properties  in 
the  2  and/or  3  directions  only.  In  other  words,  ft,  E,  v  indicate  properties  involving  the  direction  of  transverse  isotropy. 
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where  E  and  E  are  the  Young’s  moduli  in  the  1  direction  and  the  2  and  3  directions, 
respectively,  while  v  and  V  are  Poisson’s  ratios  associated  with  the  1-2  and  1-3  directions 
and  the  2-3  direction,  respectively.  In  addition,  the  relation 

E  =  2|i(l  +v)  (2.3.4) 


must  hold,  so  that  there  exist  five  independent  elastic  moduli  associated  with  a 
transversely  isotropic  linearly  elastic  material.  That  is,  we  have  (i,  jl,  E,  E,  v,  and  v 
along  with  (2.3.4).  We  can  invert  the  relations  (2.3.3)  to  give 

E(l-v2E/E) 

(1+v)(1-v-2v2E/E)  ’ 


a  = 


E(l-v)  _ 

(1-v-2v2E/E)  ’ 


(2.3.5) 


b  = 


(1  —  v 


Ev 

-2v2E/E) 


where  (2.3.4)  has  been  used. 

A  general  representation  for  the  strain  energy  function  w  for  a  linear  elastic  solid  as 
a  quadratic  form  in  E  is  given  by  [36]  for  a  material  transversely  isotropic  about  the 
direction  a  as 

2w  =  X(trE)2  +  2(tTtrE2  +  2a  (a*Ea)trE  +  4(|iL  -  (ix)a  *E2a  +  y(a»Ea)2  .  (2.3.6) 

Here  a  is  a  unit  vector  in  the  direction  of  transverse  isotropy  and  X,  (ix,  (1l»  oc,  and  y  are 
elastic  parameters.  Letting  a  =  (1,0,0)  in  (2.3.6),  so  that  the  material  is  transversely 
isotropic  about  the  1  direction,  on  regrouping  the  Ey  as  in  (2.3.1)  we  obtain  the  following 
relationships  among  the  parameters  appearing  in  (2.3.1)  and  (2.3.6): 
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a  =  X  +  2(iT  ,  a  =  X  —  2(ij  +  2 a  +  4|1l  +  Y » 
b  =  X  +  a  ,  (i  =  |iT  ,  P- =  M-l  > 

or,  on  inverting, 

a  =  b  -  a  +2(i  ,  y=a  +  a  -2b-4jl, 

X  =  a  -  2|i  ,  |iT  =  |i  ,  Hl  =  P  • 


(2.3.7) 


(2.3.8) 


Using  (2.3.8),  (2.3.5),  and  (2.3.4),  we  determine  the  parameters  appearing  in  (2.3.6)  in 
terms  of  the  physical  moduli  to  be 


.  _  Efr  +  ^E/E) 

(1+v)(1-v-2v2E/E)  ’ 

^ _  E(v  - v  +  vv - y2E/E) 
(1+v)(1-v-2v2E/E)  ’ 


(2.3.9) 


and  trivially, 


E(1-v2)  +  E(1-2v-2vv-v2E/E)  4_ 

(l+v)(l  -v-2v2E/E) 


Ht  =  1A  ,  4l  =  4- 


(2.3.10) 


Similarly,  (2.3.3)  and  (2.3.7)  give  the  physical  elastic  moduli  associated  with  the 
parameters  of  (2.3.6)  as 


E  =  4|iT  [  1  -  —  ( X  -  2px  +  2a  +  4jil  +  y)  ] , 

E  =  [  D  —  jij(  X  —  2|ij  +  2a  +  4)1l  +  y)  ]  /  (X  +  |ix ) , 

2\ij  ^ 

\  =  l--^-(X-2\iT  +  2a  +  4\iL  +  y) , 

\  =  (X  +  a)/2(X  +  \LT) , 


(2.3.11) 


where  D  in  (2.3. 1 1)  is  defined  as 
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D  =  (X  +  2p.T)(X,-2|XT  +  2a  +  4(iL  +  y)-(X  +  d)2 .  (2.3.12) 


Finally,  from  [12]  we  also  have  the  following  necessary  and  sufficient  conditions  for 
positive  definiteness  of  w: 


E>0  ,  E>0  ,  |i>0  ,  (i>0  , 

— 1  <  v  <  1  ,  l-v>2  -1-v2  . 

E 


(2.3.13) 


2.4.  Elastic  parameters:  linear,  incompressible,  transversely  isotropic  materials 

We  first  wish  to  show  how  the  constitutive  law  for  an  incompressible  linearly  elastic 
transversely  isotropic  material  can  be  recovered  from  that  for  a  compressible  material. 
From  [35]  we  write  down  the  latter  for  a  material  with  preferred  direction  along  the  unit 
vector  a  as 

Oy  =  XEjdfSy  +  2jixEjj  +  a(  atamEfonfiy  +  E^a^ )  (2  41) 

+  2(|1l  _  M-T)(aiakEkj  +  aja^E^)  +  YakamEkmaiaj  > 

where  i,j=  1,2, 3,  Cy  are  components  of  the  stress  tensor  in  the  linear  theory,  and 
summation  of  repeated  indices  is  implied.  As  indicated  in  [35],  the  idealized 
incompressible  material  may  be  regarded  as  the  limiting  case  of  the  material  with 
constitutive  equation  (2.4.1)  when 

X  — >  °°  and  Ekk  — >  0  such  that  XE^  — >  -q ,  (2.4.2) 

where  q  is  an  arbitrary  finite  constant.  Applying  (2.4.2)  to  (2.4.1),  we  obtain  the 
constitutive  equation  for  an  incompressible  linearly  elastic  transversely  isotropic  material 
with  preferred  direction  a  as 
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°ij  -  ~  ( q  ~  «akamEkm)8ij  +  2jJ.TEy  +  2(jJ.|_  -  |lT)(aiakEkj  +  aj akEki)  ^  4  3) 

"b  y akamEkmaiaj  • 

As  q  is  arbitrary,  the  term  in  (2.4.3)  multiplying  Sy  is  often  written  simply  as  -p,  p 
arbitrary. 

We  now  determine  the  corresponding  stored-energy  w  from  (2.4.3),  (2.4.2),  and 

2w  =  OyEy  (2.4.4) 

to  be 

2w  =  2|iTEjjEy  +  2(Pl  ~  M,T)(aiakEkj  4"  ajakEki)Eij  "b  YakamEkmaiajEij  »  (2.4.5) 

or  since  E  is  symmetric, 

2w  =  2|iTEijEjj  +  4((J.l  ~  pT)ajakEkiEij  +  y (ajajEij)2  .  (2.4.6) 

In  direct  tensor  notation,  (2.4.6)  is  equivalent  to 

2w  =  2(iTtrE2  +  4(|iL  -  |xT)(a  •  E2  a)  +  y  (a  •  Ea)2  .  (2.4.7) 

Letting  a  =  (1,0,0),  the  stored-energy  for  an  incompressible  linearly  elastic  material  that 
is  transversely  isotropic  about  the  1  direction  has  the  form 

2w  =  2|iT[ 2E22  +  2E?3  +  2E|>  -  (E?i  +  E^2  +  E33)  1 

7  v,2  r2N  r-,2  (2.4.8) 

+  4pL(Eu  +  E12  +  E^)  "b  YE 1 1  • 

We  can  determine  the  corresponding  physical  elastic  moduli  from  the  results  of  the 
previous  section  and  (2.4.2).  From  (2.3.10),  we  again  have 

4  =  4t  ,  4  =  14.  (2-4.9) 

as  the  shear  moduli,  while  the  Young’s  moduli  and  Poisson’s  ratios  can  be  determined 
from  (2.3.11),  (2.3.12)  in  the  limit  as  X  -»  «».  Taking  these  limits  and  using  l’Hopital’s 
rule,  the  elastic  parameters  corresponding  to  incompressibility  can  be  given  as 
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lim  E  =  4(iT[  1  - 
X— 


Mr 

4^l+Y 
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lim  E  =  4(iL  +  y-  |iT  , 

4— >oo 


lim  v  =  1  - 
X— 


4Hl+Y  ’ 


lim  v  =  lh , 

A,-*®* 


> 


(2.4.10) 


where  again  in  (2.4.9),  (2.4.10),  the  superposed  bar  refers  to  material  properties  involving 
the  direction  of  transverse  isotropy  (e.g.,  the  1  direction)  while  parameters  on  the  left 
hand  sides  of  these  equations  without  the  superposed  bar  describe  properties  independent 
of  this  direction  (e.g.  the  2  and/or  3  directions  only). 


2.5.  Implications  for  the  nonlinearly  elastic  stored-energy 

We  return  now  to  nonlinear  hyperelasticity.  As  discussed  in  Section  2.2,  there 
exists  a  stored-energy  function 

W  =  W(F)  (2.5.1) 

associated  with  any  nonlinearly  elastic  material.  For  W  to  satisfiy  (2.2.6)  for  all  proper 
orthogonal  Q,  it  can  be  shown  further  that  W  may  be  regarded  as  a  function  of  C  =  FTF, 
i.e., 

W  =  W(C).  (2.5.2) 

Thus  the  strain  energy  W  is  independent  of  superposed  rigid  body  motions  (i.e.,  satisfies 
frame-indifference)  when  (2.5.2)  holds  instead  of  (2.5.1).  (In  the  above  and  in  what 
follows  we  only  introduce  new  symbols  for  W  depending  on  different  arguments  when 
there  is  a  need  to  avoid  confusion.  The  functions  W  on  the  right  hand  sides  of  (2.5.1)  and 
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(2.5.2)  are  clearly  not  the  same,  but  no  new  notation  is  employed  in  order  to  avoid  an 
unnecessary  proliferation  of  symbols.)  In  the  case  of  a  transversely  isotropic  material, 
equation  (2.5.2)  can  be  further  simplified  so  that  W  depends  on  C  only  through  five 
invariants.  Thus,  (see  e.g.  [27]), 


W  =  W(I1,I2,I3,I4,I5), 


(2.5.3) 


where  without  loss  of  generality,  taking  the  axis  of  symmetry  (direction  of  transverse 
isotropy)  to  be  the  1  direction,  we  can  define  the  invariants  in  (2.5.3)  to  be 


Ii=trC, 

12  =  1/2[(trC)2-trC2], 

13  =  detC , 

14  =  C?2  +  C?3, 

15  —  C„  . 


(2.5.4) 


(We  note  that  for  an  isotropic  material,  the  stored-energy  (2.5.3)  depends  only  on  the  first 
three  invariants  defined  exactly  as  in  (2.5.4).)  Recalling  from  Section  2.2  that  J  =  detF  is 
a  measure  of  volume  change,  the  constraint 

J=1  =>  I3=J2  =  1  (2.5.5) 

must  hold  when  the  material  is  incompressible  so  that  volume  is  preserved.  Thus,  for  an 
incompressible  material,  (2.5.3)  becomes 

W  =  Wd1,I2,I4,I5)  (2.5.6) 

and  (2.5.5)  must  hold.  With  (2.5.6)  as  the  general  form  for  the  stored-energy  function  of 
a  transversely  isotropic  incompressible  material,  our  next  immediate  goal  is  to  determine 
a  generalized  Young’s  modulus  in  the  direction  of  transverse  isotropy  in  order  to  relate 
this  nonlinear  stored-energy  to  the  quantity  E  from  the  linear  theory  discussed  in  Section 


2.4.  The  importance  of  such  a  relationship,  aside  from  its  intrinsic  interest,  will  become 
evident  in  Chapters  3  and  4  of  this  dissertation. 

2.5.1  Generalized  moduli  and  restrictions  on  W(IX,  I2, 14, 15) 

(i)  Uniaxial  tension: 

Consider  first  a  uniaxial  tension  deformation  in  the  Xj  direction  for  a  nonlinearly 
elastic  transversely  isotropic  material  with  stored-energy  described  by  (2.5.4)  -  (2.5.6). 
This  deformation  is  examined  in  [27,  p.473].  As  in  [27],  we  write  (2.2.1)  for  this 
deformation  as 

xx  =  AXi  ,  x2=X2/^A,  x3  =  X3/ >/a  ,  (2.5.7) 

so  that  corresponding  to  (2.5.7),  the  deformation  gradient  tensor  F  defined  in  (2.2.3)  is 

A  0  0 

F  =  0  1/  yA  0  .  (2.5.8) 

0  0  i/Va 

Then 

Ta2  0  0 

C  =  FtF=  0  1/A  0  (2.5.9) 

0  0  1/A 

and 

I3  =  detC  =  1 ,  (2.5.10) 

so  that  the  incompressibility  constraint  (2.5.5)  holds.  It  is  shown  in  [27]  that  the  only 
non-zero  stress  Tn  can  be  written  as 

Tn  =  2 A  [  (A  -  -y )  C W  j  +  -f W2 )  +  A' W5  ] , 

A1  A 


(2.5.11) 
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where  the  Wk  =  5W/  3lk  are  evaluated  at  the  following  values  of  the  invariants  (defined 
by  (2.5.4))  which  correspond  to  this  deformation: 

11  =  A2  +2/ A , 

12  =  2A  +  A-2  , 

I4=0, 

Is  =  A2  . 

We  will  return  to  (2.5.11),  (2.5.12)  after  the  following  discussion  of  the  uniaxial  tension 
deformation  for  incompressible  transversely  isotropic  linearly  elastic  materials. 


Consider  now  uniaxial  tension  in  the  context  of  the  linear  theory  of  elasticity.  From 
(2.3.2)  and  (2.5.8),  the  infinitesimal  strain  tensor  E  has  the  form 

A  - 1  0  0 

■  E=  0  1/VA-l  0  .  (2.5.13) 

0  0  l/^A-l 

We  write  the  linear  stress-strain  equations  involving  the  physical  moduli  as  (see  e.g., 

[12]) 


En  =  ^[O’ll  —  V(<T22  +^33)]  > 

E 

E22  =  Trta22  -VCT33  -  -i-VOn]  , 

E  E 

c  lr  -  E-  , 
E33  =  — [033  -va22  -  — von] , 
E  E 

°12  =  a21  =  2pE!2  , 

^13  =  <*31  =  2jIEi3  , 


(2.5.14) 


ct23  “  °32  =  2(lE23  . 
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The  equation  of  interest  in  (2.5.14)  involving  Young’s  modulus  in  the  1  direction  (the 
direction  of  transverse  isotropy)  is  (2.5.14)i.  Now  On  is  the  only  non-zero  stress 
component  and  so  (2.5.13)  and  (2.5.14)!  give 

an=E(A-l).  (2.5.15) 

Thus,  if  we  can  write  the  nonlinear  stress  Tn  of  (2.5.11)  in  the  same  form  as  the  linear 
stress  an  of  (2.5.15),  i.e, 

TU=E(A)(A-1),  (2.5.16) 

then  we  call  E(A)  the  generalized  Young's  modulus  in  the  direction  of  transverse 
isotropy.  This  is  analogous  to  the  concept  of  the  generalized  shear  modulus  determined 
from  comparisons  of  the  linear  and  nonlinear  (isotropic)  simple  shear  problems  (see  e.g., 
[27],  pp.  447-453  or  [15],  p.176). 

To  arrive  at  (2.5.16),  we  first  multiply  the  right  hand  side  of  (2.5.11)  by  A2/ A2: 

Tn  =  \  [  (A3  -  1)  ( Wi  +  -f  W2 )  +  A3  W5  ] 

A  A 

=  ^-[(A-1)(A2  +  A+1)(Wi  +  ^W2)+A3W5], 

A  A 

or 

2 

Tn  =  (A  -  1)  [  \  (A2  +  A  +  1)  ( Wi  +  -f  W2 )  +  W5  ] .  (2.5.17) 

A  A  (A  —  1) 

Equation  (2.5.17)  is  of  the  form  (2.5.16)  where 

E(A)  =  \  (A2  +  A  +  1)  C Wi  +  -J-W2 )  +  77^rW5  (2.5.18) 

is  the  generalized  Young’s  modulus  in  the  direction  of  transverse  isotropy.  To  be 
consistent  with  the  classical  linear  theory  (see  e.g.,  [41]  or  [29],  p.349),  equation  (2.5.18) 
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should  be  such  that 


lim  E(A)  =  E , 

A-»  1 


(2.5.19) 


where  E  on  the  right  in  (2.5.19)  is  the  E  of  Section  2.4,  and  the  Wk  of  (2.5.18)  are 
evaluated  at  the  values  of  the  invariants  (2.5.12)  consistent  with  A  — >  1: 
It  =  I2  =  3, 14  =  0, 15  =  1.  Applying  (2.5.19)  to  (2.5.18)  we  obtain 


A2W5 

E  =  6(W1+W2)Iii=i2  =  3  +2  lim  — — f 

A — ^  1  A.  1 

14=0,15  =  1 


(2.5.20) 


The  remaining  limit  in  (2.5.20)  is  an  indeterminate  form,  since  it  will  be  seen  in  Chapter 
3  that 


W5(3,  3, 1,0, 1)  =  0 


(2.5.21) 


(see  equation  (3.2.10)  and  the  sentence  preceding  it).  On  using  l’Hopital’s  rule,  we  find 


A2W5  ,  dW5 

lim  — — f  =  lim  (A2———  +2AW5)Ii1=i2  =  3 
A  — ►  1  A  “  1  A  — ►  1  d  A 

14  =  0, i5  =  i 

dW5 

=  *lim,  ~7Z~ llj  =I2=3 
a  — ►  1  d  A 

14  =  0,  i5  =  i 

on  using  (2.5.21).  Now,  since  I3  s  1, 14  =  0,  the  chain  rule  gives 


(2.5.22) 


dw5  aw5  dix  aw5  di2  aw5  di5 

dA  =  aix  dA  +  dl2  dA  +  3l5  dA  ’ 


where,  from  (2.5.12), 


(2.5.23) 


dlj  ,  dl2  ,  dl5 

-f  =  2A-2A"2  ,  -f  =  2-2A"3  ,  -f  =  2A 
dA  dA  dA 


(2.5.24) 


Substituting  (2.5.24)  into  (2.5.23)  and  then  taking  the  limit  of  (2.5.23)  as  A— » 1,  we 


obtain 
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lim 
A— ►  1 


A2W5  dW5 

- —  =  lim  — — 

A  —  1  A-»  l  dA 


=  2W55  I  i,  =  I2  =  3  , 

14  =  0,  i5  =  i 


(2.5.25) 


where  W55  =  dWs/Bls  =  92W/9l2.  Substituting  (2.5.25)  into  (2.5.20),  we  determine  the 
relation  between  the  linear  Young’s  modulus  E  and  the  nonlinear  strain  energy  W  to  be 


E  =  [ 6(W  1  +  W2)  +  4W55]  lij  =i2  =  3  .  (2.5.26) 

14  =  0,15  =  1 

Then  (2.5.26)  and  (2.3.13)2  imply  the  following  condition  on  the  nonlinear  stored-energy 
W: 


[3(W1+W2)  +  2W55]  Iii=i2=3  >0.  (2.5.27) 

i4  =  o,  i5  =  1 

(ii)  Constitutive  inequalities: 

The  inequality  (2.5.27)  results  from  requiring  consistency  between  the  linear  and 
nonlinear  theories,  and  from  one  of  the  conditions  for  positive  definiteness  of  the  linear 
strain  energy  density.  There  are  however  six  such  conditions  for  positive  definite  w 
given  in  (2.3.13),  which  we  restate  here  for  convenience: 


E >0  , 

(2.5.28) 

E >0  , 

(2.5.29) 

|X  >  0  , 

(2.5.30) 

pi  >  0  , 

(2.5.31) 

— 1  <  v  <  1  , 

(2.5.32) 

-E_2 

-v > 2— v  . 

E 

(2.5.33) 

Constitutive  inequalities  are  of  great  importance  in  elasticity  to  ensure  physically 
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reasonable  response.  A  thorough  examination  of  such  constitutive  restrictions  is  carried 
out  in  [41]  for  compressible  nonlinearly  elastic  materials.  As  remarked  in  [41],  no 
comparable  fundamental  studies  of  inequalities  had  been  undertaken  for  incompressible 
materials,  which  to  this  author’s  knowledge,  remains  the  case  today.  It  is  therefore  of 
interest  to  determine  what  restrictions  are  placed  on  the  nonlinear  stored-energy  W  from 
all  of  (2.5.28)  -  (2.5.33).  We  have  thus  far  considered  only  (2.5.29).  Toward  this  end, 
we  first  wish  to  show  that  in  the  case  of  incompressibility,  it  is  not  necessary  to  consider 
all  of  the  conditions  (2.5.28)  -  (2.5.33)  from  the  linear  theory.  Since  we  have  already 
determined  the  necessary  restriction  placed  on  the  nonlinear  strain  energy  W  by  (2.5.29), 
we  assume  in  what  follows  that  E>0  and  show  that  the  above  conditions  are  not 
completely  independent.  We  will  show  that  it  suffices  to  consider  only  (2.5.29)  - 
(2.5.31),  and  thus  we  need  only  be  concerned  with  two  additional  conditions  analogous 
to  (2.5.27).  These  two  conditions  on  W  associated  with  (2.5.30)  and  (2.5.31)  will  be 
determined  later  in  this  Section. 


Proposition  1:  Consider  an  incompressible  linearly  elastic  transversely  isotropic  material. 
Suppose  E  >  0.  Then  E  >  0  if  p  >  0. 


Proof:  From  (2.4.10)2,  (2.4.9)i, 


E>0  =>  4pL+y-px>0 

p  >  0  =>  px  >  0 


Dividing  (2.5.34)  by  px,  we  have 


4pL+Y 

px 


-1>0 


(2.5.34) 

(2.5.35) 
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4^l+Y  ,  Mr  , 

=>  - >  1  =>  — - <  1 

Mr  4|iL+y 

=>  1  -  - - >0 . 

4^l+Y 

Thus  from  (2.4.10)!,  (2.5.35),  and  (2.5.36), 


E  =  4(iT[  1 


Mr 


4|xl  +  Y 


]>0, 


and  so  (2.5.28)  holds  whenever  (2.5.29)  and  (2.5.30)  hold. 


(2.5.36) 


(2.5.37) 


Proposition  2:  Consider  an  incompressible  linearly  elastic  transversely  isotropic  material. 
Suppose  E  >  0.  Then  -1  <  v  <  1  if  |X  >  0. 


Proof:  We  first  show  the  left  hand  inequality,  -1  <v.  As  in  Proposition  1,  (2.5.34), 
(2.5.35)  lead  to  (2.5.36).  Multiplying  (2.5.36)  by  two,  we  have 


or 


4pL+Y 


>0 


(2.5.38) 


4Ml  +  Y 


>-l . 


(2.5.39) 


Thus  from  (2.4.10)3 » we  have  the  desired  result  v  >  -1. 

We  now  show  the  right  hand  inequality  of  (2.5.32).  From  (2.5.34)  and  (2.5.35), 


and  so 


4(4.  +  Y  >  Pt  >  0 


(2.5.40) 


0< 


M-t 

4liL  +  Y  ' 


(2.5.41) 


Thus  from  (2.5.41)  and  (2.4.10)3, 
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V  =  1  — 


2Pt  1 

1 - —  <*» 

4^l  +  Y 


(2.5.42) 


and  (2.5.32)  therefore  holds  whenever  (2.5.29),  (2.5.30)  hold.  This  completes  the  proof 
of  Proposition  2. 

Consider  now  the  condition  (2.5.33).  From  (2.4.10)3,  the  left  hand  side  of  (2.5.33) 
is 


1  -  v  = 


2|xt 


4|4L  +  Y  ’ 


(2.5.43) 


while  the  right  hand  side  is,  from  (2.4.10), 


_  E-2  _ 

2— v  =2- 


4(J-t 


4Pl  +  Y-Pt 


4pL  +  Y 


4Pl  +  Y-Pt 


_2 

V 


8Mt  -2 

=  - v  . 

4|XL+Y 

Therefore,  since  v  — >  Vi  according  to  (2.4.10)4  in  the  incompressibile  case,  we  have 

E  _2  2|!t  _ 

2— v  - -  as  v  V4 .  (2.5.44) 

E  4jil  +  y 

While  equations  (2.5.43),  (2.5.44)  may  at  first  seem  to  contradict  (2.5.33),  this  situation  is 
analogous  however  to  that  encountered  in  isotropic  elasticity.  It  is  well  known  in  linear 
elasticity  that  the  conditions  which  guarantee  a  positive  definite  isotropic  stored-energy 
are  E  >  0 ,  -1  <  v  <  Vi,  but  that  v  =  Vi  corresponds  to  incompressibility  (see  e.g.,  [15]). 
Similarly,  for  transverse  isotropy  we  have  in  general  (2.5.33)  for  positive  definite  w,  but 
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*  „  E-2  E 

1  -  v  =  2— v  =—  (2.5.45) 

E  2E 

corresponds  to  incompressibility,  where  (2.4.10)4  has  been  employed  to  obtain  the  final 

equality  in  (2.5.45). 

We  thus  have  that  in  infinitesimal  deformations  of  an  incompressible  transversely 
isotropic  material,  the  condition  (2.5.33)  is  replaced  by  (2.5.45),  and  the  inequalities 
(2.5.28),  (2.5.32)  hold  whenever  (2.5.29),  (2.5.30)  hold.  As  stated  earlier,  the  inequality 
(2.5.27)  follows  from  (2.5.29).  It  therefore  remains  to  determine  what  conditions  are 
placed  on  the  nonlinear  stored-energy  by  (2.5.30),  (2.5.31),  and  (2.5.45)  only,  to 
completely  account  for  the  conditions  (2.5.28)  -  (2.5.33)  from  the  linear  theory  of 
elasticity.  We  next  address  (2.5.30)  and  (2.5.31)  in  a  manner  similar  to  that  undertaken 
in  part  (i)  of  Section  2.5.1.  That  is,  we  consider  two  simple  shear  deformations  in  order 
to  determine  generalized  shear  moduli,  p(K2)  and  |i(k2),  that  can  then  be  related  to 
jl  and  p  of  the  linear  theory  (see  below).  Finally  in  Section  2.5.1,  we  show  that  (2.5.45) 
puts  no  further  restrictions  on  W. 


and 


1  k  0 
F=  0  1  0 
0  0  1 


(2.5.47) 
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I3  =  detC  =  1  (2.5.48) 

so  that  the  incompressibility  constraint  (2.5.5)  holds.  It  can  then  be  shown  that  the  only 
non-zero  shearing  stress  Ti2  =  T2i  has  the  form 

T12  =  2k  (Wi  +  W2  +  W4 ) ,  (2.5.49) 

where  the  Wk  s  3W/  3lk  are  evaluated  at  the  following  values  of  the  invariants  (2.5.4) 
corresponding  to  this  deformation: 

I1=I2=3  +  k2, 

14=*?,  -  (2.5.50) 

Is  = 1  • 

We  now  consider  (2.5.46)  in  the  context  of  linear  theory.  From  (2.3.2)  and 
(2.5.47)i ,  the  infinitesimal  strain  tensor  has  the  form 

0  k/2  O' 

E  =  k/2  0  0  ,  (2.5.51) 

_  0  0  0. 

and  Oi2,  the  only  non-zero  stress  component,  is  thus  given  by  (2.5.14)4  and  (2.5.51)  as 

oi2  =  pK .  (2.5.52) 

On  comparing  (2.5.52)  and  (2.5.49),  we  define  the  generalized  shear  modulus  in  the  1-2 
direction  as 

P(k2)  =  2(Wx  +  W2  +  W4) ,  (2.5.53) 

where  p  as  a  function  of  k2  results  from  the  Wk  being  evaluated  at  the  values  of  the 
invariants  given  by  (2.5.50).  Then  in  the  limit  as  k2  -»  0,  we  recover  the  linear  theory, 
so  that  the  relationship  between  p  and  the  nonlinear  stored-energy  W  is 
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(I  =  lim  p(K2)  =  2(Wi  +  W2  +  W4)  1 1!  =  i2  =  3  .  (2.5.54) 

I4 = o ,  Ij  =  1 

Therefore,  the  inequality  (2.5.31)  implies  the  condition 

(Wi+W2+W4)Ii1=i2  =  3  >0.  (2.5.55) 

14  =  0, i5  =  i 

We  remark  that  (2.5.54),  (2.5.55)  could  have  also  been  obtained  by  considering  simple 
shear  in  the  1-3  direction,  so  that  starting  with  xx  =Xx  +kX3  ,  x2  =X2  ,  x3  =X3  , 
instead  of  (2.5.46)  would  have  resulted  in  (2.5.54)  on  comparison  of  T13  and  a13.  Thus, 
(2.5.53)  actually  defines  the  generalized  shear  modulus  involving  the  direction  of 
transverse  isotropy,  i.e.,  in  both  the  1-2  and  1-3  directions. 

We  proceed  now  to  determine  the  generalized  shear  modulus  |J.(k2)  in  the  2-3 
direction,  and  thus  the  condition  on  the  strain  energy  W  implied  by  (2.5.30),  by 
examining  a  second  simple  shear  deformation  defined  by 

xj  =Xj  ,  x2  =X2  +  kX3  ,  x3  =X3  .  (2.5.56) 

Then  F,  C  corresponding  to  (2.5.56)  are 

100]  100 

F=  0  1k  ,C=01  k  ,  (2.5.57) 

.0  0  lj  0  k  1+k2 

and  again,  incompressibility  holds  since  I3  =  1.  We  then  have  the  only  non- zero  shearing 
stress  T23  =  T32  as 

T23  =2k(Wi  +W2) ,  (2.5.58) 

where  the  W*  are  evaluated  at  the  values  of  the  invariants 


Il  —  I2  -  3  +  K2  , 
14=0, 

I5  =  l, 
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(2.5.59) 

associated  with  (2.5.57)2. 

In  the  linear  theory,  the  only  non-zero  stress  component  a23  is  obtained  from 
(2.5.14)6  and 

'0  0  O' 

E  =  0  0  k/2  ,  (2.5.60) 

.0  k/2  0 

as 

a23=ltK.  (2.5.61) 

Thus  from  (2.5.58)  and  (2.5.61),  the  generalized  shear  modulus  in  the  2-3  direction  is 

HOc2)  =  2(Wj  +  W2)  I  ii  =  h  =  3  +  k2  ,  (2.5.62) 

i4  =  o.i5  =  i 

so  that  the  relationship  between  p  and  the  nonlinear  stored-energy  W  is 

p=  lira  P(k2)  =  2(W1  +W2)Ii1=i2=3  .  (2.5.63) 

I4 =  0 , 15  =  1 

Thus  the  inequality  (2.5.30)  implies  the  following  condition  on  W: 

(Wi +W2)Ii1=i2  =  3  >0.  (2.5.64) 

14  =  0. i5  =  i 

We  remark  that  (2.5.62)  -  (2.5.64)  reduce  to  the  corresponding  results  for  isotropic 

materials.  That  is,  when  the  nonlinear  stored-energy  is  viewed  as  a  function  of  Ii  and  I2 

only,  and  p  represents  the  shear  modulus  in  infinitesimal  deformations  of  an  isotropic 

material,  p  =  2(W  1  +  W2)  I  ij  =  i2  =  3  >  0,  where  the  inequality  is  also  a  consequence  of  the 

a3  =  D 

well-known  Baker-Ericksen  inequalities.  The  above  is  not  surprising  in  that  a  material 
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that  is  transversely  isotropic  about  the  1  direction  may  in  some  sense  be  viewed  as 
"isotropic"  when  only  the  2-3  direction  is  involved. 

We  have  therefore  determined  three  conditions,  (2.5.27),  (2.5.55),  and  (2.5.64)  on 
the  nonlinear  stored-energy  W,  which  result  from  (2.5.28)  -  (2.5.33).  We  claim  that  these 
three  conditions  are  sufficient  to  guarantee  that  all  of  (2.5.28)  -  (2.5.32),  and  (2.5.45)  are 
satisfied.  For  this  claim  to  hold,  it  only  remains  to  show  that  (2.5.45)  puts  no  further 
restrictions  on  W.  We  do  this  in  the  upcoming  final  part  of  Section  2.5.1  by  considering 
an  equibiaxial  stress/stretch  deformation. 

(iv)  Equibiaxial  stress/stretch: 

We  finally  consider  an  equibiaxial  stress/stretch  deformation  defined  by 

xj  =  r-2xt ,  x2  =  rx2 ,  x3  =  rx3 ,  (2.5.65) 

with  corresponding  deformation  tensors 

r-2  oo  r-4  o  o 

F=  0  r  0  ,  C=  0  r2  0  .  (2.5.66) 

o  o  r  o  o  r2 

We  again  have  I3  =  1  so  that  the  deformation  is  volume  preserving.  It  can  then  be  shown 
that  on  requiring  Tn  =  0,  the  only  remaining  non-zero  stresses  are 

t = t22 = t33  =  2[  (r2  -  r-4 )  Wi  +  (r4  -  r~2  )w2  -  r-4w5  ] ,  (2.5.67) 


where  the  Wk  are  evaluated  at 
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i1=r-4  +  2r2, 
i2=r4  +  2r-2, 
i4=o, 
i5=r-4. 


(2.5.68) 


For  reasons  that  will  become  clear  in  the  next  paragraph,  we  rewrite  (2.5.67)  as 

Tn  =  pr  [  (r6  - !)' W1  +  (r8  -  r2)  w2  -  w5  ] 

=  ^-[(r3  -  i)(r3  +  i)Wi  +  r2(r3  - 1  xr3  +  i)w2  -w5] 

=  ^-[(r-i)(r2  +  r+i)(r3  +  i)(Wi+r2w2)  -w5], 

or 

„  ,  Wi  w2  w5 

T  =  2[(r2+r+l)(r3  +  l)(-pf  +  -j5-)-  r4  ](T-1).  (2.5.69) 


Consider  now  the  equibiaxial  stress/stretch  deformation  in  the  linear  theory.  Then 
E  corresponding  to  (2.5.66)  i  is 


E  = 


0 

r-i 

o 


o 

o 

r- i 


(2.5.70) 


so  that  from  (2.5.70),  (2.5.14)2,  (2.5.14)3,  an  =  0,  and  a22  =  a33  =  a,  we  have 


(r-i)=4-a-v)o 

E 


(2.5.71) 


or 


a=7ir^(r-1)-  (2.5.72) 

On  comparing  (2.5.69)  and  (2.5.72),  we  define  the  following  generalized  parameter 
E 


1  — v 


which  is  a  function  of  T  as 
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(1  -v) 


(T)  =  2[(r2+r+i)(P  +  i)( 


Wx  w2 


)- 


w5 


r4(r-i) 


],  (2.5.73) 


where  the  Wk  are  evaluated  at  the  values  of  the  invariants  given  by  (2.5.68).  We  now 
recall  the  condition  (2.5.45),  which  we  write  here  as 


2E  = 


(1-v) 


(2.5.74) 


Now  our  aim  is  to  show  that  (2.5.45),  or  equivalently,  (2.5.74),  puts  no  additional 
restrictions  on  the  nonlinear  stored-energy  W.  This  will  be  accomplished  by  showing 
that  the  generalized  parameters  (2.5.18)  and  (2.5.73)  satisfy 


2E(1)  =  .  (2.5.75) 

where  E(l)  is  given  by  (2.5.19),  (2.5.26).  Taking  the  limit  of  (2.5.73)  as  T  — »  1,  we 
obtain 


(1 


(1)=  lim 
■v)  r-»i 


(1-v) 


(r)  =  2[6(Wi  +W2)-  lim 


W5 


r->i  T4(T-1) 


]  111  =I2  =  3 
I4=0,I5  =  1 


(2.5.76) 


where,  similarly  to  what  occurred  in  part  (i)  of  this  Section,  the  remaining  limit  in 


(2.5.76)  is  as  yet  indeterminate  by  virtue  of  (2.5.21).  On  using  l’Hopital’s  rule,  we  have 


that 


lim 

r->i 


W5  dW5/dr 

— - =  lim  - r - 

T4(T-1)  r->i  T4  +  4T3(T-  1) 


dW5 
"d F 


I  ii  =  i2  =  3 
14  =  0,15  =  1 


By  the  chain  rule 


(2.5.77) 


dW5  dW5  dlk 
dr  “  aik  dr 


(2.5.78) 


where  ,  since  I3  =  1 , 14  =  0,  summation  over  the  indices  k=  1, 2,  5  is  implied.  From 


32 


(2.5.68), 

—  =-4r~5  +4r  ,  =  4r3  -  4r3  ,  =  -4I"5  ,  (2.5.79) 

dr  dr  dr 

so  that  on  substituting  (2.5.79)  into  (2.5.78)  and  letting  r  — >  1,  we  obtain 

ip-H-W  =-4W55Ii,=ii=3  .  (2.5-80) 

I4  =  0 , 15  =  1  14=0,15  =  1 

Then  (2.5.80),  (2.5.77),  and  (2.5.76)  give 

— —  (1)  =  2[6(Wi+W2)  +  4W55]Ii1=i2=3  ,  (2.5.81) 

(!-v)  14  =  0,15  =  1 

which  is  exactly  2E(1)  where  E(l)  is  given  by  (2.5.26).  Thus  the  condition  (2.5.75)  holds 

identically  for  the  nonlinear  stored- energy,  and  so  (2.5.74),  or  equivalently  (2.5.45),  puts 

no  further  restrictions  on  W. 

2.5.2  Further  considerations  of  W(Ii,  I2, 1, 0, 15) 

We  wish  now  to  show  the  connection  between  the  preceding  and  the  cavitation 
problem  of  concern  in  this  dissertation.  We  make  no  attempt  at  this  point  to  elaborate  on 
this  problem  or  on  the  determination  of  the  invariants  or  deformation  tensors  associated 
with  it.  This  will  be  addressed  in  the  next  chapter.  The  results  to  come  in  this  Section 
depend  only  on  the  assumptions  of  incompressibility  and  transverse  isotropy,  and  on  the 
ability  to  write  the  invariants  associated  with  a  particular  deformation  in  the  following 


form: 
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I1(v)  =  v"4  +  2v2, 
I2(v)  =  v4  +  2v-2, 

13  =  l. 

14  =  0, 

I5(v)  =  v-4 . 


(2.5.82) 


Thus,  Section  2.5.2  may  be  applicable  to  a  broad  array  of  problems  which  includes  the 
cavitation  problem.  For  example,  the  invariants  (2.5.68)  associated  with  equibiaxial 
stress/stretch  are  of  this  form,  and  on  defining  A  =  v-2  in  the  uniaxial  tension 
deformation  of  Section  2.5.1(i),  the  invariants  (2.5.12)  become  (2.5.82).  For  this  reason 
and  because  the  goals  of  this  Section  are  in  the  spirit  of  this  Chapter,  we  present  this 
material  here  rather  than  elsewhere. 


We  introduce  the  notation 


W  at,  h,  1,  0, 15)  =  W  at(v),  I2(v),  I5(v))  =  W(v) , 


(2.5.83) 


and  we  first  intend  to  show  that 


d2W 

dv2 


(1)  =  4E(>0) . 


(2.5.84) 


With  this  goal  in  mind,  we  define 


W !  (v)  = 


dW  =  8W  dlk 
dv  Blk  dv 


(2.5.85) 


where  summation  over  the  indices  k  =  1, 2, 5  is  assumed.  From  (2.5.82), 


=-4v-5  +4v  ,  =  4v3  -  4v~3  ,  -^-  =  -4v~5 

dv  dv  dv 


(2.5.86) 


dli  dl2 

When  evaluated  at  v  =  1, - =  — —  =  0,  so  that 

dv  dv 
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*  dW  3W 

Wi  (1)  =  -—  (1)  =  -  4  I  =  -  4WS  (3,3, 1,0, 1)  =  0 

dv  dls  v=i 


on  using  (2.5.21).  Now,  from  (2.5.85) 

d2W  ,  _  dW  d2Ik  d 
dv2  dlk  dv2  dv 

so  that,  by  the  arguments  leading  to  (2.5.87), 


O' 


r 


dw 


dik 

V. 


dl* 

dv 


dv 


dw  d2Ii  aw  d2I2  (  d 

dw 

dl5' 

3li  dv2  51 2  dv2  dv 

di5 

L  J 

dv 

V=1 


From  (2.5.86), 


d2Ii  .  d2I2  ,  , 

— =  20v-6  +  4  ,  — f  =  12v2  +  12v-4 
dv2  dv2 


On  evaluating  (2.5.90),  (2.5.86)3,  and  (2.5.82)  at  v=l  and  using  (2.5.83), 
(2.5.89)  becomes 

~f-(l)  =  24  [Wj  (3, 3, 1, 0, 1)  +  W2(3, 3, 1, 0, 1)  ]  -  4-^- 1  , 

dv  dv  v=i 

where  the  last  term  in  (2.5.91)  is 

dW5  dW5  dlk 

-4—^1  =-4-r— 1  =  16W55  • 

dv  v= 1  dlk  dv  v=  1 

On  substituting  (2.5.92)  into  (2.5.91),  we  obtain 
d2W 


dv2 


(1)  =  4[  6(W!  +  W2)  +  4W55]  li!  =i2  =3 

13  =  1,14  =  0,15  =  1 


(2.5.87) 

(2.5.88) 

(2.5.89) 

(2.5.90) 
equation 

(2.5.91) 

(2.5.92) 

(2.5.93) 


Comparing  (2.5.93)  with  (2.5.26)  and  noting  (2.5.27),  we  obtain  the  result  claimed  in 
(2.5.84),  i.e., 
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5~(1)  =  4E  ( >  0) ,  (2.5.94) 

dv 

where  E  is  Young’s  modulus  in  the  direction  of  anisotropy  for  infinitesimal  deformations 
of  an  incompressible  transversely  isotropic  material. 

Consider  now  the  stored-energy  w  for  such  a  linearly  elastic  material  given  by 
(2.4.8).  The  deformation  gradient  tensor  F  corresponding  to  the  invariants  (2.5.82)  is 
written  as 


F  = 


v“2  0  0 
0  v  0 
0  0  v 


(2.5.95) 


Thus  the  infinitesimal  strain  tensor  E  whose  components  Ey  appear  in  (2.4.8)  is 
determined  for  the  types  of  deformation  of  concern  from  (2.5.95),  (2.3.2)  as 


E  = 


v-2  - 1  0  0 

0  v-1  0 

0  0  v-1 


(2.5.96) 


so  that  (2.4.8)  takes  the  form 

w  =  w(v)  =  \lt  [(V2  -  l)2  +  2(v  -  l)2]  +  2(pL  -  \iT)  (v_ 2  -  l)2  +  Vi  y  (v~2  -  l)2  .  (2.5.97) 
A  simple  calculation  gives 

£(v)  =  \iT  [-4v-3(v-2  -  1)  +  4(v  -  1)]  -  4[2(hl  -  pT)  +  Vi  y]  w~3(v~2  -  1) 

=  4px  [~v-5  +  v-3  +  v  -  1]  -  4[  2(hl  -  (ij)  +  Vi  y]  (v“5  -  v~3 ) ,  (2.5.98) 

so  that 

^(v)  =  4pT  [5v-6  -  3V"4  +  1  ]  -  4[  2(|iL  -  pT)  +  Viy]  (-5v"6  +  3v"4) . 


(2.5.99) 
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On  letting  v  =  1  in  (2.5.99),  we  obtain 


d2w 

dv2 


(1)  =  4(3(J.t)  +  8[  2(|iL  -  Ji.T)  +  */2  y] 


or 


d2w 

dv2 


(1)  =  4  (4|xl  +  Y  -  M-t)  • 


(2.5.100) 


On  comparing  (2.5.100)  with  (2.4.10)2,  we  obtain 


d2w 

dv2 


(1)  =  4E , 


(2.5.101) 


so  that  consistency  between  the  nonlinear  result  (2.5.94)  and  the  linear  (2.5.101)  gives 


‘TT'(i) =  T^(i)  ’ 

dv2  dvz 


(2.5.102) 


or  in  terms  of  the  generalized  Young’s  modulus,  (2.5.102)  is  the  consistency  requirement 


E(l)  =  E 


(2.5.103) 


(see  e.g.,  equation  (2.5.19)).  We  remark  that  (2.5.102)  would  also  follow  from 
considering  v=l  +  e,  0<e<l,  and  equating  the  first  non-zero  coefficients  in  the 
expansions  of  W(v)  and  w(v)  in  powers  of  e. 

2.5.3  Implications  for  a  special  class  of  nonlinearly  elastic  materials 

Consider  an  incompressible  transversely  isotropic  nonlinearly  elastic  material  with 
stored-energy 


W=-|[ai-3)  +  g(I2-3)  +  f(I5)], 


(2.5.104) 


where  p.  is  a  constant,  and  g  and  f  are  as  yet  arbitrary  functions  of  their  arguments.  The 
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material  model  (2.5.104)  is  a  generalization  to  transversely  isotropic  materials  of  the 
Rivlin-Saunders  material  for  incompressible  isotropic  nonlinear  elasticity 

W  =  Cl  (Ii  -  3)  +  g(l2  -  3) ,  Ci  constant .  (2.5.105) 

We  note  that  the  Rivlin-Saunders  material  reduces  to  the  classic  Mooney-Rivlin  and 
neo-Hookean  materials  when  g  =  c2(l2  -3),  c2  constant,  and  g  =  0,  respectively.  The 
proposed  strain  energy  (2.5.104)  has  no  dependence  on  I3  by  virtue  of  incompressibility, 
and  no  dependence  on  I4  since,  as  stated  earlier,  our  concern  is  with  those  deformations 
for  which  the  invariants  Ii,  •  •  •  ,  I5  can  be  given  by  (2.5.82). 

From  (2.5.104),  we  see  that  we  must  have 

g(0)  +  f(l)  =  0  (2.5.106) 

for  the  stored-energy  to  be  zero  in  the  undeformed  state 

Ii  =  I2  =  3  ,  I5  =  1 .  (2.5.107) 

In  addition,  the  further  normalization  condition  (2.5.21)  implies 

f'(l)  =  0,  (2.5.108) 

where  the  prime  notation  denotes  differentiation  with  respect  to  argument.  For  the 
material  (2.5.104),  the  conditions  (2.5.26),  (2.5.27)  imply 

E  =  £L[3  +  3g'(0)  +  2f "(1) ]  >  0 ,  (2.5.109) 

or  on  employing  (2.5.82)  to  write  (2.5.104)  in  terms  of  v,  the  result  (2.5.94)  gives 

d2W 

-5-^(1)  =  4fi [3  +  3g'(0)  +  2f"(l)]  > 0 .  (2.5.110) 

dvz 


In  addition,  the  conditions  (2.5.54),  (2.5.55),  and  (2.5.63),  (2.5.64)  imply  that  for  the 
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material  (2.5.104), 

U  =  fL  =  pi[H-g'(0)]  >0.  (2.5.111) 

Therefore,  the  conditions  on  the  material  (2.5.104)  for  normalization  and  physical 
reasonability  are  (2.5.106),  (2.5.108),  (2.5.109),  and  (2.5.111).  We  note  from  (2.5.111) 
that  if  g  =  0,  or  if  simply  g'(0)  =  0,  then 

fj.  =  p;  =  fL  (2.5.112) 

for  the  material  (2.5.104). 

A  special  case  of  the  material  model  (2.5.104)  with  g  =  0  is  employed  in  Chapters  3 
and  4  of  this  dissertation.  This  is  given  by 

W  =  -^  [  (Ii  -  3)  +  a(l|  -  2I5  +  1)]  ,  (2.5.1 13) 

where  we  have  |X  >  0  by  virtue  of  (2.5.1 12)  and  the  comment  preceding  it,  and  the  choice 

f(I5 )  =  ads  -  2I5  +  1)  (a:  constant)  (2.5.1 14) 

can  easily  be  seen  to  satisfy  the  normalization  conditions  (2.5.106)  and  (2.5.108).  Thus 
the  only  remaining  condition  that  (2.5.113)  must  satisfy  is  (2.5.109),  which  implies 

E  =  |X  (3  +  4a)  >  0 .  (2.5.115) 

By  (2.5.109),  (2.5.110),  and  (2.5.115),  we  also  have 

-^-^(1)  =  4\i  (3  +  4a)  (2.5.116) 

dvz 

for  the  material  model  (2.5.113).  Finally,  from  (2.5.112),  (2.5.115),  we  determine  that 


we  must  have 
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a>-3/4.  (2.5.117) 

We  remark  that  when  a  =  0  in  (2.5.113),  we  recover  the  neo-Hookean  (isotropic) 


material. 
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Chapter  3 

HOMOGENEOUS  ANISOTROPIC  SPHERE 


3.1.  Introduction 

As  stated  in  Chapter  1,  our  aim  in  this  dissertation  is  to  examine  the  effects  of 
anisotropy  on  the  problem  of  cavitation  in  nonlinearly  elastic  incompressible  materials. 
This  problem  is  examined  for  a  homogeneous  sphere  in  this  Chapter.  In  Section  3.2,  we 
formulate  the  basic  boundary  value  problem  of  concern.  We  consider  a  sphere  composed 
of  an  incompressible  anisotropic  nonlinearly  elastic  material,  which  is  transversely 
isotropic  about  the  radial  direction.  The  sphere  is  subjected  to  a  prescribed  uniform 
radial  tensile  dead-load  po  on  its  boundary.  In  Section  3.3,  it  is  shown  that  one  solution 
to  this  problem,  for  all  values  of  po,  corresponds  to  a  trivial  homogeneous  state  in  which 
the  sphere  remains  undeformed  but  stressed.  However,  for  sufficiently  large  values  of 
Po,  one  has  in  addition  other  possible  radially  symmetric  configurations  involving  an 
internal  traction-free  cavity.  Such  solutions  bifurcate  from  the  homogeneous  solution  at 
Per,  a  critical  value  of  p0.  The  existence  of  these  bifurcated  solutions  depends  on  the 
constitutive  law  of  the  material  (Section  3.4).  In  Section  3.5,  a  specific  anisotropic 
material  model  (which  may  be  viewed  as  a  generalization  of  the  classic  isotropic  neo- 
Hookean  model)  is  considered  and  explicit  results  are  obtained  for  the  relationship 
Po  =  Po(c)  between  the  applied  tensile  load  p0  and  the  deformed  cavity  radius  c.  In 
contrast  to  the  situation  for  a  neo-Hookean  sphere,  where  bifurcation  is  always  to  the 
right,  bifurcation  here  may  occur  locally  either  to  the  right  (supercritical)  or  to  the  left 
(subcritical)  depending  on  the  degree  of  anisotropy.  In  Section  3.6,  the  stability  of  the 
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foregoing  solutions  is  examined  using  an  energy  minimization  approach.  The  cavitation 
solutions  are  shown  to  be  the  only  stable  solutions  for  sufficiently  large  loads.  When 
bifurcation  at  po  =  pCT  occurs  locally  to  the  right,  a  smooth  cavitation  occurs  with  cavity 
radius  increasing  continuously  (from  zero)  as  po  increases.  When  bifurcation  is  to  the 
left,  a  cavity  of  finite  radius  appears  at  a  transition  load  which  may  be  less  than  Pc-.  Such 
a  discontinuous  change  in  stable  equilibrium  configurations  is  called,  as  in  [20],  "snap 
cavitation",  by  analogy  with  snap-through  buckling  in  structural  mechanics.  (In  Sections 
3.5  and  3.6,  for  simplicity  we  focus  mainly  on  results  pertaining  to  the  specific  material 
model  described  in  Section  3.5.  In  the  next  Chapter,  a  more  general  treatment  will  be 
given  for  the  problem  of  cavitation  in  a  composite  anisotropic  sphere.  As  noted  in 
Chapter  1,  general  results  for  the  homogeneous  sphere  may  thus  be  recovered  from  those 
in  Chapter  4  by  the  proper  choice  of  the  parameters  of  Chapter  4.)  In  Section  3.7,  the 
stress  distribution  is  described.  When  the  cavity  is  nucleated  smoothly  (as  described 
above),  the  dominant  stress  variation  is  confined  to  a  narrow  boundary  layer  near  the 
cavity  wall.  Such  boundary  layers  have  previously  been  encountered  in  cavitation 
problems  for  isotropic  materials  ([6, 21]). 

3.2.  Problem  formulation 

We  are  concerned  with  a  sphere  composed  of  a  homogeneous  incompressible 
anisotropic  nonlinearly  elastic  material.  We  denote  the  interior  of  the  sphere  in  its 
undeformed  configuration  by  D0  =  {(R,  0,  <b)l  0  <  R  <  B,  0  <  0  <  27t,  0  <  <  7t}.  The 

sphere  is  subjected  to  a  prescribed  uniform  radial  tensile  dead-load  of  magnitude  p0  on 
its  boundary  R  =  B.  We  assume  that  the  resulting  deformation,  which  takes  the  point  in 
D0  with  spherical  polar  coordinates  (R,  0,  O)  to  the  point  (r,  0,  <|>)  in  the  deformed 
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configuration,  is  radially  symmetric.  Thus  the  deformation  has  the  form  0  =  ©,<(>  =  4>, 

r  =  r(R)  >  0,  0  <  R  <  B;  r(0f)>0,  (3.2.1) 


where  r(R)  is  to  be  determined. 

The  deformation  gradient  tensor  F  associated  with  (3.2.1),  referred  to  spherical 
polar  coordinates,  is  given  by 


F  = 


(R),  r(R)/R,  r(R)/R 


(3.2.2) 


where  ?  (R)  = 


Incompressibility  then  requires  that  J  =  detF=l,  which  upon 


integration  yields 


r(R)  =  (R3+c3)*, 


(3.2.3) 


where  c  >  0  is  a  constant  to  be  determined.  If  it  is  found  that  c  =  0,  (3.2.3)  implies  that 
the  body  remains  a  solid  sphere  in  the  current  configuration.  However,  if  c  is  found  to  be 
positive,  then  r(0+)  =  c  >  0  and  so  there  is  a  cavity  of  radius  c  centered  at  the  origin  in  the 
current  configuration.  In  this  event,  the  cavity  surface  is  assumed  to  be  traction-free. 

The  strain-energy  density  per  unit  undeformed  volume  for  an  elastic  material  that  is 
transversely  isotropic  about  the  Xi  direction  is  given  by  (see  e.g.  [27]) 

W  =  W(I1,I2,I3,I4,I5)  (3.2.4) 

where 


Ii=trC, 

12  =  1/2[(trC)2-trC2], 

13  =  det  C  , 

14  =  C22  +  C23  , 

15  =  Cn  , 


(3.2.5) 
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C  is  the  right  Cauchy-Green  deformation  tensor  C  =  FT  F,  and  Xj  (i  =  1,  2,  3)  are 
coordinates  of  a  point  in  the  deformed  configuration  associated  with  an  orthonormal 
basis.  For  incompressible  materials,  we  must  have  I3  =  J2  =  1.  The  corresponding 
response  equation  for  the  Cauchy  stress  tensor  T  for  transversely  isotropic 
incompressible  materials  is  [27] 

T  =  -pi  +  2(Wx  B  -  W2  B_1  +  W4  M  +  W5  N),  (3.2.6) 

where  B  is  the  left  Cauchy-Green  deformation  tensor  B  =  F  FT,  1  is  the  unit  tensor,  p  is 
the  unknown  hydrostatic  pressure  associated  with  the  incompressibility  constraint  J  =  1, 
and  Wq  =3W/8lq  (q  =  1,  2,  4,  5).  The  deformation  tensors  M,  N  associated  with  the 
anisotropy  are  given  by 

My  =  Cla  (Fia  FjX  +  Fja  Fji )  (a  =  2,  3)  (3.2.7) 

and 

Ny=FilFjl,  (3.2.8) 

where  i,  j  =  1,  2,  3,  and  summation  over  repeated  indices  is  used.  We  shall  assume  that 
the  strain-energy  W  vanishes  in  the  undeformed  state  where  Ii  =  3, 12  =  3, 13  =  1, 14  =  0, 
and  I5  =  1,  so  that  we  have  the  normalization  condition 

W(3,  3,  1,  0,  1)  =  0.  (3.2.9) 

We  also  assume  that  in  the  undeformed  state  the  initial  stress  is  a  hydrostatic  pressure,  so 
that  (3.2.6)  gives  the  further  normalization  condition  (cf.  [11],  p.  292) 

W5(3,  3,  1,  0,  1)  =  0 .  (3.2.10) 

In  (3.2.5)  -  (3.2.8)  we  now  let  (r,  0,  §)  be  associated  with  the  indices  (1,  2,  3)  so  that 
the  sphere  is  transversely  isotropic  about  the  radial  direction.  Corresponding  to  the 
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deformation  field  (3.2.1),  we  thus  obtain 


C  =  B  =  diag  [r2,  r2/R2,  r2/R2j  , 

(3.2.11) 

M  =  0 , 

(3.2.12) 

N  =  diag  ^r2,  0,  oj  , 

(3.2.13) 

lx  =  r  2  +  2  r2/R2  , 

(3.2.14) 

I2  =  r4/  R4  +  2r  2  r2/  R2  , 

(3.2.15) 

I3  =  r  2  r4/  R4  =  1 , 

(3.2.16) 

* 

O 

II 

(3.2.17) 

I5=?2. 

(3.2.18) 

Substitution  from  (3.2.11)  -  (3.2.13)  into  (3.2.6)  yields  the  nonzero  components  of  the 
Cauchy  stress  T  as 

T„  =  -p  +  2  (r  2  Wx  -  r  "2  W2  +  r  2  W5)  (3.2.19) 

and 

Tee  =  Th  =  -p  +  2(r  -1  Wx  -  r  W2) ,  (3.2.20) 

where  the  Wt  (t  =  1,  2,  5)  are  evaluated  at  the  values  of  the  invariants  given  by  (3.2.14)  - 
(3.2.18). 

The  equilibrium  equations,  in  the  absence  of  body  forces,  are 

div  T  =  0  ,  (3.2.21) 

which  in  the  present  case  reduce  to 

aTn-  2 

_iL  +  A(Trr_Tee)  =  0, 


(3.2.22) 
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_J_  dTee  _ 
r  ae 

1 

rsinG  3<j) 


(3.2.23) 

(3.2.24) 


Equation  (3.2.4),  in  view  of  equations  (3.2.14)-(3.2.18),  implies  that  W,  and  thus  Wq,  are 
independent  of  the  angular  variables.  Thus  (3.2.20),  together  with  (3.2.23)  and  (3.2.24), 
show  that  p  =  p(r)  only.  Since  r  =  r(R),  it  is  convenient  for  us  to  consider  T  =  T(R)  rather 
than  the  more  conventional  T  =  T(r).  The  remaining  equilibrium  equation  (3.2.22),  on 
using  the  chain-rule,  becomes 


dTg 

dR 


Tee)  =  0, 


0<R<B . 


(3.2.25) 


In  view  of  (3.2.3),  (3.2.19),  (3.2.20),  equation  (3.2.25)  is  a  first-order  nonlinear  ordinary 
differential  equation  for  the  pressure  p(R). 


The  dead-load  boundary  condition  now  requires  that 


Tjt(B)  =  p0 


(3.2.26) 


where  the  constant  po  >  0  is  prescribed.  We  note  that  the  boundary  conditions  of 
vanishing  shear  tractions  are  satisfied  identically. 

Thus,  the  boundary  value  problem  to  be  solved  is  the  following:  For  a  prescribed 
value  of  the  dead-load  traction  po  >  0,  we  seek  a  pressure  field  p(R)  and  a  constant  c>0 
such  that  (3.2.25),  (3.2.26)  are  satisfied  where  Tg,  Too  are  given  by  (3.2.19),  (3.2.20), 
(3.2.3).  In  addition,  ifc  >  0,  then  the  condition  for  a  traction-free  cavity  surface 

Tg(0)  =  0  (3.2.27) 


must  also  hold. 
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3.3.  Solutions 

On  using  the  normalization  condition  (3.2.10),  it  may  readily  be  shown  that  one 
solution  of  the  foregoing  problem  for  all  values  of  po  is 

p(R)  =  2[Wi  (3,  3,  1,  0,  1)  -  W2(3,  3,  1,  0,  1)]  -  p0  ,  c  =  0 .  (3.3.1) 

This  corresponds  to  the  trivial  homogeneous  state  of  deformation 

r(R)  =  R ,  (3.3.2) 


with  corresponding  stresses  Tn.  =  T0e  =  =  p0,  so  that  the  sphere  remains  undeformed 

under  a  state  of  hydrostatic  pressure. 

We  next  describe  solutions  for  which  c  >  0,  corresponding  to  the  presence  of  a 
traction-free  cavity  at  the  origin.  For  this  purpose,  it  is  convenient  to  introduce  the 
notation  (cf.  Ball  [5]), 


v  =  v(R)  =  r(R)/R  =  (1  + 

R 

-2  • 
v  =r 


(3.3.3) 


where  the  second  equation  in  (3.3.3)  results  from  the  definition  of  v  and  (3.2.16).  Using 
(3.3.3),  the  invariants  (3.2.14)  -  (3.2.18)  can  be  written  as 


11  =  v-4  +  2v2  , 

12  =  v4  +  2v-2  , 

13  =  l, 


I4  =  0, 
I5=v-4. 


(3.3.4) 


On  using  (3.2.19),  (3.2.20),  we  adopt  an  approach  developed  in  [20]  and  rewrite  the 
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differential  equation  (3.2.25)  in  the  form 

—  J  2  [  v-4  Wi  (v-4  +  2v2,  v4  +  2v-2,  1,  0,  v-4)  -v4  W2(v~4  +  2v2,  v4  +  2v-2,  1,  0,  v-4) 
dR 

k. 

*« 

+  v-4  W5  (v-4  +  2v2,  v4  +  2v-2,  1,  0,  v-4)  ]  -  p(R)  - 

+  [(v-4  -  v2)  Wi  (v-4  +  2v2,  v4  +  2v~2,  1,  0,  v-4)  (3.3.5) 

R 

+  (v-2  -  v4)  W2  (v-4  +  2v2,  v4  +  2v-2,  1,  0,  v-4) 

+  v-4  W5  (v-4  +  2v2,  v4  +  2v~2,  1,  0,  v-4)]  =  0on0<R<B. 

On  integration  of  (3.3.5),  we  have 


P(R)-P(0)  =  2|v-4(R)(Wi+W5)Ir-v4(R)W2Ir|  +  4J(R),  0<R<B,  (3.3.6) 

where 


J(R)  =  J  j^|v_7(s) -  v_1(s)  W!ls  +  [v-5(s)-v(s)  W2Is+v-7(s)W5|J  y  .  (3.3.7) 


In  (3.3.6),  (3.3.7)  the  W/  (/  =  1,  2,  5)  are  evaluated  at  the  values  of  the  invariants  given 
by  (3.3.4),  and  for  simplicity  we  have  introduced  the  notation 

W/  I  =  W /  (v— 4 (s)  +  2v2(s),  v4(s)  +  2v-2  (s),  1,  0,  v_4(s)) .  (3.3.8) 

S 

Using  (3.3.3),  (3.3.6),  the  radial  stress  (3.2.19)  becomes 


Trr(R)  =  -p(0)-4J(R),  0<R<B  . 


(3.3.9) 


The  traction-free  cavity  surface  condition  (3.2.27),  together  with  (3.3.9)  and  J(0)  =  0, 
now  yield 


so  that 


P(0)  =  0 


(3.3.10) 
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T„ (R)  =  -4J(R) ,  0  <  R  <  B  .  (3.3. 1 1) 

Finally,  by  virtue  of  (3.3.11),  (3.3.3),  we  see  that  the  dead-load  boundary  condition 
(3.2.26)  at  the  outer  surface  R  =  B,  is  satisfied  if 

-  4J(B)  =  p0  [v(B)]“2  .  (3.3.12) 

The  condition  (3.3.12)  may  be  written  in  a  compact  fashion  on  utilizing  the  change 

of  variables  s  ->  v  in  the  integral  (3.3.7).  From  (3.3.3)  it  is  seen  that  this  change  of 
variables  is  one-to-one  and  invertible  if  and  only  if  c  >  0.  Introducing  the  function 

W  (x)  =  W(x-4  +  2x2,  x4  +  2x-2,  1,  0,  x"4) ,  (3.3.13) 

and  adopting  the  notation 

Wi  (x)  =  -f-W(x),  (3.3.14) 

dx 

it  is  shown  in  Appendix  A  that  one  may  write  (3.3.12)  as 

7  Wx  (v) 

“5 - dv,  c > 0 . 

J3  (3.3.15) 

B3  ’ 

We  remark  that  (3.3.15),  here  established  for  the  anisotropic  sphere,  has  formally  the 
same  structure  as  the  analogous  relation  for  the  isotropic  sphere  (see  equation  (5.18)  of 
[5]  or  equation  (22)  of  [6]).  For  a  given  dead-load  po,  solutions  involving  a  traction-free 
internal  cavity  of  radius  c  exist  provided  that  c  is  a  positive  root  of  (3.3.15).  Substitution 
from  (3.3.10)  into  (3.3.6)  gives  the  associated  pressure  field  as 

P(R)  =  2[v-4(R)(W1+W5)Ir-v4(R)W2Ir]+4J(R),  0<R<B.  (3.3.16) 
In  summary,  we  have  seen  that  for  all  values  of  the  applied  dead-load  traction  p0, 
one  obtains  the  trivial  solution  (3.3.1)  corresponding  to  the  homogeneous  state  of 
deformation  (3.3.2).  Moreover,  if  positive  roots  c  of  (3.3.15)  exist,  then  one  obtains  the 
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additional  solutions  involving  a  traction-free  internal  cavity  described  above. 


3.4.  The  critical  load  and  bifurcation 

The  critical  load  pcr  at  which  an  internal  cavity  may  be  initiated  is  found  by 
formally  letting  c  ->  0+  in  (3.3.15),  and  so 


7  Wi(v) 

Pc  =  |  dv .  (3.4.1) 

Since  the  integral  (3.4.1)  is  improper,  pcr  may  or  may  not  be  finite,  and  thus  cavitation 
may  or  may  not  take  place. 

As  regards  the  lower  limit  in  (3.4.1),  it  is  shown  in  Chapter  2  (see  equations 


(2.5.87),  (2.5.94))  that 


W!(l)  =  --^(1)  =0,  -W(1)  =4E, 
dv  dv 


(3.4.2) 


where  E  is  Young’s  modulus  in  the  radial  direction  associated  with  infinitesimal 
deformations  of  the  transversely  isotropic  material  at  hand.  Thus  by  1’  Hopital’s  rule,  the 
limit  of  the  integrand  in  (3.4.1)  is  finite  as  v  — >  1.  Thus,  as  in  the  isotropic  case,  the 
finiteness  of  pCT  depends  on  the  behavior  of  W(v)  for  large  values  of  the  stretch  v. 
Suppose  that  the  strain-energy  density  per  unit  undeformed  volume  for  the  transversely 
isotropic  elastic  material  can  be  written  in  the  form 


W(v)  =  £L_m  v  m  +  —  +  a_jv  1 +ao +  ai  v  + —  +  anvn  ,  m,  n>0,  m  +  n>0,  (3.4.3) 


so  that 


Wj  (v)  =  -m  a_mv  m  1 - a_jv  2  +  ai  +  2a2V H - f  nanvn  1  .  (3.4.4) 
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From  (3.4.4)  we  see  that  for  large  values  of  the  stretch  v,  the  integral  (3.4.1)  will 
converge  if 

vn-4  <  v— 1  .  (3.4.5) 


Thus,  for  a  strain-energy  density  which  can  be  written  in  the  form  (3.4.3),  we  must  have 


n<3 . 


(3.4.6) 


A  similar  condition  was  previously  determined  for  the  incompressible  isotropic  sphere 
(see  equations  (28),  (29)  of  [6]).  For  more  general  W,  the  finiteness  of  pcr  in  (3.4.1)  also 
requires  restrictions  on  the  rate  of  growth  of  W  for  large  stretch. 

We  now  consider  the  local  character  of  the  bifurcation  at  po  =  pcr  by  analyzing  the 
curve  Po  =  Po(c),  given  by  (3.3.15),  for  small  values  of  c.  A  Taylor  expansion  of 
(3.3.15)  about  c  =  0  shows  that 


where 


Po  =  pcr+k 


B 


+  o  (c3)  asc-»0 


(3.4.7) 


k  =  lim  dpo/dc3  . 
c  — >0 


(3.4.8) 


From  (3.3.15)  one  obtains 


i  2  1  r 

k  =  -7  Per  -  bm 

3  3  v->l 


dW/dv 
v3-!  ' 


(3.4.9) 


On  using  l’Hopital’s  rule,  we  write  k  as 


k 


2 

3 


1  d2W(l) 
Pcr  6  dv2 


(3.4.10) 


where  from  (3.4.2)2  the  quantity  d2W(l)/dv2  is  finite  and  can  be  written  in  terms  of  the 


51 


infinitesimal  Young’s  modulus  in  the  radial  direction.  We  see  from  (3.4.7)  that  if  k  >  0, 
bifurcation  is  locally  to  the  right  (supercritical)  while  if  k  <  0,  bifurcation  is  locally  to  the 
left  (subcritical).  (As  an  illustration  of  these  two  cases,  see  Figures  3.3,  3.4,  respectively, 
for  values  of  p  =  c/  B  close  to  zero.) 

3.5.  Illustrative  example 

3.5.1.  A  class  of  anisotropic  incompressible  elastic  materials 

We  now  consider  a  specific  material  model  and  provide  an  explicit  solution  for  the 
bifurcation  problem  discussed  in  Sections  3.2-3.4.  We  recall  from  (3.2.17)  that  the 
invariant  I4  vanishes  identically  for  the  radially  symmetric  deformations  of  concern  here. 
Thus  one  of  the  simplest  strain-energy  densities  modeling  the  anisotropy  of  interest  in 
this  dissertation  may  be  taken  as 

W  =  -^[(I1-3)  +  f(I5)],  (3.5.1) 

where  jx  >  0  is  a  constant  and  f(Is)  is  an  arbitrary  function  of  the  fifth  invariant  subject 
only  to  the  restrictions  arising  from  the  normalization  conditions  (3.2.9),  (3.2.10),  namely 

f(l)  =  0,  f/(l)  =  0,  (3.5.2) 

where  the  prime  notation  denotes  differentiation  with  respect  to  argument.  In  the  special 
case  where  f  =  0,  (3.5.1)  reduces  to  the  well-known  neo-Hookean  isotropic  material 
model. 

Rather  than  proceeding  with  a  general  function  f  in  (3.5.1),  we  now  restrict  f  to  be  a 
polynomial  in  I5.  To  satisfy  the  conditions  (3.5.2),  it  can  be  seen  that  the  simplest 
nontrivial  representation  for  f  is  a  quadratic  function  of  I5,  f(I5)  =  alf  +  W5  +  c,  where  a. 
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b,  c  are  constants.  On  employing  (3.5.2)  this  representation  for  f  becomes 

f(I5)  =  a(I?-2I5  +  l),  (3.5.3) 

where  the  constant  a  >  0  is  a  dimensionless  parameter  which  serves  as  a  measure  of  the 
degree  of  anisotropy  of  the  material.  Thus,  from  (3.5.1)  and  (3.5.3)  we  obtain  the  strain- 
energy  density  function 

W=f[(II-3)  +  a(I?-2I5  +  l)],  (3.5.4) 

which  we  use  as  an  illustrative  example  for  the  remainder  of  this  paper.  When  a  =  0  in 
(3.5.4),  one  recovers  the  neo-Hookean  (isotropic)  material.  We  note  that  the  above 
condition  a>0  is  somewhat  more  restrictive  than  a >-%  determined  in  (2.5.117),  and 
that  p  in  (3.5.1)  and  (3.5.4)  is  equivalent  to  the  shear  moduli  in  infinitesimal 
deformations  of  the  transversely  isotropic  material  at  hand  (see  Section  2.5.3  and 
equation  (2.5.112)). 

The  response  of  the  material  (3.5.4)  to  certain  basic  pure  homogeneous 
deformations  will  now  be  discussed. 

(i)  Uniaxial  Tension: 

For  uniaxial  tension  in  the  Xj  -  direction  (the  axis  of  anisotropy)  we  have 

T22  =  T33  =  0,  Tu  #0,  =  X>1,  X2=X3  =  X"*  ,  (3.5.5) 

where  the  Xj  (i=  1,  2,  3)  are  the  principal  stretches.  For  this  deformation,  the  stress- 
stretch  relationship  for  the  material  (3.5.4)  is  given  by 

Tn(X)  =  \i  [(X2  -  X"1)  +  2a  X2(X2  -  1)] .  (3.5.6) 


It  is  readily  verified  that  dT^/dX  >  0  for  X  >  1  and  so  Tn(X)  is  monotone  increasing  in  X, 
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1  <  X  <  for  all  values  of  the  anisotropy  parameter  a  >  0.  The  relationship  (3.5.6)  is 
plotted  in  Figure  3.1  for  various  values  of  a. 

(ii)  Equibiaxial  Stress  (Stretch): 

We  next  consider  the  deformation 

Tn  =  0,  T22(X)  =  T33(X) * 0,  X2=X3=X,  Xi  =r2.  (3.5.7) 

We  note  that  for  the  bifurcation  problem  of  concern  in  this  paper,  we  have 

Tu  =  0,  Tee  =  T,^  *  0,  =  v>  =  y-2  (3.5.8) 

on  the  cavity  boundary,  so  that  locally,  a  state  of  equibiaxial  stress  (stretch)  exists  at  the 
cavity  wall.  For  the  deformation  (3.5.7),  the  stress-stretch  relationship  for  the  material 
(3.5.4)  is  given  by 

T(?0  ■  T2200  =  T33  OO  =  li  [(A2  -  r4)  -  2a  r4^"4  -1)] .  (3.5.9) 

This  relationship  is  plotted  in  Figure  3.2  for  various  values  of  the  anisotropy  parameter  a. 
It  can  be  seen  that  T(A)  is  monotone  increasing  in  A  for  1  <  A  <  «*»  provided  the  constant  a 
is  such  that  a  <  5.0  (approximately).  It  is  assumed  in  what  follows  that  0  <  a  <  5.  It  will 
be  seen  that  interesting  effects  of  anisotropy  on  cavitation  will  occur  for  values  of  a  well 
within  this  range. 

We  note  here  that  for  both  uniaxial  tension  and  equibiaxial  stress,  the  stress-stretch 
curves  involving  the  nominal  stress  S  =  (detF)F_1T  are  analogous  to  Figures  3.1,  3.2, 
respectively.  Finally,  we  also  remark  that  the  response  of  the  material  (3.5.4)  in  simple 
shear  is  exactly  that  of  the  neo-Hookean  material. 
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3.5.2.  Cavitation  solutions 

We  now  return  to  the  cavitation  problem  for  the  material  model  (3.5.1),  (3.5.2).  By 


virtue  of  (3.3.4),  (3.3.13),  we  rewrite  (3.5.1)  as 

W  =  it  lYv-4  4-  ?v2  - 


W  (v)  =  [(v-4  +  2v2  -  3)  +  f(v~4)] , 

Zr 


so  that,  recalling  the  notation  (3.3.14), 


Wi(v)  =  [-4v-5  +  4v  -  4v-5  f'(v~4)] 

Zr 


Substituting  (3.5.11)  into  (3.3.15),  we  obtain  the  relationship 


v-v— 5  v-5  nv"4) 

v3  -  1  v3-! 


f  * 

„3 

% 

1  +  -^r 

f 

B3 

k.  - 

C3  1 

O  +  TT)*  L 

(3.5.10) 


(3.5.11) 


(3.5.12) 


so  that  positive  roots  c  of  (3.5.12)  correspond  to  solutions  involving  a  traction-free 
internal  spherical  cavity.  We  rewrite  (3.5.12)  as 


where 


Po  =  2jx  (1  +  ~j)%  (Jt  +  J2)  > 


(3.5.13) 


Ji=  f 


v-v 
„3  _ 


(1  +  c3/B3)1/4  V  1 


(3.5.14) 


(1  +  c3/B3)m 


•y— sf/(v~4) 
v3  - 1 


(3.5.15) 


Upon  integration  of  (3.5.14),  (see  Appendix  C  for  evaluation  of  this  and  all  subsequent 
integrals),  we  obtain 
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Jl 


=  d  + 


sLx-* 

B3  ’ 


+  7»  +  F^ 


(3.5.16) 


We  now  consider  the  special  case  (3.5.4)  of  the  material  model  (3.5.1),  so  that  from 
(3.5.3),  (3.3.4)  we  have 

f'(v-4)  =  2a(v-4-l).  (3.5.17) 

Substituting  (3.5.17)  into  (3.5.15)  we  obtain 

“  -5  _  -9 

J2  =  2a  f  V  3  V  -  dv .  (3.5.18) 

(1  +  c3/B3)m  *  * 


Evaluation  of  (3.5.18)  (see  Appendix  C)  gives 


J 2  — 


+  T0  +  Fr"  +  i(1  +  F)r* 


(3.5.19) 


1 

V3 


arctan 


2(l+c3/B3)1/3  +  l 


V3 


+  arctan 


V3  (l  +  c3/B3)* 
2  +  (1  +  c3/B3)^ 


Thus,  for  the  material  (3.5.4),  the  relationship  p0  =  po(c)  between  the  prescribed  dead¬ 
load  traction  po  and  the  cavity  radius  c  is  given  by  (3.5.13),  (3.5.16),  (3.5.19). 

We  now  recall  the  discussion  in  Section  3.4  regarding  the  critical  load  pCT  at  which 
an  internal  cavity  may  be  initiated.  First  we  use  (3.3.4)  to  write  (3.5.4)  as 

W  (v)  =  [(v-4  +  2v2  -  3)  +  a  (v-8  -  2v"4  +  1)] .  (3.5.20) 

Comparing  (3.5.20)  with  (3.4.3),  we  have  n  =  2  and  thus  (3.4.6)  holds,  so  that  the  integral 
(3.4.1)  is  finite  for  the  material  (3.5.4).  (In  fact,  since  I5  =  v-4,  the  integral  (3.4.1)  will  be 
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finite  for  the  material  (3.5.1)  whenever  f(I5)  is  a  polynomial  in  I5  of  any  (non-negative) 
degree  and  so  the  anisotropy  does  not  affect  the  finiteness  of  pCT.)  On  letting  c  — >  0+  in 
(3.5.13),  (3.5.16),  (3.5.19),  it  can  be  shown  that 


Per  “  11 


40ft-5lV3 
30V3  J 


a 


(2.5  + 0.72a). 


(3.5.21) 


When  a  =  0  in  (3.5.4),  (3.5.21),  we  recover  the  result  for  the  neo-Hookean  material  (see 
e-g-  [5],  [6]). 

Thus,  if  we  consider  a  quasi-static  loading  process  in  which  the  solid  sphere  is 
subjected  to  a  tensile  dead-load  p0  that  increases  slowly  from  zero,  cavity  formation  and 
growth  is  described  by  the  relationship  p0  =  po(c)  given  by  (3.5.13),  (3.5.16),  (3.5.19). 
For  various  values  of  the  parameter  a,  the  critical  load  pcr  given  in  (3.5.21)  is  the  value  at 
which  the  curve  p0  =  po(c)  bifurcates  from  the  straight  line  c  =  0  corresponding  to  the 
trivial  solution. 


In  what  follows,  it  is  convenient  to  introduce  the  dimensionless  quantities 


P 


Po 

11 


(3.5.22) 


so  that  (3.5.22),  (3.5.13),  (3.5.16),  (3.5.19)  yield 
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p  =  P(p)  =  2(1  +  p3)%  J  (1  +  p  3fA  +  J  (1  +  p3)"473  +  2a  j  -  (1  +  p3)-^ 


+ --  a + p3fA  -  j  a + pV* + j  a + p3r5/3 + j  a + p3r8/3  (3.5.23) 


i  +arctan  Jii±e^ 

VJ  VJ  2  +  (l  +  p3)/5 


On  letting  p  — »  0+  in  (3.5.23)  one  obtains 


5  [4071-51^3 


p  _  ^  1 

"  2  +  30^3 


(3.5.24) 


which  is,  of  course,  exactly  (3.5.21).  By  analyzing  the  curve  (3.5.23)  in  the  (P,  p)  -  plane 
for  small  values  of  p,  we  can  determine  the  local  character  of  the  bifurcation  at  p  =  0. 
Letting 


P*(p3)  =  P(p), 


(3.5.25) 


the  Taylor  expansion  (3.4.7)  for  the  material  (3.5.4)  can  now  be  written  as 


P~(pJ)  =  Pa-  +  Kp3  +  o(p3)  as  p  — »  0 , 


(3.5.26) 


where  Pc-  is  given  by  (3.5.24),  and  (3.4.10),  (3.5.20),  (3.5.24)  give 


^  1  4071-  131^3 

K  =  —  +  - = -  a, 

3  45^ 


(3.5.27) 


Thus  the  curve  P  =  P*  (p3)  branches  from  the  line  p  =  0  ,  P  >  0,  at  the  point 
p  =  0,  P  =  PCT,  and  (3.5.26),  (3.5.27)  show  that  this  branching  (bifurcation)  is  to  the  right 


(see  Fig.  3.3)  if  K  >  0  and  is  to  the  left  (see  Fig.  3.4)  if  K  <  0.  From  (3.5.27),  these 
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conditions  on  K  may  be  written  asa<a<)Ora>ao  respectively,  where 

15^3 

ao  = - -  =  0.257.  (3.5.28) 

131V3 -40* 

(Note  that  ao  is  well  within  the  range  of  values  for  a  assumed  earlier  in  this  section.) 
Thus,  for  sufficiently  small  values  of  the  anisotropy  parameter  a,  bifurcation  is  locally  to 
the  right  as  for  the  neo-Hookean  material.  As  the  degree  of  anisotropy  increases  beyond 
a  =  ao,  the  bifurcation  character  changes  to  being  locally  to  the  left.  In  Figs.  3.3  and  3.4 
we  have  plotted  the  curve  (3.5.23)  for  two  representative  choices  of  the  parameter  a.  In 
Fig.  3.4,  there  exists  a  turning  point  (Pt,pt)  of  the  curve  at  which 

^  =  0.  (3.5.29) 

dp 

The  point  (Pt,  pt)  is  thus  determined  by  solving  (3.5.29)  for  the  root  pt  and  then  using 
(3.5.23)  to  give  Pt  =  P(pt). 

To  conclude  this  section,  we  examine  the  curve  (3.5.23)  for  large  values  of  p.  A 
detailed  analysis  given  in  Appendix  B  shows  that 

P  ~  2p  as  p  -» <» ,  (3.5.30) 

so  that  the  curve  (3.5.23)  is  asymptotic  to  a  straight  line  with  slope  =  2,  for  all  values  of 
the  parameter  a.  Thus  the  initial  character  of  the  bifurcation  at  (P,  p)  =  (Pc,.,  0)  does  not 
affect  the  large  p  asymptotic  behavior  of  P(p). 

3.6.  Energy  and  stability  of  solutions 

To  examine  the  stability  of  the  cavitation  solutions,  we  now  carry  out  an  energy 
analysis.  Returning  to  the  considerations  of  Sections  3.2  -  3.4  for  general  W,  we  find  that 
the  total  energy  associated  with  any  equilibrium  configuration  of  the  body  is  given  by 
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B 


E(c)  =  4jc  j  W  (r 2  +  2r2/R2,  r4/R4  +  2r  r  2/R2,  1,  0,  r  2)  R2  dR 
0 

-4tiB2  Po  [r(B)  -  B] . 


(3.6.1) 


The  first  term  on  the  right  in  (3.6.1)  is  the  total  strain-energy  while  the  second  term  is  the 
work  done  by  the  dead-load  tractions.  For  the  trivial  solution  (3.3.1)  with  deformation 
(3.3.2),  we  see  that 

E(0)  =  0  (3.6.2) 

on  using  the  normalization  (3.2.9).  We  will  see  in  this  section  that  equilibrium  solutions 
with  an  internal  cavity  of  radius  c  >  0  may  lead  to  either  positive  or  negative  values  for 
E(c).  Following  the  approach  of  [20],  we  now  extend  our  interpretation  of  the  energy 
function  E(c),  given  by  (3.6.1),  to  hold  for  all  values  of  c  >  0  and  not  just  for  those  values 
associated  with  equilibrium  solutions.  We  shall  seek  values  of  c  >  0  that,  for  a  given 
po  >  0,  minimize  the  extended  function  E(c). 

Recalling  the  notation  (3.3.13),  (3.3.3),  we  may  write  equation  (3.6.1)  as 


E(c)  =  4 to3  J  v2  (v3  -  l)"2  W(v-4  +  2v2,  v4  +  2v-2,  1,  0,  v-4) 

n  _L^3/n3\Vi 


dv 


(l+crVB^) 


-  47tB2p0  [r(B)  -B],  c  >  0  . 


(3.6.3) 


An  analogous  expression  has  been  obtained  in  Ball  [5]  for  the  isotropic  homogeneous 
sphere,  and  also  by  Horgan  and  Pence  [20]  for  the  isotropic  composite  sphere.  (See 
equations  (5.5)  and  (5.19)  of  [5]  and  equation  (2.27)  of  [20].) 

In  what  follows,  we  confine  attention  to  the  illustrative  example  described  in 
Section  3.5.  For  the  transversely  isotropic  sphere  with  strain-energy  density  given  by 
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(3.5.4)  (or  equivalently,  by  (3.5.20)),  the  expression  (3.6.3)  for  E(c)  yields 


E(c)  =  4nc3  ■£ 


f  (v3  -  1)  2  [(v  2  +  2v4  -  3v2)  +  a(v  6  -  2v  2  +  v2)]  dv 
+  c^3)*  (3.6.4) 


(l  +  cJ/BJ) 

-  4tcB2po  [(B3  +  c3)‘/3  -  B] ,  c  >  0  . 


Upon  integration  of  (3.6.4)  (see  Appendix  C  )  and  introduction  of  the  dimensionless 
quantity 


we  obtain 


2  = 


E 

(4/3)  tcB3  p.  ’ 


(3.6.5) 


2(p)  =  a(p)  +  Pp(p),  p>0 


(3.6.6) 


where 


a  (p)  =  3(1  +  p3)*  - 1  (1  +  p3)"14  -  ~  +  a-j 


K  ,  5  +  8p3  n  3r5/3 

9  10  1  v 


+  2p3  (1  +  p3r2/3  +  3(1  +  p3)"*  -  4(1  +  p3)%  +  j 


(3.6.7) 


4  3 


arctan 

2  +  (l  +  p)  V3 


and 


P  (p)  =  -  3[(1  +  p3)'A  -  1] .  (3.6.8) 

Thus  the  energy  minimization  problem  for  the  material  (3.5.4)  can  be  stated  in  a  purely 
algebraic  setting:  for  fixed  P  >  0 ,  we  seek  those  values  of  p  >  0  that  minimize  2(p). 

We  begin  by  examining  the  end-point  p  =  0  where  2(0)  =  0.  Using  a  superposed 
dot  to  denote  differentiation  with  respect  to  p,  it  can  be  shown  that 


2  (Of)  =  2  (Of)  =  0,  2  (Of)  =  6(Pcr  -  P) , 


(3.6.9) 


where  we  recall  that  for  the  material  (3.5.4),  PCT  is  given  by  (3.5.24).  Thus  p  =  0  is  a 
local  minimum  if  P  <  PCT  and  a  local  maximum  if  P  >  PCT.  Observe  that 

p  =  0,  P>0,  (3.6.10) 

corresponds  to  the  trivial  solution  (3.3.1). 

We  now  consider  the  possible  existence  of  internal  extrema  and  thus  seek  values  of 
p  >  0  for  which 

i(p)  =  a  (p)  +  P  P(p)  =  0 .  (3.6.11) 

From  (3.6.8),  P(p)  <  0  for  p  >  0  so  that  (3.6.1 1)  may  be  written  as 

P  =  -4^I,  p  >0.  (3.6.12) 

Rp) 

We  observe  that  the  expression  for  P(p)  given  by  (3.6.12),  derived  from  the  preceding 
energy  considerations,  can  be  seen  to  be  identical  to  the  expression  (3.5.23),  which  was 
obtained  from  the  solution  (3.3.15)  determined  directly  in  Section  3.3.  Thus,  as 
expected,  the  extremals  provided  by  (3.6.12)  correspond  to  the  equilibrium  solutions  with 
an  internal  cavity. 

Analyzing  the  dimensionless  energy  X(p)  given  by  (3.6.6)  -  (3.6.8)  near  p  =  0  and 
using  (3.5.26),  we  obtain 

X(p)  =  -Vi  Kp6  +  0(p9)  as  p  ->  0 ,  (3.6.13) 

where  K  is  given  by  (3.5.27).  Thus,  recalling  the  results  at  the  end  of  Section  3.5,  where 
the  curve  (3.6.12)  (or  equivalently  (3.5.23))  was  shown  to  bifurcate  from  the  trivial 
solution  (3.6. 10)  at  the  point  p  =  0,  P  =  PCT  ,  (3.6.13)  shows  that  X  <  0  on  this  curve  near 
p  =  0  if  the  bifurcation  is  to  the  right  (K  >  0  or  a  <  ao)  and  X  >  0  if  the  bifurcation  is  to 
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the  left  (K  <  0  or  a  >  ao).  Furthermore,  examining  (3.6.6)  -  (3.6.8)  for  large  values  of  p, 
and  using  (3.5.30),  an  analysis  similar  to  that  of  Appendix  B  shows  that 

I —  3p2  as  p  ~ ,  (3.6.14) 

so  that  X  will  always  be  negative  on  the  curve  (3.5.23)  for  sufficiently  large  values  of  p, 
independent  of  the  values  of  a.  Thus  the  initial  character  of  the  bifurcation  at 
(P,  p)  =  (PCT,  0)  does  not  affect  the  large  p  asymptotic  behavior  of  E  on  the  curve 
(3.5.23). 

The  results  obtained  thus  far  in  this  section  for  the  transversely  isotropic  sphere  with 
strain-energy  density  given  by  (3.5.4)  are  analogous  to  those  established  by  Horgan  and 
Pence  [20]  for  an  isotropic  composite  neo-Hookean  sphere.  (See  Section  3  of  [20]).  We 
now  state  two  properties  established  in  [20]  which  will  be  utilized  in  our  discussion  of 
stability  of  the  equilibrium  solutions  (3.6.10),  (3.5.23).  These  properties  may  be  verified 
using  arguments  analogous  to  those  given  in  [20,  pp.  72-74]. 

Property  A:  The  curve  P  =  P(p)  is  a  locus  of  local  minima  of  E(p)  for  fixed  P 
whenever  dP/dp  >  0,  and  is  a  locus  of  local  maxima  of  Z(p)  for  fixed  P  whenever 
dP/dp  <0. 

Property  B:  On  the  curve  P  =  P(p),  the  values  of  X(p)  are  decreasing  with  p 
whenever  dP/dp  >  0  and  are  increasing  with  p  whenever  dP/dp  <  0. 

We  begin  our  discussion  of  stability  by  considering  the  case  where  the  function  P(p) 
in  (3.5.23)  is  monotonically  increasing  so  that  the  bifurcated  curve  (3.5.23)  is  as  shown 
in  Fig.  3.3.  As  seen  from  Fig.  3.3,  the  number  of  equilibrium  solutions  depends  on  the 
value  of  P.  For  P  <  PCT,  the  only  solution  is  the  trivial  solution  p  =  0.  For  P  >  PCT,  there 
are  exactly  two  solutions,  namely  the  trivial  solution  p  =  0  and  in  addition  a  single 
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bifurcated  solution  involving  an  internal  cavity.  As  was  shown  previously  in  this  section, 
the  trivial  solution  minimizes  £  among  all  possible  (radially  symmetric)  configurations 
when  P  <  Per.  For  P  >  PCT,  the  trivial  solution  provides  a  local  maximum  for  £  among 
this  same  class  of  configurations,  whereas  the  bifurcated  solution  yields  a  local  minimum 
from  Property  A,  and  hence  the  absolute  minimum.  Thus  we  say  that  the  trivial  solution 
is  stable  for  P  <  Pcr.  For  P  >  Per,  the  trivial  solution  is  unstable  while  the  bifurcated 
solution  is  stable.  Thus  for  all  values  of  P  there  exists  exactly  one  stable  solution.  We 
now  consider  this  result  in  the  context  of  the  quasi-static  process  in  which  the  sphere  is 
subjected  to  a  tensile  dead-load  surface  traction  p0  =  P|i  which  increases  slowly  from 
zero.  As  po  is  increased  from  zero  to  Pcr  |i,  the  stable  equilibrium  configuration  is  that 
produced  by  the  homogeneous  deformation  (3.3.2).  At  po  =  Per  It  an  internal  void  is 
nucleated,  and  as  p0  is  further  increased  from  PCT  |i  ,  the  cavity  radius  c  =  Bp  increases 
from  zero  in  a  continuous  fashion. 

Analogous  results  to  the  above,  as  well  as  to  the  remainder  of  this  section,  have 
been  obtained  and  discussed  for  the  isotropic  composite  neo-Hookean  sphere  in  [20]. 
(See  Sections  3  and  4  of  [20]).  Our  subsequent  discussion  follows  closely  that  of  [20]. 

We  now  consider  the  case  where  the  function  P(p)  in  (3.5.23)  is  monotone 
decreasing  for  0  <  p  <  pt,  and  monotone  increasing  for  p  >  pt  as  in  Fig.  3.4.  As  seen 
from  Fig.  3.4,  the  number  of  equilibrium  solutions  depends  on  the  values  of  P  in  the 
following  way:  for  0  <  P  <  Pt,  the  only  solution  is  the  trivial  solution 
p  =  0;  for  Pt  <  P  <  Pcr,  there  are  exactly  three  solutions,  namely,  the  trivial  solution  p  = 
0  and  in  addition  two  bifurcated  solutions  each  involving  an  internal  cavity;  finally,  for 
P>Pc,.,  there  are  exactly  two  solutions,  namely  the  trivial  solution  and  a  single 
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bifurcated  solution  involving  an  internal  cavity.  As  noted  earlier,  the  trivial  solution 
provides  a  local  minimum  for  £  and  hence  is  stable  for  0  <  P  <  PCT,  while  for  P  >  PCT  the 
trivial  solution  provides  a  local  maximum  for  £  and  thus  is  unstable.  By  virtue  of 
Property  A,  the  bifurcated  curve  (3.5.23)  provides  a  local  maximum  for  £  if  0  <  p  <  pt 
and  hence  the  associated  solution  is  unstable,  while  for  p  >pt,  this  curve  provides  a  local 
minimum  for  £  and  hence  this  solution  is  stable.  Thus  if  0  <  P  <  Pt  or  if  P  >  PCT,  there  is 
exactly  one  stable  equilibrium  solution,  while  if  Pt  <  P  <  PCT  there  are  two  stable 
equilibrium  solutions.  To  differentiate  between  these  two  stable  solutions,  we  say  that  a 
solution  is  absolutely  stable  (AS)  if  it  provides  the  absolute  minimum  for  £  whereas  it  is 
said  to  be  metastable  (MS)  if  it  provides  a  local  minimum  for  £  which  is  not  the  absolute 
minimum.  To  address  the  transition  between  metastability  and  absolute  stability,  we 
refer  to  Fig.  3.5,  where  pt  and  Pt  are  as  before,  and  pc  denotes  the  value  of  p  >  0  for 
which  P(p)  =  Per-  Then  there  exists  a  unique  value  of  P,  namely  P==PN>  obeying 
Pt  <  Pn  <  Pcr,  such  that  for  P  <  Pn  the  trivial  solution  is  absolutely  stable,  while  for 
P  >  Pn  the  absolutely  stable  solution  occurs  on  the  bifurcated  branch.  Since  £  =  0  for  the 
trivial  solution  p  =  0,  the  preceding  result  follows  from  the  existence  of  a  unique  value  of 
p,  namely  p  =  Pn  »  satisfying  pt  <  pN  <  pc,  such  that  £  (p)  >  0  for  pt  <  p  <  Pn  and 
£  (p)  <  0  for  Pn  <  p  <  pc  on  the  curve  P  =  P(p)  given  by  (3.5.23).  The  existence  of  Pn 
follows  immediately  from  Properties  A  and  B. 

We  again  consider  the  previously  envisioned  quasi-static  loading  process  for  the 
transversely  isotropic  sphere  whose  equilibrium  solutions  give  rise  to  a  bifurcation 
diagram  of  the  type  shown  in  Fig.  3.5.  If  one  supposes  that  the  body  always  assumes  an 
absolutely  stable  configuration,  the  trivial  homogeneous  deformation  (3.3.2)  persists  for 
0  <  po  <  Pn  It-  As  po  passes  through  PN  |i,  the  foregoing  notion  of  stability  predicts  the 
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sudden  appearance  of  a  cavity  of  finite  radius  c  =  Bp^-  As  po  is  further  increased,  the 
cavity  radius  increases  in  a  continuous  fashion.  On  the  other  hand,  one  may  suppose  that 
the  trivial  homogeneous  deformation  persists  until  it  becomes  unstable  at  po  =  PCT  (i,  in 
which  case  the  initial  cavity  radius  is  predicted  to  be  c  =  Bpc,  with  subsequent  growth 
occurring  in  a  continuous  fashion  as  p0  is  further  increased.  Regardless  of  the  stability 
criterion  adopted,  one  obtains  a  discontinuous  transition  from  the  stable  trivial  solution 
p  =  0,  P  <  Pc-  to  the  stable  bifurcated  solution  P  =  P(p),  p  >  pt.  A  similar  discontinuous 
transition  was  also  obtained  in  [20],  [21]  for  a  composite  isotropic  incompressible  sphere 
and  by  Antman  and  Negron-Marrero  [4]  for  compressible  anisotropic  spheres  and 
cylinders.  Such  behavior  was  called  "snap  cavitation"  in  [20],  [21]  by  analogy  with  the 
snap-through  buckling  phenomenon  of  structural  mechanics. 

3.7.  Stress  distribution 


We  proceed  now  to  discuss  the  stress  distribution  in  the  anisotropic  sphere.  We 
begin  by  recalling  some  results  obtained  in  Sections  3. 2-3. 4  for  general  W.  At  the 
beginning  of  Section  3.3,  we  saw  that  the  stresses  corresponding  to  the  trivial  solution 
(3.3.1)  were  given  as 


Tn- =  Tee  =  =  po  .  (3.7.1) 

Thus,  prior  to  cavitation,  we  have  simply  the  constant  stress  distribution  (3.7.1). 
Subsequent  to  cavitation,  (3.3.11),  (3.3.7),  and  the  considerations  of  Appendix  A  give  the 
radial  stress  as 


oo 

T„(R)=  J 

(1  +  C3/B3)W 


Wi  (v) 
v3  - 1 


dv. 


(3.7.2) 
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where  we  recall  the  notation  (3.3.13),  (3.3.14).  Observe  that  for  fixed  R,  0  <  R  <  B, 
letting  c  — »  0+  formally  in  (3.7.2)  results  in  the  right-hand  side  of  (3.4.1),  so  that 

lim  Tff  (R)  =  pCT  ,  0<R<B.  (3.7.3) 

c  -»0+ 


We  shall  use  (3.7.3)  later  on  in  this  section.  Now,  from  (3.2.19),  (3.2.20),  (3.3.3),  we  can 
write  the  hoop  stresses  as 

Tee  (R)  =  T#  (R)  =  2  [(v2  -  v'4)  Wt  I R  -  (v~2  -  v4)W2 1 R  -  V4W5  I R]  +  Tn (R) ,  (3.7.4) 

where  T^R)  is  given  by  (3.7.2)  and  we  recall  the  notation  (3.3.8). 

In  what  follows,  we  confine  attention  to  the  illustrative  example  of  Section  3.5.  On 
comparing  (3.7.2)  and  (3.3.15),  it  follows  immediately  from  (3.5.13),  (3.5.16),  (3.5.19) 
that  for  the  material  (3.5.4)  we  have 


Tn-  (R)  =  2|i 
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2(1  +c3/R3)^  + 1 

€ 


+  arctan 


^3  (l+c3/R3)>/3 


2  +  (1+c3/R3)1/3 


Now,  from  (3.5.4)  it  is  readily  verified  that 


Wx  =  y  ,  W2=0,  W5  =  p a (I5  -  1) , 


(3.7.6) 


so  that  (3.7.6),  (3.7.4),  (3.3.4),  (3.3.3)  combine  to  give  the  hoop  stresses  for  the  material 
(3.5.4)  as 
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Tee(R)  -  T^R)  -  (i 
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+  T„(R), 


(3.7.7) 


where  ^(R)  is  given  by  (3.7.5). 

The  stresses  (3.7.5),  (3.7.7)  are  plotted  in  Figure  3.6  for  a  representative  value  of  the 
anisotropy  parameter  a  <  ao-  We  observe  from  Fig.  3.6  that  Tn.  is  a  monotone  increasing 
function  of  R  on  0  <  R  <  B,  while  T00  =  are  monotone  decreasing  functions  of  R  on 
this  range.  Notice  also  from  Fig.  3.6  that  Tee(R)=T^(R)  >  T„(R)  for  c  >  0.  This  latter 
fact  may  also  be  verified  analytically,  for  all  values  of  a,  as  follows:  from  (3.7.7)  we 
obtain 


T0e(R)-Trr(R)  =  T(X), 


(3.7.8) 


where  T  (X)  is  given  by  (3.5.9),  with  X  =  v(R)  = 
3.2,  T(X)  >  0  for  X  >  1,  so  that  (3.7.8)  implies  that 


1  + 


R3 


vs 


.  Thus,  as  seen  from  Fig. 


Tee(R)  =  T*(R)  >  T„(R)  for  c  >  0 .  (3.7.9) 

Similar  results  were  obtained  for  the  composite  neo-Hookean  isotropic  sphere  in  [21,  Fig. 
4]. 

Another  interesting  feature  concerning  the  stresses  Tn.(R),  T00(R),  T^(R) 
immediately  after  cavitation  is  the  presence  of  a  boundary  layer  near  the  cavity  wall 
when  smooth  cavitation  (i.e.,  K  >  0  or  a  <  ao)  has  taken  place.  Observe  from  (3.7.3), 
(3.5.21)  that  for  the  material  (3.5.4), 
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lim  Tit(R)  =  |I 
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d 

0  <  R  <  B  . 


(3.7.10) 


Since  the  boundary  condition  (3.2.27)  of  a  traction-free  cavity  surface  guarantees 
T„(0)  =  0,  the  radial  stress  thus  experiences  a  rapid  growth  in  a  narrow  region  near  the 
cavity  wall  for  applied  dead-loads  po  slightly  larger  than  pCT.  From  (3.7.7)  it  is  clear  that 
a  similar  boundary  layer  exists  in  the  stress  components  Tee,  T^.  These  boundary 
layers  in  the  stresses  are  illustrated  in  Figures  3.7  and  3.8  for  a  representative  value  of  the 
anisotropy  parameter  a.  Similar  boundary  layer  behavior  in  the  stresses  was  observed  in 
[21],  [6].  (See  Sec.  4  and  Fig.  5  of  [21],  Figs.  3  and  4  of  [6].)  As  noted  in  [21],  such 
severe  boundary  layers  for  the  stresses  do  not  occur  in  the  case  of  snap  cavitation 
(K  <  0  or  a  >  ao)  since  in  this  case  the  deformed  cavity  radius  jumps  to  a  finite  value 
upon  cavity  formation. 


70 


T(X) 

Figure  3.2.  Nondimensional  stress  — —  vs.  stretch  X  for  equibiaxial  stress  (stretch)  for 

ft 

the  anisotropic  material  with  strain-energy  density  given  by  (3.5.4). 


Figure  3.3.  Variation  of  the  deformed  cavity  radius  p  =  c/B  with  applied  dead-load 
traction  P  =  p0/\i  for  the  anisotropic  sphere  with  strain-energy  density  given  by  (3.5.4) 
with  a  =0.125. 
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Figure  3.4.  Variation  of  the  deformed  cavity  radius  p  =  c/B  with  applied  dead-load 
traction  P  =  p0/\i  for  the  anisotropic  sphere  with  strain-energy  density  given  by  (3.5.4) 
with  a  =  2.0. 


’„.(/?),  Tqq(R )  with  undeformed  radius  R 
material  with  strain-energy  density  given  by 
=  2.644,  p  =  c/B=  0.3. 
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Figure  3.7.  Variation  of  the  radial  stress  T^iR)  with  undeformed  radius  R  subsequent  to 
cavitation  for  the  anisotropic  material  with  strain-energy  density  given  by  (3.5.4)  with 
a  =  0.2.  Here,  Pcr  =  2.644. 
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Figure  3.8.  Variation  of  the  hoop  stress  Tqq(R)  with  undeformed  radius  R  subsequent  to 
cavitation  for  the  anisotropic  material  with  strain-energy  density  given  by  (3.5.4)  with 
a  =  0.2.  Here,  Pcr  =  2.644. 
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Chapter  4 

COMPOSITE  ANISOTROPIC  SPHERE 


4.1.  Introduction 

In  the  previous  chapter,  the  effects  of  material  anisotropy  on  cavitation  for 
incompressible  nonlinearly  elastic  spheres  was  examined.  That  investigation  was 
motivated  by  the  work  of  Antman  and  Negron-Marrero  [4]  on  compressible  transversely 
isotropic  nonlinearly  elastic  solids  and  by  the  studies  of  Horgan  and  Pence  [20-22]  on 
composite  incompressible  nonlinearly  elastic  spheres  with  isotropic  phases.  The  purpose 
of  the  present  Chapter  is  to  investigate  the  combined  effects  of  anisotropy  and  material 
inhomogeneity  on  cavitation.  A  more  complicated  sequence  of  cavitational  instabilities 
might  be  expected  to  occur  due  to  the  interaction  of  these  effects  and  we  show  that  this  is 
indeed  the  case. 

In  Section  4.2,  we  formulate  the  basic  boundary  value  problem  of  concern.  We 
consider  a  composite  sphere  composed  of  two  incompressible  anisotropic  nonlinearly 
elastic  phases,  each  of  which  is  transversely  isotropic  about  the  radial  direction.  The 
sphere  is  subjected  to  a  prescribed  uniform  radial  tensile  dead-load  po  on  its  boundary. 
In  Section  4.3,  it  is  shown  that  one  solution  to  this  problem,  for  all  values  of  po, 
corresponds  to  a  trivial  homogeneous  state  in  which  the  sphere  remains  undeformed  but 
stressed.  However,  as  seen  previously  for  the  homogeneous  anisotropic  sphere,  for 
sufficiently  large  values  of  p0,  one  has  in  addition  other  possible  radially  symmetric 
configurations  involving  an  internal  traction-free  cavity.  Again,  such  solutions  bifurcate 
from  the  homogeneous  solution  at  pCT,  a  critical  value  of  po-  The  possibility  for  these 
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bifurcated  solutions  to  exist  depends  only  on  the  constitutive  law  for  the  material  at  the 
core  of  the  composite  sphere  (Section  4.4).  In  Section  4.4,  we  also  give  conditions  to 
determine  whether  bifurcation  is  supercritical  or  subcritical.  In  Section  4.5,  the  stability 
of  the  foregoing  solutions  is  examined  using  an  energy  minimization  approach.  For  a 
composite  sphere  described  by  a  general  transversely  isotropic  strain-energy  density 
function,  W,  the  cavitation  solutions  are  shown  to  be  the  only  stable  (radially  symmetric) 
solutions  for  sufficiently  large  loads.  In  Section  4.6,  we  present  an  example  which 
illustrates  explicitly  the  preceding  results.  The  specific  anisotropic  material  model 
described  in  Chapter  3,  Section  5,  is  considered  for  each  phase  of  the  composite,  and 
explicit  results  are  obtained  for  the  relationship  p0  =  po(c)  between  the  applied  tensile 
load  po  and  the  deformed  cavity  radius  c.  We  also  obtain  an  explicit  result  for  the 
relationship  between  the  total  energy  L  and  the  cavity  radius.  We  recall  that  for  the 
isotropic  composite  sphere  discussed  in  [20],  and  the  homogeneous  anisotropic  sphere 
described  in  Chapter  3,  Section  5,  when  bifurcation  at  p0  =  pcr  occurs  locally  to  the  right, 
a  smooth  cavitation  takes  place  with  cavity  radius  increasing  continuously  (from  zero)  as 
Po  increases.  Here  however  the  situation  is  quite  different.  We  exhibit  three  possible 
configurations  the  sphere  may  assume  when  bifurcation  occurs  locally  to  the  right.  First, 
a  smooth  cavitation  may  take  place.  Second,  a  snap  cavitation  may  occur  in  which  there 
is  a  cavity  of  finite  radius  upon  first  appearance  at  a  transition  load  which  is  less  than  pcr. 
Third,  the  sphere  may  undergo  a  smooth  cavitation  initially,  and  then  experience  a 
discontinuous  jump  in  cavity  radius.  That  is,  the  cavity  radius  "snaps"  from  one  finite 
Value  to  another  (larger)  finite  value  at  a  load  po  greater  than  pcr.  The  last  two 
phenomena  were  not  seen  for  the  spheres  discussed  in  [20]  and  Chapter  3.  When 
bifurcation  is  to  the  left,  the  sphere  experiences  a  snap  cavitation  so  that  again  a  cavity  of 
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finite  radius  appears  at  a  transition  load  which  is  less  than  pcr.  In  Section  4.7,  the  stress 
distribution  is  described.  When  a  cavity  is  nucleated,  in  the  case  of  smooth  cavitation,  a 
predominant  stress  variation  exists  in  a  narrow  boundary  layer  near  the  cavity  wall. 
Also,  the  effect  of  cavitation  at  the  center  on  possible  interface  debonding  is  considered. 
It  is  assumed  that  debonding  occurs  uniformly  whenever  the  normal  stress  at  the 
interface  reaches  a  threshold  value.  A  condition  on  the  strain-energy  is  given  so  that  in  a 
quasi-static  loading  process,  cavitation  relieves  the  interfacial  normal  stress  so  that 
subsequent  interfacial  debonding  is  precluded.  The  boundary  layers  in  the  stresses  as 
well  as  the  stress  relaxation  in  the  normal  stress  at  the  interface  subsequent  to  cavitation 
are  illustrated  for  the  specific  material  model  of  Section  4.6. 

4.2.  Problem  formulation 

We  are  concerned  with  a  sphere  composed  of  an  incompressible  anisotropic 
nonlinearly  elastic  material.  As  in  Chapter  3,  we  denote  the  interior  of  the  sphere  in  its 
undeformed  configuration  by  Dq  =  {(R,  0,  <t>)  I  0  <  R  <  B,  0  <  O  <  2tt,  0  <  d>  <  7t}.  The 
sphere  is  subjected  to  a  prescribed  uniform  radial  tensile  dead-load  of  magnitude  po  on 
its  boundary  R  =  B.  We  assume  that  the  resulting  deformation,  which  takes  the  point  in 
D0  with  spherical  polar  coordinates  (R,  0,  <&)  to  the  point  (r,  0,  <]))  in  the  deformed 
configuration,  is  radially  symmetric.  Thus  the  deformation  has  the  form  0  =  ©,  <J>  =  O, 

r  =  r(R)  >  0,  0  <  R  <  B;  r(0f)>0,  (4.2.1) 

where  r(R)  is  to  be  determined. 

As  in  Chapter  3,  the  deformation  gradient  tensor  F  associated  with  (4.2.1),  referred 
to  spherical  polar  coordinates,  is  given  by 


F  =  diag  [  r  (R),  r(R)/R,  r(R)/R] 
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(4.2.2) 

where  r  (R)  =  .  Incompressibility  then  requires  that  J  =  detF  =  l,  which  upon 

dR 

integration  yields 

r(R)  =  (R3  +  c3)'^,  (4.2.3) 

where  c  >  0  is  a  constant  to  be  determined.  As  in  Chapter  3,  if  it  is  found  that  c  =  0, 
(4.2.3)  implies  that  the  body  remains  a  solid  sphere  in  the  current  configuration. 
However,  if  c  is  found  to  be  positive,  then  r(0+)  =  c  >  0  and  so  there  is  a  cavity  of  radius 
c  centered  at  the  origin  in  the  current  configuration.  In  this  event,  the  cavity  surface  is 
assumed  to  be  traction-free. 

The  incompressible  material  of  concern  here  is  assumed  to  be  transversely  isotropic 
in  the  radial  direction,  and  possibly  inhomogeneous  in  the  radial  direction  as  well.  The 
strain-energy  density  per  unit  undeformed  volume  for  such  an  elastic  material  is  denoted 
by 

W  =  W(I1,I2,l3,l4,I5;R)  (4.2.4) 

where  the  explicit  dependence  of  W  on  the  radial  coordinate  R  reflects  the  possible 
material  inhomogeneity.  We  shall  assume  that  the  strain-energy  W  vanishes  in  the 
undeformed  state  where  Ii  =  3, 12  =  3, 13  =  1, 14  =  0,  and  I5  =  1,  so  that  we  have  the 
normalization  condition 

W(3,  3,  1,  0,  1;  R)  =  0  V  R  >  0 .  (4.2.5) 

As  in  Chapter  3,  we  also  assume  that  in  the  undeformed  state  the  initial  stress  is  a 
hydrostatic  pressure,  so  that  we  have  the  further  normalization  condition 


81 


W5(3,  3,  1,  0,  1;  R)  =  0  V  R  >  0 .  (4.2.6) 

We  shall  be  concerned  in  what  follows  with  the  case  of  a  composite  sphere  composed  of 
two  different  transversely  isotropic  incompressible  materials  perfectly  bonded  across  the 
interface  R  =  A  ( <  B  ).  The  material  properties  in  each  individual  phase  are  assumed  to 
have  smooth  radial  dependence.  Thus  W(Ii,  I2, I3, I4, 15 ;  R )  is  assumed  to  be  smooth 
on  0<R<A,A<R<B  while  suffering  a  jump  discontinuity  at  the  interface  R  =  A. 
Observe  that,  by  virtue  of  (4.2.3),  the  deformation  field  is  continuous  at  this  interface.  In 
the  remainder  of  this  section,  we  proceed  formally  and  assume  that  W  possesses 
sufficient  regularity  properties  to  permit  the  subsequent  analysis. 

The  response  equation  for  the  Cauchy  stress  tensor  T  for  transversely  isotropic 
incompressible  materials  is,  as  given  previously, 

T  =  -pi  +  2(W!  B  -  W2  B"1  +  W4  M  +  W5  N),  (4.2.7) 

where,  as  before,  B  is  the  left  Cauchy-Green  deformation  tensor  B  =  F  FT,  1  is  the  unit 
tensor,  p  is  the  unknown  hydrostatic  pressure  associated  with  the  incompressibility 
constraint  J  =  1,  and  Wq  =3W/3lq  (q  =  1,  2,  4,  5).  The  deformation  tensors 
B,  C  =  Ft  F,  M,  and  N,  and  the  invariants  Ij,  (i=l,...,5),  are  given  in  Chapter  3,  Section 
2,  in  general,  as  well  as  for  the  particular  deformation  and  transverse  isotropy  of  concern. 
Thus,  for  the  sphere  which  is  transversely  isotropic  about  the  radial  direction,  we  have, 
corresponding  to  the  deformation  field  (4.2.1), 


C  =  B  =  diag 


r2,  r 2/R2,  i^/R2 


(4.2.8) 


M  =  0 , 


(4.2.9) 
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N  =  diag  ^r2,  0,  oj  , 

(4.2.10) 

Ix  =?2  +  2r2/R2  , 

(4.2.11) 

I2  =  i4/  R4  +  2r  2  r2/  R2  , 

(4.2.12) 

I3  =  ?  2  r4/  R4  =  1 , 

(4.2.13) 

I4  =  o, 

(4.2.14) 

T  -?2 

I5  =r  . 

(4.2.15) 

Substitution  from  (4.2.8)  -  (4.2.10)  into  (4.2.7)  yields  the  nonzero  components  of  the 
Cauchy  stress  T  as 

T„  =  -p  +  2  (r  2  W!  -  r “2  W2  +  r  2  W5)  (4.2.16) 


Tee  =  T#  =  -p  +  2(r  Wx  -  r  W2) ,  (4.2.17) 

where  the  Wt  (t  =  1,  2,  5)  are  evaluated  at  the  values  of  the  invariants  given  by  (4.2.11)  - 
(4.2.15),  and  we  note  that  Wt  is  also  dependent  on  the  radial  coordinate  R.  The 
derivatives  in  (4.2.16),  (4.2.17)  are  understood  to  be  one-sided  derivatives  at  the  material 
interface  R  =  A. 


The  equilibrium  equations,  in  the  absence  of  body  forces,  are 

div  T  =  0 ,  (4.2.18) 


which  in  the  present  case  reduce  to 


dTn 

3r 


+ 


2 

r 


(Tn.  -  Tee)  =  0 , 


1  dTee 


(4.2.19) 


r  38 


(4.2.20) 
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1  3T* 


rsin0  3(f> 


=  0. 


(4.2.21) 


As  seen  in  Chapter  3,  Section  2,  we  again  have  Wq  independent  of  the  angular  variables, 
and  thus  it  follows  from  (4.2.17)  and  (4.2.20)  -  (4.2.21)  that  the  pressure  p=p(r)  only. 
There  may  be  a  jump  discontinuity  at  the  interface  R  =  A  in  p(r).  Elsewhere  the  pressure 
is  assumed  sufficiently  smooth.  Since  r  =  r(R)  ,  we  again  consider  T  =  T(R)  (and 
p  =  p(R)).  The  remaining  equilibrium  equation  (4.2.19),  on  using  the  chain-rule, 
becomes 


— ^-  +  —  (Tir-Tee)  =  0,  on  0<R<A,  A<R<B.  (4.2.22) 
dR  r 

Equation  (4.2.22)  is  a  first-order  nonlinear  ordinary  differential  equation  for  the  pressure 
p(R).  Also  for  equilibrium,  the  normal  traction  component  at  the  material  interface  must 
be  continuous,  i.e., 


Trr(A-)  =  T„(A+) . 

The  dead-load  boundary  condition  now  requires  that 

12 


Trr(B)  =  p0 


B 


r(B) 


(4.2.23) 


(4.2.24) 


where  the  constant  po  >  0  is  prescribed.  We  note  that  the  boundary  conditions  of 
vanishing  shear  tractions  are  satisfied  identically. 

Thus,  the  boundary  value  problem  to  be  solved  is  the  following:  For  a  prescribed 
value  of  the  dead-load  traction  po  >  0,  we  seek  a  pressure  field  p(R)  and  a  constant  c>0 
such  that  (4.2.22)  -  (4.2.24)  are  satisfied  where  Tn.,  Tqq  are  given  by  (4.2.16),  (4.2.17), 
(4.2.3).  In  addition,  if  c  >  0,  then  the  condition  for  a  traction-free  cavity  surface 
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must  also  hold. 


Tff  (0)  =  0 


(4.2.25) 


4.3.  Solutions 

On  using  the  normalization  condition  (4.2.6),  it  may  readily  be  shown  that  one 
solution  of  the  foregoing  problem  for  all  values  of  po  is 


p(R)  =  2[W1(3,3,  1,0,  1;  R)  -  W2(3,  3,  1,0,  l;R)]-p0,  c  =  0.  (4.3.1) 

This  corresponds  to  the  trivial  homogeneous  state  of  deformation 

r(R)  =  R ,  (4.3.2) 

with  corresponding  stresses  Tff  =  Tee  =  =  p0,  so  that  the  sphere  remains  undeformed 

under  a  state  of  hydrostatic  pressure.  Note  that  even  though  the  pressure  p  given  in 
(4.3.1)  is  in  general  discontinuous  across  the  interface  R  =  A,  the  corresponding  stresses 
are  continuous. 


We  next  describe  solutions  for  which  c  >  0,  corresponding  to  the  presence  of  a 
traction-free  cavity  at  the  origin.  We  again  employ  the  notation  used  in  Chapter  3  for 
convenience: 


v  =  v(R)MR)/R  =  (1  +  -T3-) 

R 

-2  • 
v  1  =r 


(4.3.3) 


and  thus  rewrite  the  invariants  (4.2.11)  -  (4.2.15)  as 
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Ii  =  v“4  +  2v2  , 
I2=v4+2v-2, 

I3  =  l, 


I4  =  0 , 
I5=v"4. 


(4.3.4) 


Following  the  techniques  of  Chapter  3,  Section  3,  we  rewrite  the  differential  equation 
(4.2.22)  in  the  form 


_d_ 

dR 


|2  [v“4  Wi  (V4  +  2v2,  v4  +  2v~2,  1,  0,  v"4;  R)  -  v4  W2(v"4  +  2v2,  v4  +  2v"2,  1,  0,  v"4;  R) 
+  v"4  W5  (v-4  +  2v2,  v4  +  2v-2,  1,  0,  v"4;  R)]  -  p(R)  ► 

J 

+  [(v-4  -  v2)  Wi  (v"4  +  2v2,  V4  +  2v~2, 1, 0,  v-4 ;  R)  (4.3.5) 

R 


+  (v~2  -  V4)  W2  (v~4  +  2v2,  v4  +  2v“2, 1, 0,  v"4; R) 

+  V"4  W5(v"4  +  2v2,  v4  +  2v-2,  1,  0,  v_4;R)]  =  0on0<R<  A,  A<R<B  . 


On  integration  of  (4.3.5),  we  have 


P(R)  -  P(0)  =  2 

and 


v-4^)  (Wi  +  W5)l  -  v4(R)  W2  I  r  +  4J(R),  0<R<A,  (4.3.6) 

R  R 


p(R)  -  [p]£  -  p(0)  =  2 i  v-4(R)  (W!  +  W5)  I R  -  v4(R)  W2  I  ] 


-  v-4(A)  [Wi  +W5]a^  +v4(A)  [W2]a^ 


(4.3.7) 
)-  +  4J(R),  A<R<B  , 


where 


J(R)  =  j  [[v"7(s)  -  v-1  (.)]  Wi  I  g  +  [v-5(s)  -  v(s)j  W2 1  g  +  v-7(s)W5  I J  y ,  (4.3.8) 

on  0  <  R  <  B,  and  [f(s)] a-  =  f(A+)  -  f(A-).  In  (4.3.6)  -  (4.3.8)  the  W/  (/  =  1,  2,  5)  are 
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evaluated  at  the  values  of  the  invariants  given  by  (4.3.4),  and  for  simplicity  we  retain  an 
obvious  modification  of  the  notation  introduced  in  Chapter  3,  i.e., 

W,  I  =  W i  (v-4 (s)  +  2v2(s),  v4(s)  +  2v~2  (s),  1,  0,  v"4(s);  s) .  (4.3.9) 

S 

On  substitution  into  (4.2.16),  and  using  (4.3.3),  the  radial  stress  becomes 


Tir(R)  =  -p(0)-4J(R),  0 <R < A , 


Tff(R)  =  -p(0)  -  [p]£  +  2\ v'4(A)  [Wx  +  W5]£  -  v4(A)  [W2]£ 


A+  ,,4/ 


4J(R),  A<R<B. 


(4.3.10) 

(4.3.11) 


The  interface  condition  (4.2.23)  now  shows  that 

-  [p] a-  +  2-  V4(A)  [Wj  +  W5]£  -  v4(A)  [W2]a-  '  =  0,  (4.3.12) 

so  that  (4.3.11)  becomes 

T,r(R)  =  -p(0)  -  4J(R),  A  <  R  <  B  .  (4.3.13) 

The  traction-free  cavity  surface  condition  (4.2.25),  together  with  (4.3.10)  and  J(0)  =  0, 
now  yield  p(0)  =  0,  so  that 

T„(R)  =  -4J(R) ,  0  <  R  <  B  .  (4.3. 14) 

Finally,  by  virtue  of  (4.3.14),  (4.3.3),  we  see  that  the  dead-load  boundary  condition 
(4.2.24)  at  the  outer  surface  R  =  B  is  satisfied  if 

-  4J(B)  =  p0  [v(B)]“2  •  (4.3.15) 

As  in  Chapter  3,  Section  3,  the  condition  (4.3.15)  may  be  written  compactly  on 
utilizing  the  s  -»  v  change  of  variables  in  the  integral  (4.3.8).  Introducing  the  function 


(4.3.16) 


W  (x;  y)  =  WCx-4  +  2x2,  x4  +  2x  2,  1,  0,  x  4;  y) , 
and  adopting  the  notation 

Wi  (x;  y)  =  W(x;  y)  (4.3.17) 

as  in  [20],  it  follows  similarly  from  the  arguments  contained  in  Appendix  A  that  one  may 
write  (4.3.15)  as 


Po=  1  +  ^3 


<1  +  £>* 


Wi  (v;  ,  ,,  ) 

(v3-!)73 

v3  - 1 


dv,  c>0.  (4.3.18) 


We  remark  that  (4.3.18),  here  established  for  the  anisotropic  composite  sphere,  has 
formally  the  same  structure  as  the  analogous  relation  for  the  isotropic  composite  sphere 
(see  equation  (2.23)  of  [20]).  As  discussed  in  [20],  the  integrand  in  this  equation  is,  in 

c3  y 

general,  discontinuous  at  the  interface  v  =  (l-l — =-)  •  Equation  (4.3.18)  was  first 

AJ 

established  in  [5]  for  the  case  of  the  homogeneous  isotropic  sphere  (see  equation  (5.18) 
of  [5]).  Thus,  for  a  given  dead-load  po,  solutions  involving  a  traction-free  internal  cavity 
of  radius  c  exist  provided  that  c  is  a  positive  root  of  (4.3.18).  The  associated  pressure 
field  is  given  by 

P(R)  =  2  [v"4(R)  (Wx  +  W5)  I  r  -  v4(R)  W2 1  r]  +  4J(R)  (4.3.19) 

onO<R<A,  A<R<B. 

In  summary,  we  have  seen  that  for  all  values  of  the  applied  dead-load  traction  p0, 
one  obtains  the  trivial  solution  (4.3.1)  corresponding  to  the  homogeneous  state  of 
deformation  (4.3.2).  Moreover,  if  positive  roots  c  of  (4.3.18)  exist,  then  one  obtains  the 
additional  solutions  involving  a  traction-free  internal  cavity  described  above. 


88 


4.4.  The  critical  load  and  bifurcation 

Henceforth  we  confine  attention  to  a  composite  sphere  composed  of  two 
homogeneous  anisotropic  phases.  To  remove  the  explicit  radial  dependence  from  the 
strain-energy,  we  now  write  W  =  W(  Ij ,  I2, 13 , 14, 15 ;  R)  for  the  composite  sphere  as 

W  =  W1(I1,I2,I3,I4,l5)  ,  0<R< A 
W  =  W2(I1,I2, 13, 14, 15)  ,  A<R<B.  (  '  ' 

Using  the  definitions  in  (4.3.16),  (4.3.17),  and  (4.4.1),  we  rewrite  (4.3.18)  as 


-1  (1  +  T5-),/4  -2 

~  W,  (v)  ,  f  W1(v)J 

(3  s  ^rrdV'C>0  (4.4.2) 

(i  +  4r)/4  (1  +  -^-)m 


(i  +  Ar)* 


The  critical  load  pcr  at  which  an  internal  cavity  may  be  initiated  is  then  found  by 
formally  letting  c  -»  0+  in  (4.4.2),  and  so 


wl(v) 

v3  —  1 


dv . 


(4.4.3) 


As  seen  in  Chapter  3,  pCT  may  or  may  not  be  finite,  and  thus  cavitation  may  or  may  not 

take  place  since  the  integral  (4.4.3)  is  improper.  We  see  from  (4.4.3)  that  the  finiteness 

of  per  and  its  value  depend  only  upon  the  material  at  the  core  of  the  composite  sphere. 

Such  a  result  was  also  obtained  in  [21]  for  a  composite  sphere  with  isotropic  phases  (see 

also  [34]).  In  Section  4.5,  we  derive  an  alternative  formula  for  the  critical  load  involving 

*  1 

the  strain-energy  function  W  (v)  instead  of  its  derivative,  as  in  (4.4.3). 

As  discussed  in  Chapter  3,  Section  4,  and  shown  in  Chapter  2,  the  integrand  in 

*  1 

(4.4.3)  is  finite  as  v  — >  1.  Thus,  as  before,  the  behavior  of  W  (v)  for  large  values  of  the 
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stretch  v  determines  whether  cavitation  will  or  will  not  take  place.  That  is,  restrictions 

^  1 

on  the  rate  of  growth  of  W  (v)  for  large  stretch,  similar  to  those  described  in  Chapter  3, 
Section  4,  are  required. 

We  now  consider  the  local  character  of  the  bifurcation  at  po  =  pcr  by  analyzing  the 
curve  po  =  Po(c)>  given  by  (4.3.18),  for  small  values  of  c.  We  first  introduce  the 


dimensionless  parameters 


c  B  , 

p  =  —  ,  a=  — ,  (col) 


and  define 


(4.4.4) 


Po(c)  =  Po(P' 


w!  (v)  .  W?(v) 


,2/  r  «l  v;  .  r  ”1 

)  =  (l  +  p3f  f  -1— -dv+  J  -3— 

(1  +  aV)*  ^_1  d  +  P3)M  ^-1 


(4.4.5) 


A  Taylor  expansion  of  (4.4.5)  about  p3  =  0  shows  that 


Po  =  Per  +  k  P3  +  o  (p3)  as  p3  ->  0 


where 


From  (4.4.5)  one  obtains 


* /A*  3 


k=  lim  dp0/dp 


p3  -»o 


(4.4.6) 


(4.4.7) 


dW2/dv 


,  2  1  3 dW  /dv  ,  1  /  3  dW  /& 

k  =  —  pcr  -  —  or  lim  — 5 - +  —  (a  -  1)  lim  — ; — 7 

3  3  v->  1  v3-l  3  v->i  v3-! 


(4.4.8) 


On  using  l’Hopital’s  rule  and  defining 


f  =  a-3  ,  0  <  f  <  1 


(4.4.9) 


as  the  volume  fraction  of  the  core  material  to  the  total  material,  we  write  k  as 


,  _  2  n  1  !  d2WX(l)  1  !  d2W  (1) 

k-3  P“‘6r  6(r  _I)  dv2 


(4.4.10) 
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where,  from  the  considerations  of  Chapter  3,  Section  4,  and  Chapter  2,  the  quantities 
d2W  (l)/dv2,  (i=l,2),  are  positive  and  finite,  and  each  can  be  written  in  terms  of  the 
associated  infinitesimal  Young’s  moduli  in  the  radial  direction.  We  see  from  (4.4.6)  that 
if  k  >  0,  bifurcation  is  locally  to  the  right  (supercritical)  while  if  k  <  0,  bifurcation  is 
locally  to  the  left  (subcritical).  Equations  corresponding  to  (4.4.6),  (4.4.10)  were  also 
determined  in  [21]  for  the  case  of  an  isotropic  composite  sphere,  as  well  as  in  Chapter  3 
for  the  case  of  a  homogeneous  transversely  isotropic  sphere.  (To  recover  the  latter,  set 
the  volume  fraction  f  =  1  in  (4.4.10)  and  compare  (4.4.6),  (4.4.10)  with  (3.4.7),  (3.4.10).) 
We  now  make  the  further  definitions 


so  that 


T,  k  d2W2(l)/dv2 

K  =  — r^rr - -  ,  CO  —  -~i  - -  , 

d2W  (1)  /  dv2  d2W  (l)/dv2 


a  _  Per  /v  __  P0 

d2w\l)/  dv2  d2w\l)  /  dv 


2  ’ 


fc-jr1  +  -g-(r'-i)o) 


(4.4.11) 


(4.4.12) 


Letting  f  1  in  (4.4.12),  we  obtain  the  value  of  K  appropriate  to  the  homogeneous 
transversely  isotropic  sphere  composed  of  the  core  material,  namely 

Kh«|(p„-{).  (4.4.13) 

In  contrast  to  the  case  for  a  homogeneous  isotropic  sphere  for  which  Kh  >  0  so  that 
bifurcation  is  only  to  the  right,  (see  §  5.2  of  [5]  or  eqn.  (25)  of  [21] ),  we  saw  in  Chapter 
3  that  for  the  anisotropic  sphere,  bifurcation  can  occur  either  to  the  right  or  to  the  left. 
Thus,  for  the  anisotropic  sphere,  Kh  can  be  either  positive  or  negative. 
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Returning  to  (4.4.12),  we  now  discuss  the  circumstances  which  determine  whether 
the  bifurcation  for  the  composite  sphere  is  locally  to  the  right  or  to  the  left  It  is 
convenient  to  treat  these  cases  with  reference  to  the  schematic  diagrams  in  Figures  4.1 
and  4.2.  Figure  4.1  corresponds  to  the  case  in  which  bifurcation  for  the  homogeneous 
sphere  composed  of  the  core  material  would  be  to  the  right,  and  Figure  4.2  to  the  case  in 
which  such  a  sphere  would  bifurcate  to  the  left.  In  these  figures,  the  semi-infinite  strips 
0<f<l,0<co<oo,are  divided  into  two  regions,  I,  II  by  the  curve 


f(co)  = 


to-  1 

C0-6pcr 


(4.4.14) 


on  which  K  =  0.  When  Kh  >  0  (Fig.  4.1),  the  curve  (4.4.14)  is  monotone  decreasing  on 
the  region  of  interest,  and  when  Kh  <  0  (Fig.  4.2),  the  curve  (4.4.14)  is  monotone 
increasing  with  f(co)  —>  1  as  to  — »  °°  .  In  view  of  the  remark  made  above  following 
(4.4.13),  we  note  that  there  is  no  diagram  analogous  to  Fig.  4.2  for  the  composite 
isotropic  sphere  considered  in  [20, 21].  In  both  Figures  4.1  and  4.2,  if  (to,f)  lies  in  region 
I,  then  K  >  0  and  so  bifurcation  for  the  composite  sphere  is  to  the  right,  while  if  (co,f) 
lies  in  region  n,  then  K  <  0  and  so  bifurcation  is  to  the  left. 

As  in  [21],  we  may  discuss  these  diagrams  with  reference  to  a  composite  sphere  of 
either  fixed  geometric  properties  together  with  varying  material  properties  or  vice  versa. 
We  first  consider  Figure  4.1  so  that  Kh  >  0  ( =>  pcr  >  1/6 ).  Taking  the  former  view,  we 
see  that  for 


f> 


1 

6pCT 


(4.4.15) 


bifurcation  is  always  to  the  right.  Thus,  if  the  volume  fraction  of  the  core  to  the  total 
volume  exceeds  l/6pCT,  bifurcation  is  locally  to  the  right  regardless  of  the  material 
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properties  of  the  outer  material.  Taking  the  latter  view,  we  also  see  that  for  CO  >  1,  i.e., 
for 


d2W2(l)  >  d2^!) 
dv2  dv2 


(4.4.16) 


bifurcation  is  always  to  the  right.  It  is  shown  in  Chapter  2  that  for  the  transverse  isotropy 
of  concern,  one  has 


d2W(l) 

dv2 


=  4E 


(4.4.17) 


where  E  is  Young’s  modulus  in  the  direction  of  anisotropy,  i.e.,  the  radial  direction, 
associated  with  infinitesimal  deformations  of  the  transversely  isotropic  material  at  hand. 
We  thus  make  the  following  definition  for  the  composite  sphere, 


Ex :  E  for  core  material, 

E2  :  E  for  surrounding  material. 


(4.4.18) 


Thus,  from  (4.4.16)  -  (4.4.18) ,  if  the  material  in  the  surrounding  shell  is  stiffer  in  tension 
in  the  radial  direction  than  the  core  material,  (  E2  ^  Ej  ),  bifurcation  is  to  the  right 
irrespective  of  geometry. 

On  the  other  hand,  the  condition  for  bifurcation  to  the  left  requires  both  geometric 
and  material  restrictions.  These  may  be  written  as 

„  CD-I 

f  < - 77-  (4.4.19) 

CO  6pcr 

and 

(D  <  1 .  (4.4.20) 


Thus,  the  occurrence  of  bifurcation  to  the  left  in  a  composite  sphere  corresponding  to 
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Figure  4.1  requires  a  sufficiently  small  core  surrounded  by  a  shell  of  sufficiently 
compliant  material. 

We  next  consider  similar  implications  of  Figure  4.2,  so  that  now  the  core  material 
alone  would  undergo  bifurcation  to  the  left  (Kjj  <  0  =>  pCT  <  1/6).  In  contrast  to  the  case 
of  Fig.  4.1,  we  cannot  determine  a  single  condition  for  bifurcation  to  the  right  to  occur. 
Thus,  we  must  have 

r  co- 1 

f< - (4.4.21) 

co  6pcr 

and 

co  >  1 .  (4.4.22) 

In  other  words,  the  occurrence  of  bifurcation  to  the  right  in  a  composite  sphere 
corresponding  to  Fig.  4.2  requires  a  sufficiently  small  core  surrounded  by  a  shell  of 
sufficiently  stiff  material. 

The  situation  for  bifurcation  to  the  left  is  much  different.  Although  we  can  draw  no 
conclusions  based  on  geometric  properties  alone,  we  can  deduce  that  for 

co<l  (=>E2<E!  ),  (4.4.23) 

bifurcation  is  always  to  the  left.  Thus,  if  the  material  in  the  surrounding  shell  is  more 
compliant  in  tension  in  the  radial  direction  than  the  core  material,  bifurcation  is  to  the 
left  irrespective  of  geometry. 

Finally  from  Fig.  4.2,  if 

co  >  1  ( =>  E2  >  Ei  )  (4.4.24) 

and 
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,  CO-  1 

f> - —  >  (4.4.25) 

co  -  6pcr 

we  also  have  a  bifurcation  to  the  left  occuring.  Thus,  bifurcation  to  the  left  in  a 
composite  sphere  corresponding  to  Figure  4.2  is  also  possible  provided  there  is  a 
sufficiently  large  core  surrounded  by  a  shell  of  sufficiently  stiff  material 

We  conclude  this  section  by  again  emphasizing  that  the  behavior  of  the  composite 
sphere  is  highly  dependent  upon  the  core  material.  Whether  cavitation  will  or  will  not 
take  place  depends  solely  on  the  core  material.  The  value  of  the  critical  load  depends 
only  on  the  core  material.  Furthermore,  material  properties  of  the  core  alone  determine 
which  of  the  schematic  diagrams  (Fig.  4.1  or  Fig.  4.2)  applies. 

4.5.  Energy  and  stability  of  solutions 

To  examine  the  stability  of  the  cavitation  solutions,  we  now  carry  out  an  energy 
analysis  for  general  transversely  isotropic  strain-energy,  W.  As  in  Chapter  3,  Section  6, 
the  total  energy  associated  with  any  radially  symmetric  equilibrium  configuration  of  the 
body  is  given  by 

B 

E(c)  =  4rc  J  W  (? 2  +  2r2/R2,  r4/R4  +  2r  i^/R2,  1,  0,  r  2;  R)  R2  dR 

°  (4.5.1) 

-  4jiB2  po  [r(B)  -  B] . 

For  the  trivial  solution  (4.3.1)  with  deformation  (4.3.2),  we  have 

E(0)  =  0  (4.5.2) 


on  using  the  normalization  (4.2.4).  Following  the  approach  taken  in  Chapter  3,  Section  6, 
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we  now  consider  (4.5.1)  to  hold  for  all  values  of  c  >  0,  and  thus  seek  values  of  c  >  0  that, 
for  a  given  po  >  0,  minimize  the  extended  energy  function  E(c). 

Recalling  the  notation  (4.3.16),  (4.3.3),  we  may  write  equation  (4.5.1)  as 


v2  (v3  -  I)-2  W(v; 


E(c)— 47tC3  f  »  v*  a  \a 

(l+c3/B3)w  (v3-!) 


)  dv 


(4.5.3) 


4jtB2p0  [r(B)  -  B],  c  >  0 . 


Henceforth,  for  simplicity  of  notation,  we  omit  the  explicit  dependence  of  W  on  its 
second  argument  Upon  introduction  of  the  dimensionless  quantity 

^  E 

2= - ,  - T  ’  (4.5.4) 

(4/3)  TtB3  d2W  (1)  /  dv2 

and  using  (4.4.1),  (4.4.4),  and  (4.2.3),  we  rewrite  (4.5.3)  as 


d2W  (1) 
dv2 


(1  +  a3p3)w 


X(p)  =  3p3  <  J  V2  (v3  -  l)"2  W  (v)  dv 


(i  +  P3)* 


+  J  v2  (v3  -  1)  2  WJ(v)  dv 


(1  +  a3p3)M 


(4.5.5) 


h-  3Po  [(i  +  p3)1/s  ~  1]  • 


We  now  let 


^(P)  =  2  m  1 - -[^>(p)  +  p0 'P(p)],  p>0  (4.5.6) 

d2W  (1)  /  dv2 


where 
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$(p)  =  3p3 


(l  +  a3p3)M  2  \ir2/  \  00  2bVm 

f  ilM  dv  +  f  y.  w  (v)  dv 

(1  Jp3)-  (V3-1)2  (1  +  «W 


(4.5.7) 


and 


T(p)  =  -3[(l+p3)'-4-l] 


(4.5.8) 


We  thus  state  the  energy  minimization  problem  as:  for  fixed  po  >  0,  we  seek  those  values 
ofp  >  0  that  minimize  X(p). 

Proceeding  as  in  Chapter  3,  we  first  examine  the  end-point  p  =  0  where  X(0)  =  0.  It 
is  shown  in  Appendix  D  that 

t  (0+)  =  2  (0+)  =  0,  2  m  =  6(pCT  -  p) ,  (4.5.9) 

where  the  superposed  dot  again  denotes  differentiation  with  respect  to  p,  and  pcr ,  p  are 
given  by  (4.4.1 1),  (4.4.3).  Thus  p  =  0  is  a  local  minimum  if  p  <  pcr  (p0  <  pcr)  and  a  local 
maximum  if  p  >  pcr  (Po  >  Per)-  Observe  that 

p  =  0 ,  po  >0,  (4.5.10) 

corresponds  to  the  trivial  solution  (4.3.1). 

We  now  consider  the  possible  existence  of  internal  extrema  and  thus  seek  values  of 
p  >  0  for  which 

t(p)  =  0  and  so  «f>(p)  +  p0 'F(p)  =  0.  (4.5.11) 

From  (4.5.8),  'P(p)  <  0  for  p  >  0  so  that  (4.5.1 1)  may  be  written  as 

Po=-4^,  p  >0 .  (4.5.12) 

^P) 

We  again  note  that  the  expression  (4.5.12)  for  po(p)  derived  from  energy  considerations 
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can  be  shown  to  be  identical  to  the  expression  (4.3.18)  determined  directly  in  Section  4.3. 
Thus,  as  expected,  the  extremals  provided  by  (4.5.12)  correspond  to  equilibrium 
solutions  with  an  internal  cavity. 

We  next  analyze  the  dimensionless  energy  £(P)  given  by  (4.5.6)  -  (4.5.8)  near 
p3  =  0.  A  straightforward  calculation  shows  that 

¥(p)  =  Y*(p3)  =  -  p3  +  yP6  +  0(p9)  as  p3  0 .  (4.5.13) 


On  using  (4.5.13)  and  (4.4.6),  the  last  term  in  (4.5.6)  has  the  asymptotic  expression 


Po(pmp)  =  -  PcrP3  +  ( jPcr  -  k)  P6  +  0(P9)  as  P3  ->  0  , 


(4.5.14) 


where  k  is  given  by  (4.4.10).  We  next  consider  the  Taylor  expansion 


«*>(p)  =  ®V)  =  <E>'(0)  +  ^j-(0)p"  +  “^2  (°)p°  +  O(P  )  as  p  ->  0 


dp 


2  d(p3)2 


It  is  easily  seen  that 


(4.5.15) 


O  (0)  =  0. 


(4.5.16) 


This  follows  from  the  fact  that  E(0)  =  0  and  so  from  (4.5.6),  <P  (0)  +  po  (0)  =  0  ; 
however,  (4.5.8)  implies  that  VP  (0)  =  0 ,  so  that  (4.5.16)  must  hold. 

Now,  differentiating  (4.5.7)  we  obtain 


d<f>*  = 
dp3 


(1  +  a3p3)M 
(1  +  P3)W 


v2  W2(v) 
(v3  - 1)2 


dv  + 


(l  +  a3p3)w 


v2  W*(v) 
(v3  -  l)2 


dvr 


v(A)=(1-hx3p3)h 


-W2(v)l 


"I 

v(B)=(1+P3)M  f  ’ 


J 


(4.5.17) 
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so  that 


d<J>*  7  v2  M 

— ^-(0)  =  3  f  -V — T  W  (v)  dv 
dp3  i  (v3  -  l)2 

(4.5.18) 

,  mu  «~3[  W2((1+«3P3)‘'3)-W1((1-hi3p3)''3)1- W2((l+p3)w) 

p3->0  p3 


Integrating  the  first  term  on  the  right  hand  side  by  parts  gives 


f  3v2  (v3  -  1)~2  w\v)  dv  =  lim^-^-  -  lim  +  j  3—  dv  ,  (4.5.19) 

{  v — >1  v3  -  1  v— v3  -  1  1  v3  -  1 


Wi(v) 


where  we  note  from  (4.4.3)  that  the  last  term  in  the  above  equation  is  pCT.  Since 
WX(1)  =  0  (i=l,2)  from  (4.2.5)  and  (4.3.16),  using  L’Hopital’s  rule  and  (3.4.2)i  shows 
that  the  first  term  on  the  right  in  (4.5.19)  is  zero.  Also,  it  follows  from  considerations 
similar  to  those  of  Section  4.4  that  for  a  finite  pCT,  i.e.,  for  cavitation  to  take  place,  we 
must  have 

W1  (v)  =  0(vn)  ,  n  <  3  ,  as  v  ->  ~ .  (4.5.20) 


Thus,  the  second  term  on  the  right  in  (4.5.19)  is  also  zero  so  that  we  arrive  at  the 
following  alternative  formula  for  the  critical  load  (when  the  strain-energy  satisfies  the 
growth  assumption  (4.5.20)): 

3v2  -1  7  w!(v)  J 

^177  w  w = |  777 dv  “i p»  •  <4-5-21) 

Now,  the  limit  in  (4.5.18)  is  zero  by  arguments  similar  to  those  used  following  equation 
(4.5.19).  Thus,  on  using  the  formula  (4.5.21)  for  the  critical  load,  (4.5.18)  reduces  to 


d<D* 

dp3 


(0)  =  Per  • 


(4.5.22) 
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d20>* 


We  next  differentiate  (4.5.17)  and  obtain 

r  a  2  a  i  i 

Wx^-WjCv)] 


d(p3) 


3n2  3  < 


v(A>=(l+a3P3)W 


p  (1+a  p  ) 


~  2 

W 1  (v)  I  v(B)=(l+P3)‘/i 

p3  (i+p3)34 


(4.5.23) 


so  that  taking  the  limit  as  p3  -¥  0  and  applying  L’Hopital’s  rule  gives 


r  r 


d(p  ) 


a 


V.  V. 


d2w2  d2W 

-(I)--Vt- d) 


.2.1 


dv 


dv 


r.  *  2 

d2W 

dv2 


(1) 


(4.5.24) 


Combining  (4.5.15),  (4.5.16),  (4.5.22),  (4.5.24),  and  (4.4.9)  then  gives 


r  r  ^  9 


as  p3  — »  0. 


**,  3  .  1  j  f_i  dzW  dzW 

®  (p  ) = Pcrp3  +  -jg  •(  r1  — r^-(i)  -  — — (i) 


I 


dv2 


dv 


,*2 

d2W 

dv2 


(1)^P6  +  0(p9) 

(4.5.25) 


Finally,  using  (4.5.14)  and  (4.5.25)  in  (4.5.6),  we  obtain 

£(p)  =  I*(p3)  =  -  Vi  Kp6  +  0(p9)  as  p3  0 ,  (4.5.26) 

where  K  is  given  by  (4.4.12).  The  above  result,  which  holds  for  a  composite  sphere 

consisting  of  two  general  homogeneous  transversely  isotropic  phases,  generalizes  those 

found  previously  in  [20],  (see  eq.  (3.17))  for  a  composite  sphere  composed  of  two  neo- 

Hookean  phases,  and  in  Chapter  3,  (see  equation  (3.6.13)),  for  the  homogeneous 

transversely  isotropic  sphere  composed  of  the  specific  material  described  in  Chapter  3, 

Section  5.  Equations  (4.5.21)  and  (4.5.26)  of  course  hold  for  the  general  problem 

»2  ~  1 

considered  in  Chapter  3  (i.e.,  a  =  1 ,  W  =  0).  In  this  case,  W  represents  the  stored- 
energy  function  for  the  homogeneous  transversely  isotropic  sphere,  and  K  appearing  in 
(4.5.26)  is  given  by  (4.4.13),  (see  (4.4.12)  and  the  sentence  following  (4.4.12)). 
Returning  to  the  problem  at  hand  and  recalling  the  results  at  the  end  of  Section  4.4, 
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where  the  curve  (4.5.12)  (or  equivalently  (4.4.2))  was  shown  to  bifurcate  from  the  trivial 
solution  (4.5.10)  at  the  point  p  =  0,  p0  =  pCT  ,  (4.5.26)  shows  that  X  <  0  on  this  curve 
near  p  =  0  if  the  bifurcation  is  to  the  right  (K  >  0 )  and  X  >  0  if  the  bifurcation  is  to  the 
left  (K  <  0).  We  have  thus  established,  in  general,  all  of  the  energy  results  which  were 
essential  for  the  specific  material  models  in  [20]  (and  in  Chapter  3,  Section  6)  to  establish 
Properties  A  and  B,  which  we  now  restate  for  convenience: 

Property  A:  The  curve  p  =  p(p)  is  a  locus  of  local  minima  of  X(p)  for  fixed  p 
whenever  dp/dp  >  0,  and  is  a  locus  of  local  maxima  of  X(p)  for  fixed  p  whenever 
dp/dp  <  0. 

Property  B:  On  the  curve  p  =  p(p),  the  values  of  X(p)  are  decreasing  with  p 
whenever  dp/dp  >  0  and  are  increasing  with  p  whenever  dp/dp  <  0. 

With  the  above  results  we  have  that,  in  general,  where  the  curve  p  =  p(p)  is 
monotone  decreasing,  the  cavitated  solution  is  unstable,  and  where  p  =  p(p)  is  monotone 
increasing,  the  cavitated  solution  is  stable.  If  there  exist  two  stable  solutions 
corresponding  to  a  particular  load,  we  again  say  that  a  solution  is  absolutely  stable  (AS) 
if  it  provides  the  absolute  minimum  for  X  and  is  metastable  (MS)  if  it  provides  a  local 
minimum  for  X  which  is  not  the  absolute  minimum.  In  cases  such  as  this,  there  may  be  a 
discontinuous  transition  in  cavity  radius  in  moving  from  one  stable  solution  to  another. 
The  transition  between  metastability  and  absolute  stability  will  be  addressed  more 
completely  in  the  illustrative  example  considered  in  the  following  section. 


4.6.  Illustrative  example 


In  this  section,  we  illustrate  the  foregoing  results  for  a  specific  material  model. 
Consider  a  composite  sphere  whose  phases  are  described  by  the  transversely  isotropic 
strain-energy  density  function  discussed  in  Section  3.5.1.  Explicitly,  we  have 

F 

-^-[(I1-3)  +  aias-2I5  +  l)]  ,  0<R< A 

|i2  [  2  ]  (4.6.1) 

y  [(I1-3)  +  a2(li-2I5  +  l)J  ,  A<R<B 

where  Ii  ,  I5  are  given  by  (4.2.11),  (4.2.15),  the  m  (i  =  1,2)  are  positive  constants,  and 
the  with  0  <  ai  <  5.0  (i=l,2)  are  dimensionless  parameters  which  serve  as  measures  of 
the  degree  of  anisotropy  in  the  respective  phases.  When  ai  =  0 ,  =  0,  one  recovers  the 

neo-Hookean  composite  sphere  treated  in  [20,  21].  Recalling  the  notation  (4.3.3), 
(4.3.16),  and  using  (4.3.4),  we  rewrite  (4.6.1)  as 

f 

W1  (v)  =  y  [(v-4  +  2v2  -  3)  +  ai  (v~8  -  2V-4  +  1)]  ,  0  <  R  <  A 

W(v)=<  .  (4.6.2) 

W  (v)  =  y  [(v-4  +  2v2  - 3)  +  a2  (v-8  - 2v“4  +  1)]  ,  A<R<B 

I. 

where,  from  the  considerations  of  Chapter  2, 

d2Wl(l)/dv2  =  4Ei  =  4ni(  3  +  4aj )  >  0  (i  =  1,2)  (4.6.3) 

for  this  material.  We  recall  that  the  Ej  are  Young’s  moduli  in  the  radial  direction. 

On  substitution  from  (4.6.2)  into  (4.5.7)  and  from  (4.6.3)  with  i  =  1  into  (4.5.6),  and 


making  the  definitions 


Hi 

we  find  that  (4.5.6)  -  (4.5.8)  and  the  considerations  of  Appendix  C  give  the  total  energy 
for  the  composite  sphere  as 


2(p)  = 


1 


4(3  +  4ai) 


*(£i  +  *L  y(p) 

Hi  Hi 


(4.6.5) 


where 


$(p)  =  3  I 
Hi  2 


-a-3  (3-1)  [  2z2(ap)  -  z"1  (ap)  - 1]  +  p  [  2z2(p)  -  z"1  (p)  -  1 


-(pa2  -ax) 


Via  3z  5(ap)  +  -r-p 


8.3 


-JZ  5(ap)  +  yz  2(ap) 


1  .  V3z(ap)  1  .  2z(ap)  + 1 

—=r  arctan- - - —  arctan — - 

V3  2  +  z(ap)  V3 


+  a  3(  2z  x(«p) ~  -r-zz(ap)  +  Vi ) 


8  2/ 


+  Pa2 


Viz-5(p)  +  |p3 


■~z  5(p)  +  ~z  2(p)  — p  arctan  ^  — p  arctan  2z^P^_+  1 

5  VK'  2  VF/  y/3  2  +  z(p)  ^3  V3 


+  2z"1(p)--|z2(p)  +  1/3 


20jt  3 
+  -^ra!p3> 


9  V3 


(4.6.6) 


and  'F(p)  is  given  by  (4.5.8).  Furthermore,  from  either  (4.5.12)  or  (4.4.5),  and  making 
the  definitions 


F[  z(x)]  =  z"1(x)  +  1/4z-4(x), 

G[  z(x)  ]  =  F[  z(x) ,  -  ■/,  z-2(x)  -  jz-5(x)  -  M  z-8(x) 


.-4/ 


V3 


2z(x)  + 1  .  .  ^  z(x) 

arctan — —  +  arctan- - 7-7- 

V3  2  +  z(x) 


(4.6.7) 


we  obtain  the  corresponding  load/cavity  relationship  as 
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P0(P)  „  nJ2 


Pi 


2zz(p) 


^  1  -  pj  F  [  z(ap)  ]  +  P  F  [  z(p)  ] 

+  2[ Pa2  -  ai]  G  [  z(ap)  ]  -  2Pa2G  [  z(p)  ]  + 


(4.6.8) 


For  small  values  of  p  we  have  (see  (4.4.6),  (4.5.26)  respectively) 

P(P)  =  Per  +  Kp3  +  o(p3)  as  p3  ->  0 , 
S(p)  =  -Vz  Kp6  +  o(p6)  as  p3  0 , 

where,  for  the  material  at  hand, 


and 


-  1 
Pcr  4(3  +  4a0 


5 

40  71-51^3 

2  + 

30 <43 

L  J 

ai 


so  that 


Per 

Pi 


5  , 

40tc-51V3 

2  + 

30V3" 

ai 


K  = 


2 

3 


Pc-jr1  +  |(r>-i)p 


(3  +  4a2) 
(3  +  4a,) 


(4.6.9) 

(4.6.10) 


(4.6.11) 


(4.6.12) 


We  now  examine  the  curve  (4.6.8)  and  corresponding  energy  for  the  material  (4.6.1) 
for  large  values  of  p.  A  detailed  analysis  of  (4.6.8)  similar  to  that  discussed  in  Appendix 
B  for  the  homogeneous  sphere,  shows  that  for  large  values  of  p, 


■5^  ~  2  [  P(1  -  a-1)  +  a-1  ]  p  as  p  — »  eo ,  (4.6.13) 

Pi 

while  the  expression  for  the  energy  on  this  curve  for  large  p  follows  from  (4.6.5),  (4.6.6), 
(4.6.13),  and  (4.5.8)  as 
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Mei  — 3  [  p(i  _  a-i)  +  a”1  ]  P2  as  p  co .  (4.6.14) 

Pi 


In  (4.6.14),  £(p)  /  (J-i  is  defined  to  be  the  quantity  in  brackets  in  (4.6.5).  The  expressions 
(4.6.13),  (4.6.14)  are  identical  to  those  found  in  [20]  for  the  composite  neo-Hookean 
sphere.  (See  eqs.  (3.18),  (3.19);  the  additional  factor  of  two  in  the  equations  in  [20]  is 
due  to  the  authors’  use  of  a  shear  modulus  one-half  that  of  the  actual  infinitesimal  shear 
modulus).  Thus,  the  anisotropy  considered  here  does  not  affect  the  large  p  asymptotic 
behavior  of  the  bifurcation  curve,  nor  the  energy  on  it.  Since  a  >  1 ,  P  >  0  we  see  that 
P(1  -  a-1)  +  a-1  >0.  Thus,  the  curve  (4.6.8)  is  asymptotic  to  a  straight  line  with 
positive  slope,  while  the  energy  will  always  be  negative  on  this  curve  for  sufficiently 
large  values  of  p.  Thus,  as  seen  in  [20]  and  in  Chapter  3,  the  initial  character  of  the 
bifurcation  does  not  affect  the  large  p  asymptotic  behavior.  Note  that  we  recover  the 
results  (3.5.30),  (3.6.14)  for  the  homogeneous  anisotropic  sphere  on  setting  a  =  1  in 
(4.6.13),  (4.6.14)  respectively. 

In  Figures  4.3  -  4.6  we  have  plotted  the  curve  (4.6.8)  for  four  different  pairs  of 
values  of  the  anisotropy  coefficients  ai  ,  a2  ,  and  the  parameter  co  defined  in  (4.4.11). 
(Note  that  for  the  material  (4.6.1),  we  have 


co  =  p 


(3  +  4a2) 
<3  +  4aj) 


(4.6.15) 


so  that  choosing  any  three  of  the  parameters  in  (4.6.15)  specifies  the  last.)  All  graphs  in 
Sections  4.6  and  4.7  correspond  to  a  representative  value  f  =  0.1  for  the  volume  fraction 
of  the  core  material  to  the  total  material.  The  situations  corresponding  to  Figures  4.3  and 
4.4  were  discussed  extensively  in  Chapter  3.  We  will  treat  these  briefly  in  what  follows, 
and  pay  closer  attention  to  the  more  novel  phenomena  displayed  in  Figures  4.5  and  4.6. 
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In  the  case  of  Fig.  4.3,  the  composite  sphere  undergoes  a  bifurcation  to  the  right  and  a 
smooth  cavitation  occurs.  Here,  an  internal  cavity  is  nucleated  at  the  critical  load,  with 
its  radius  increasing  continuously  from  zero  as  po  is  further  increased.  In  Fig.  4.4,  the 
composite  sphere  undergoes  a  bifurcation  to  the  left.  There  exists  a  turning  point 
( pt ,  pt )  of  the  curve  at  which 


dpo 

dp 


=  0, 


(4.6.16) 


and  a  cavity  of  finite  radius  is  nucleated  at  a  transition  load  p,  with  pt  <  p  <  Per-  The 
existence  of  a  unique  value  p  =  p  was  established  in  [20]  for  the  composite  neo-Hookean 
sphere,  and,  as  discussed  in  Chapter  3,  follows  from  the  Properties  A  and  B  given  in 
Chapter  3  and  at  the  end  of  Section  4.5.  (The  nondimensional  analog  of  p  was  called  PN 
in  Chapter  3).  Stability  arguments  of  the  type  carried  out  in  Chapter  3  again  show  that  a 
"snap  cavitation"  occurs  in  which  a  discontinuous  change  in  stable  equilibrium 
configurations  takes  place. 


We  consider  now  the  situation  in  which  the  curve  (4.6.8)  bifurcates  to  the  right  from 
the  trivial  solution  (4.5.10)  at  the  point  ( Per ,  0 ),  and  has  two  turning  points, 
( Pti  ,  Pti  )  >  ( Pt2  >  Pt2  )>  at  which  (4.6.16)  holds.  This  is  the  case  in  Figures  4.5(a)  and 
4.6(a).  In  each  of  these  figures,  the  curve  (4.6.8)  is  mononotically  increasing  for 
0  <  p  <  pt2  and  p  >  pti,  and  is  monotonically  decreasing  for  pt2  <  P  <  Pti-  Note  that  in 
Figure  4.5(a),  pti  and  p^  both  fall  to  the  right  of  pCT,  while  in  Figure  4.6(a),  pcr  falls 
between  pti  and  pt2. 


We  first  consider  Figure  4.5(a),  for  which  the  number  of  equilibrium  solutions 
depends  on  the  values  of  po  in  the  following  way:  for  0  <  po  <  pcr,  the  only  solution  is 
the  trivial  solution  p  =  0;  for  pCT  <  Po  <  Pti  >  there  are  exactly  two  solutions,  namely,  the 
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trivial  solution  p  =  0  and  an  additional  bifurcated  solution  involving  an  internal 
cavity;  for  ptl  <  po  <  Pt2>  there  are  exactly  four  solutions,  namely,  the  trivial  solution  p 
=  0  and  in  addition  three  bifurcated  solutions  each  involving  an  internal  cavity;  finally, 
for  Po  >  Pt2»  there  are  exactly  two  solutions,  namely  the  trivial  solution  and  a  single 
bifurcated  solution  involving  an  internal  cavity.  Again,  the  trivial  solution  provides  a 
local  minimum  for  £  and  hence  is  stable  for  0  <  po  <  pcr,  while  for  p0  >  pCT  the  trivial 
solution  provides  a  local  maximum  for  £  and  thus  is  unstable.  By  virtue  of  Property  A, 
the  bifurcated  curve  (4.6.8)  again  provides  a  local  maximum  for  £  if  p^  <  p  <  pti  and 
hence  the  associated  solution  is  unstable,  while  for  0  <  p  <  pt2  and  p  >  pti ,  this  curve 
provides  a  local  minimum  for  £  and  hence  this  solution  is  stable.  Thus  if  0  <  po  <  Pti  or 
if  po  >  Pt2>  there  is  exactly  one  stable  equilibrium  solution,  while  if  ptl  <  po  <  Pt2  there 
are  two  stable  equilibrium  solutions,  each  involving  an  internal  cavity. 

We  now  address  the  transition  between  metastability  and  absolute  stability  in  the 
region  ptl  <  p0  ^  Pt2-  We  again  refer  to  Fig.  4.5(a).  Observe  that  the  stable  branches  on 
the  region  of  interest  correspond  to  the  intervals  ptl  <  p  <  p^  and  ptl  <  p  <  p^,  where 
Pti  >  Pt2  arc  the  points  at  which  Po(Pti)  =  Pti  >  Po(Pt2)  =  Pt2  >  respectively,  and  where 
(4.6.16)  is  not  satisfied.  From  Property  B  and  £(0)  =  0,  we  have  that  £(p)  decreases  from 
zero  on  the  interval  0  <  p  <pt2-  Property  B  and  (4.6.14)  give  that  £(p)  is  decreasing  with 
p  for  p  >  ptl ,  and  £(p)  is  negative  for  sufficiently  large  p.  (Note  that  £(ptl )  can  be  either 
positive  or  negative.)  We  now  define  a  transition  load  po  =  Ptr,  satisfying  ptl  <  p^  <  pt2, 
such  that  for  ptl  <  p0  ^  p^,  we  have  £(p)  on  ptl  <  p  <  p^  less  than  £(p)  on  ptl  <  p  <  pte, 
and  for  p^  <  p0  ^  pt2>  we  have  £(p)  on  p^  <  p  <  p^  greater  than  £(p)  on  ptt  <  p  <  pt2, 
where  p*  ,  pfr  (p^  >  p^)  are  the  values  of  p  occurring  on  the  stable  branches  for  which 
Po(p)  =  ptr-  In  other  words,  Pu-  is  the  value  of  the  dead-load  at  which  the  energy  on  the 
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upper  portion  of  the  bifurcation  curve  becomes  less  than  that  on  the  lower  portion  of  this 
curve,  and  thus  the  value  after  which  the  upper  branch  remains  the  absolutely  stable 
solution.  The  preceding  is  sketched  in  Fig.  4.5(b).  Note  that  if  Ptr  =  Pti  »  then  ptt  =  ptl 
and  the  upper  branch  is  absolutely  stable  for  p  >  Pti,  and  if  Ptr  =  Pt2>  then  pte  =  Pa  and 
the  lower  branch  is  the  absolutely  stable  solution  for  0  <  p  <  pt2- 

Consider  now  the  quasi-static  loading  process  in  which  the  composite  sphere  is 
subjected  to  a  tensile  dead-load  surface  traction  po  which  increases  slowly  from  zero. 
We  again  refer  to  Fig.  4.5(b).  The  trivial  homogeneous  deformation  (4.3.2)  persists  as  po 
is  increased  from  zero,  until  the  critical  load,  po  =  Per*  is  reached.  Then  an  internal  void 
is  nucleated,  with  cavity  radius  increasing  continuously  from  zero  as  po  increases,  until 
Po  =  Ptr-  At  this  point  there  are  two  possibilities.  First,  if  one  supposes  that  the  body 
always  assumes  an  absolutely  stable  configuration,  then  as  po  passes  through  p^,  a 
discontinuous  growth  in  cavity  radius  from  c  =  Bpt,.  to  c  =  Bpfr  is  predicted,  with  cavity 
radius  increasing  continuously  as  po  is  further  increased.  On  the  other  hand,  if  one 
supposes  that  continuous  deformation  of  the  body  persists  as  long  as  possible,  then  the 
cavity  radius  grows  continuously  as  po  is  increased  until  po  =  Pt2  and  then  experiences  a 
jump  from  c  =  Bpt2  to  c  =  Bpt2  as  po  passes  through  pt2.  The  cavity  radius  then 
increases  continuously  as  p0  is  further  increased.  Note  that  the  foregoing  two 
possibilities  coincide  if  Ptr  =  Pt2- 

We  next  consider  Figure  4.6(a),  for  which  the  number  of  equilibrium  solutions 
depends  on  the  values  of  po  in  the  following  way:  for  0  <  Po  <  Pti  >  the  only  solution  is 
the  trivial  solution  p  =  0;  for  ptj  <p0  <  pcr,  there  are  exactly  three  solutions,  namely, 
the  trivial  solution  p  =  0  and  in  addition  two  bifurcated  solutions  each  involving  an 
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internal  cavity;  for  pCT  <  Po  <  Pt2>  there  are  exactly  four  solutions,  namely,  the  trivial 
solution  p  =  0  and  in  addition  three  bifurcated  solutions  each  involving  an  internal 
cavity;  finally,  for  po  >  pt2,  there  are  exactly  two  solutions,  namely  the  trivial  solution 
and  a  single  bifurcated  solution  involving  an  internal  cavity.  As  noted  earlier,  the  trivial 
solution  provides  a  local  minimum  for  X  and  hence  is  stable  for  0  <  po  <  per,  while  for 
Po  >  per  the  trivial  solution  provides  a  local  maximum  for  X  and  thus  is  unstable.  Again 
by  Property  A,  the  bifurcated  curve  (4.6.8)  provides  a  local  maximum  for  X  if 
Pt2  <  P  <  Pti  and  hence  the  associated  solution  is  unstable,  while  for  0  <  p  <  and 
p  >  Pti ,  this  curve  provides  a  local  minimum  for  X  and  hence  this  solution  is  stable. 
Thus  if  0  <  po  <  Pti  or  if  po  >  Pt2>  there  is  exactly  one  stable  equilibrium  solution,  while 
if  Pti  <  Po  <  Per  there  are  two  stable  equilibrium  solutions,  one  involving  an  internal 
cavity,  and  if  pcr  <  Po  <  Pt2  there  are  two  stable  equilibrium  solutions,  each  involving  an 
internal  cavity. 

To  address  the  transition  between  metastability  and  absolute  stability,  we  refer  to 
Figs.  4.6(b),(c),  where  pti  ,  Pa  >  Pti  ,  Pt2  are  as  before,  and  pci  denotes  the  value  of  p  > 
0  occurring  on  the  stable  branch  for  which  po(p)  =  pCr>  while  pC2  denotes  the  value  of 
p  >  0  such  that  po(p)  =  Per  occurs  on  the  unstable  branch.  In  Figure  4.6(c)  we  have 
graphed  the  energy,  S(p)  /  Pi  =  [  d>(p)  +  po(p)  'P(p)  ]  /  Pi ,  for  p  >  0  on  the  curve  shown 
in  Fig.  4.6(b).  Note  that  the  following  analysis  depends  only  on  the  sign  of 
X(p)  =  8  /  [4pi  (3  +  4ai )]  which  has  the  same  sign  as  8,  so  that  all  necessary  information 
on  X  can  be  extracted  from  Fig.  4.6(c).  We  shall  now  show  that  there  exists  a  unique 
value  of  po,  namely  p0  =  pi,  satisfying  ptl  <  pi  <  p^,  such  that,  if  pi  <  p,*,  the  trivial 
solution  is  absolutely  stable  for  po  <  pi ,  while  for  po  >  pi ,  the  absolutely  stable  solution 
occurs  on  the  bifurcated  branch.  (We  note  here  that  we  shall  confine  attention  to  the  case 
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Pi  <  Pc,,  in  what  follows;  stability  for  the  case  pi  >  pCT  is  analogous  to  the  situation 
previously  described  for  Figure  4.5(b)).  The  desired  result  is  established  by  the 
following  argument:  First,  since  2(0)  =  0,  we  have  2(pt2)  <  0  from  Property  B.  Second, 
from  Property  A,  we  know  that  for  p0  =  pCT  fixed,  p  =  pC2  corresponds  to  a  local 
maximum  of  the  energy  while  p  =  0  and  p  =  pcl  correspond  to  local  minima.  Thus,  in 
particular,  2(pc2)  >  2(0)  =  0.  [We  remark  here  that  the  energy  curve  for  fixed  p0  is 
distinctly  different  from  that  displayed  in  Fig.  4.6(c).  In  Fig.  4.6(c),  every  value  of  p 
corresponds  to  a  particular  maximum  or  minimum  of  the  energy  associated  with  a 
corresponding  value  of  the  load  p0.  To  see  this,  compare  Fig.  4.6(c)  to  Fig.  4.6(d),  which 
represents  the  energy  curve  2(p)  vs.  p  for  fixed  Po  =  Per-  Figure  4.6(d)  indeed  shows  that 
Pci  .  Pc2  correspond  to  extrema  of  the  energy  for  p0  fixed  at  p0  =  pcr-  Note  also  that 
although  Fig.  4.6(d)  implies  2(pcl)  <  2(0)  =  0,  this  is  not  necessarily  true  for  all  curves  of 
the  type  shown  in  Fig.  4.6(a).]  Returning  now  to  our  argument,  we  have  from 
£(Pt2)  <  0  >  2(pc2)  >  0,  and  Property  B,  that  there  exists  a  unique  value  of  p  =  p2, 
satisfying  Pt2  <  p2  <  Pc2>  for  which  2(p2)  =  0.  (See  Fig.  4.6(c)).  We  now  consider  the 
value  of  the  load  associated  with  p2,  defined  to  be  po(p2)  =  P2-  Again  from  Property  A 

//  ff 

for  Po  =  p2  fixed,  0  =  Z(p2)  >2(p2),  where  p2  is  the  largest  value  of  p  such  that 

Po(P)  =  P2-  (See  Figs.  4.6(b),(c)).  We  also  observe  that  the  preceding  along  with 

Property  B  give  2(ptl)>0.  Therefore,  by  Property  B,  there  exists  a  unique  value  of 

// 

p  =  Pi,  satisfying  ptl  <p!  <p2,  for  which  2(pi)  =  0.  The  associated  load  is  then 
Po(pi)  =  Pi-  (See  Figs.  4.6(b), (c)).  Thus,  the  energy  on  the  bifurcated  curve  is  such  that 
2(p)  <  0  for  0  <  p  <  p2  and  for  p  >  Pi,  while  2(p)  > 0  for  p2  <  p  <  px.  This  completes 
the  argument  since  the  energy  associated  with  the  trivial  solution  is  everywhere  zero. 
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Consider  now  the  previously  envisioned  quasi-static  loading  process.  We  thus  refer 
to  Figure  4.6(e),  where,  as  remarked  earlier,  we  have  pi  <  pCT.  If  one  supposes  that  the 
body  always  assumes  an  absolutely  stable  configuration,  then  the  trivial  homogeneous 
deformation  (4.3.2)  persists  for  0  <  po  <  Pi-  As  po  passes  through  pl5  the  foregoing 
notion  of  stability  predicts  the  sudden  appearance  of  a  cavity  of  finite  radius  c  =  Bpi.  As 
Po  is  further  increased,  the  cavity  radius  increases  in  a  continuous  fashion.  On  the  other 
hand,  one  may  suppose  that  the  trivial  homogeneous  deformation  (4.3.2)  persists  until  it 
becomes  unstable  at  po  =  Per-  In  this  case,  there  are  two  possibilities  as  po  passes 
through  pCT,  similar  to  those  discussed  for  Fig.  4.5(b).  First,  if  one  supposes  that  the 
body  then  assumes  an  absolutely  stable  configuration,  the  sudden  appearance  of  a  cavity 
of  finite  radius  c  =  Bpci  is  predicted,  with  cavity  radius  increasing  in  a  continuous 
fashion  as  po  is  further  increased.  On  the  other  hand,  one  may  suppose  that  continuous 
deformation  of  the  body  persists  as  long  as  possible  so  that  an  internal  void  is  nucleated 
at  Po  =  Per  and  the  cavity  radius  grows  continuously  from  zero  until  p0  =  pt2  and  then 
experiences  a  jump  from  c  =  Bpt2  to  c  =  Bpt2,  as  po  passes  through  p^.  The  cavity 
radius  then  increases  continuously  as  po  is  further  increased. 

We  have  seen  in  the  foregoing  that  there  are  four  distinct  types  of  bifurcation  that 
can  occur  (Figs.  4.3  -  4.6).  Note  that  these  four  possibilities  exist  regardless  of  which  of 
the  schematic  diagrams  (Fig.  4.1  or  Fig.  4.2)  applies.  In  each  of  the  diagrams,  points  in 
region  n  give  rise  to  only  one  type  of  bifurcation  (i.e.  that  shown  in  Fig.  4.4),  while 
points  in  region  I  have  three  types  of  bifurcation  associated  with  them  (i.e.  those  shown 
in  Figs.  4.3,  4.5,  4.6).  We  also  remark  that  although  there  corresponds  a  different 
diagram  resembling  either  Fig.  4.1  or  Fig.  4.2  for  every  value  of  sl\,  these  two  diagrams 
exhaust  the  possibilities.  For  <ao,  we  have  a  figure  similar  to  Fig.  4.1,  and  for 
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ai  >  ao,  we  have  a  figure  similar  to  Fig.  4.2.  We  recall  from  Chapter  3  that  for  a 
homogeneous  sphere  consisting  of  the  inner  material  only,  ao  is  the  transition  value  of 
the  anisotropy  parameter,  a,  which  governs  whether  bifurcation  for  this  sphere  is  to  the 
right  (a  <  ao)  or  to  the  left  (a  >  ao).  The  numerical  value  of  ao  is  given  by  (3.5.28). 
Finally,  we  observe  that  regardless  of  the  stability  criterion  adopted,  a  discontinuous 
growth  in  cavity  radius  is  predicted  in  three  of  the  four  types  of  bifurcation  discussed. 

We  conclude  this  section  by  addressing  more  closely  how  the  anisotropy  considered 
here  results  in  the  four  different  bifurcations  we  have  seen.  As  stated  above,  when  the 
parameter  pair  (co ,  f)  lies  in  region  II  in  either  of  the  schematic  diagrams  Figs.  4.1,  4.2, 
we  have  seen  only  one  type  of  bifurcation.  However,  when  (co ,  f)  lies  in  region  I,  there 
are  three  possibilities.  We  now  examine  the  role  of  the  anisotropy  parameters  ,  a2,  in 
determining  which  of  the  three  types  of  bifurcation  associated  with  regions  I  occurs.  For 
this  purpose,  we  refer  to  Figures  4.5(a),  4.7,  and  4.8.  These  three  figures  correspond  to 
the  same  value  o/a1;  the  same  parameter  pair  (co ,  f),  and  different  values  of  a2.  Note 
that  the  value  of  P  necessarily  changes  if  one  holds  ai  and  co  fixed,  and  varies  a2  (see  the 
comment  following  (4.6.15)).  As  the  value  of  a2  is  decreased  from  a2  =  3.7  in  Fig. 
4.5(a),  to  a2  =  1.475  in  Fig.4.7,  the  value  of  pti  approaches  that  of  pt2,  so  that  the  region 
in  which  there  is  more  than  one  stable  solution,  i.e.,  pti  <  po  ^  pt2»  becomes  much 
smaller.  When  a2  is  further  decreased  to  a2=0.14  in  Fig.  4.8,  the  curve  (4.6.8)  is 
monotone  increasing,  so  that  bifurcation  here  is  of  the  type  displayed  in  Fig.  4.3.  We 
also  note  that,  if  we  held  the  values  of  ai  ,  co,  and  f  associated  with  Fig.  4.6(a)  fixed  and 
varied  a2  as  above,  this  would  give  rise  to  figures  similar  to  those  shown  in  Figs.  4.7,  4.8. 
Thus,  we  conclude  that  it  is  the  magnitude  of  the  anisotropy  parameter  a2  associated 
with  the  surrounding  shell  that  is  important  in  determining  the  region  in  which  multiple 
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stable  solutions  occur  and  thus  the  monotonicity  of  the  curve  (4.6.8).  For  large  values  of 
a2  in  the  range  0  <  a2  £  5,  one  obtains  figures  like  Fig.  4.5(a)  or  Fig.  4.6(a),  while  for 
small  a2  one  obtains  figures  like  Fig.  4.8.  In  addition,  we  observe  from  Figs.  4.5(a), 
4.6(a)  that  the  magnitude  of  the  anisotropy  parameter  a\  associated  with  the  core 
material  is  also  crucial  in  determining  which  of  these  two  more  novel  types  of 
bifurcation  occurs.  These  two  figures  correspond  to  the  same  values  of  a2  ,  P,  and  f,  but 
different  values  of  a\  (and  again,  necessarily  different  values  of  ©).  The  above  statement 
follows  from  the  fact  that  the  point  at  which  the  curve  (4.6.8)  bifurcates  from  the  trivial 
solution,  i.e.,  at  p0  /  Hi  =  pCT  / Hi,  is  solely  dependent  upon  the  value  of  ai  (see  (4.6.11)2 
as  well  as  the  last  paragraph  of  Section  4.4).  Thus,  the  magnitude  of  ai  determines 
whether  pCT  <  ptl  or  pcr  >  pti,  which  in  turn  implies  whether  we  have  a  figure  like  Fig. 
4.5(a)  or  Fig.  4.6(a),  respectively.  Finally,  we  observe  that  for  a  composite  sphere 
corresponding  to  a  diagram  like  that  of  Fig.  4.1,  so  that  a^  <  a<>  =  0.257,  the  changes  in 
character  of  the  bifurcations  are  not  as  drastic  as  those  associated  with  Fig.  4.2  since 
variations  in  the  magnitude  of  ai  are  limited. 

4.7.  Stress  distribution 

We  conclude  this  chapter  with  a  discussion  of  the  stress  distribution  in  the 
composite  anisotropic  sphere.  Recall  that  at  the  beginning  of  Section  4.3,  we  saw  that 
the  stresses  corresponding  to  the  trivial  solution  (4.3.1)  were  given  as 

Ttr  =  Tee  =  =  po  .  (4.7.1) 

Thus,  prior  to  cavitation,  we  have  simply  the  constant  stress  distribution  (4.7.1). 
Subsequent  to  cavitation  ( c  >  0 ),  (4.3.14),  (4.3.8),  and  the  considerations  of  Appendix  A 
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give  the  radial  stress  as 


Tn(R)  =  J  (v3  -  l)"1  Wx  (v;  /- )  dv ,  on  0  <  R  <  B  ,  (4.7.2) 


where  we  recall  the  notation  (4.3.16),  (4.3.17).  On  using  (4.4.2),  we  may  write  (4.7.2) 
for  the  composite  anisotropic  sphere  as 


Trr(R)  =  { 


j  (v3  -  I)-1  W  J  (v)  dv  ,  on  0  <  R  <  A , 

(1  +  c3/R3)m 

(i  +  c3/b3)1/4  ^  (4.7.3) 

p0(  1  +  c3/B3  f/3  +  J  (v3  -  l)_1Wi  (v)  dv ,  on  A  <  R  <  B  . 

(1  +  c3/R3)1/S 


Now,  from  (4.2.16),  (4.2.17),  (4.3.3),  we  can  write  the  hoop  stresses  as 

Tee  (R)  =  T#  (R)  =  2 [(v2  -  v'4)  Wt  I R  -  (V2  -  v4)W2 1  r  -  V4W5  I R]  +  T„(R) ,  (4.7.4) 

where  Tn(R)  is  given  by  (4.7.2)  or  (4.7.3),  and  we  recall  the  notation  (4.3.9).  From 
(4.7.3),  (4.7.4)  and  our  assumptions  on  the  smoothness  properties  of  W,  we  see  that  TM  is 
a  continuous  function  of  R  on  0  <  R  <  B  ,  with  a  discontinuity  in  slope  at  the  interface 
R  =  A,  while  Tee  (R)  =  (R)  suffers  a  jump  discontinuity  at  the  interface.  The  stress 

distributions  (4.7.3),  (4.7.4)  are  given  explicitly  (see  (4.7.12),  (4.7.13)  below)  for  the 
material  model  (4.6.1)  and  are  illustrated  in  Figures  4.9  and  4.10.  Further  discussion  of 
these  figures  will  be  carried  out  below. 

Another  interesting  feature  of  the  stresses  immediately  after  cavitation  is  the 
presence  of  a  boundary  layer  near  the  cavity  wall  in  the  case  of  smooth  cavitation.  To 
see  this,  we  observe  that  for  fixed  R,  0  <  R  <  B,  letting  c  ->  0+  formally  in  (4.7.3)  results 
in  the  right-hand  side  of  (4.4.3),  so  that 
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lim  Tn-  (R)  =  pCT  ,  0  <  R  <  B  . 

c  — »  0+ 


(4.7.5) 


Since  T^O)  =  0,  (4.7.5)  shows  that,  when  smooth  cavitation  takes  place  (so  that  the 
cavity  radius  increases  continuously  from  zero  in  a  neighborhood  of  Per,)  the  radial 
stress  suffers  a  rapid  growth  near  the  cavity  wall  for  applied  dead  loads  po  slightly  larger 
than  pCT.  This  boundary  layer  behavior  is  shown  in  Figure  4.11  for  a  particular  case  of 
the  material  model  (4.6.1).  From  (4.7.4),  it  is  clear  that  a  similar  boundary  layer  exists  in 
the  stress  components  Tee  (R) ,  T^  (R) . 

We  next  consider  the  radial  stress  at  the  interface.  From  (4.7.3),  we  have 


Trr(A)=  J  (v3  -  l^wj  (v)  dv .  (4.7.6) 

(1  +  c3/A3)* 

We  wish  now  to  determine  a  sufficient  condition  on  the  strain-energy  so  that  the 
maximum  interfacial  stress  is  reached  at  the  point  at  which  cavitation  occurs.  This  will 
be  the  case  if  T„(A)  is  monotone  decreasing  with  respect  to  cavity  radius  c  >  0 .  It  is 
readily  verified  that 


dTJX(A)  i 


*  1 


(4.7.7) 


so  that  for  c  >  0 ,  dT^A)  /  dc  <  0  if  and  only  if 


W  i  (v)  >  0  for  all  v  >  1 


(4.7.8) 


Now,  from  (4.7.1),  we  have  prior  to  cavitation, 


Tn-(A)  =  po  ,  c  =  0 ,  (4.7.9) 

and  from  (4.7.5),  subsequent  to  cavitation, 

lim  Tn  (A)  =  pCT  =  lim  p0  (c)  ,  c  >  0 . 

c  Of  c  — ►  0+ 


(4.7.10) 
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There  are  now  three  possibilities  to  consider.  First,  in  the  case  of  smooth  cavitation,  the 
radial  stress  at  the  interface  T^A)  increases  linearly  with  p0  until  po  =  Per »  T^A)  is 
continuous  at  p0  =  pcr  *  and  then  decreases  for  Po  >  Per  •  (See  Figure  4.12).  Second,  for 
the  case  of  snap  cavitation,  where  a  cavity  of  finite  radius  c  =  Co  >  0  appears  at  a 
transition  load  p0  <  pcr  >  (Po  is  analogous  to  the  quantities  p  or  pj  discussed  in  Section 
4.6  in  connection  with  Figs.  4.4  or  4.6,  respectively)  equation  (4.7.10)  is  no  longer 
relevant.  The  interfacial  stress  T„(A)  increases  linearly  with  p0  until  Po  =  Po  »  at  which 
point  it  suffers  a  jump  discontinuity.  If  condition  (4.7.8)  holds,  it  is  easily  seen  that  the 
integrand  in  (4.7.6)  is  positive  on  the  range  of  integration,  so  that  the  jump  in  Tff(A)  is 
negative.  Thus  upon  appearance  of  a  cavity  with  finite  radius,  (c  =  Co  >  0 ,  po  =  Po  )>  the 
value  of  T„.(A)  necessarily  decreases.  (See  Figure  4.13).  Third,  consider  the  radial  stress 
at  the  interface  for  the  case  in  which  the  sphere  undergoes  a  bifurcation  of  the  type 
portrayed  in  Fig.  4.5.  Here  T^A)  increases  linearly  with  p0  until  Po  =  Per*  Tff(A)  is 
continuous  at  po  =  Per*  and  decreases  until  po  reaches  a  transition  value  (greater  than  pcr) 
where  it  again  suffers  a  jump  discontinuity.  (See  Figure  4.14). 

As  observed  in  [21],  the  foregoing  considerations  have  immediate  implications 
relating  to  the  issue  of  possible  debonding  at  the  material  interface  R  =  A.  (We  remark 
that  in  [21]  a  figure  of  the  type  displayed  in  Fig.  4.14  did  not  arise).  Suppose  that  the 
interface  bond  is  sustained  only  so  long  as  the  normal  stress  at  the  interface  remains  less 
than  a  threshold  value,  Td,  a  measure  of  the  strength  of  the  interface  bond.  Then 
debonding  would  occur  if 

Tjt(A)  =  Td  .  (4.7.11) 


Consider  again  a  quasi-static  loading  process  in  which  po  increases  slowly  from  zero.  If 
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Td  <  Per  for  smooth  cavitation  (Td  <  p0  for  snap  cavitation)  then  debonding  occurs  when 
Po  =  Td  and  cavitation  is  not  relevant.  However,  if  Td  >  pcr  for  smooth  cavitation 
(Td  >  po  for  snap  cavitation),  then  when  (4.7.8)  holds,  the  resulting  stress  relief  at  the 
interface  precludes  the  criterion  (4.7.11)  from  being  met  and  thus  eliminates  the 
possibility  of  interface  debonding. 

We  conclude  this  section  by  returning  to  the  material  model  (4.6.1).  Recalling  the 
notations  (4.6.4)  and  (4.6.7),  the  stresses  (4.7.3),  (4.7.4)  for  a  composite  sphere  composed 
of  such  a  material  can  be  written  as 


r 


Trr(R)  . 

1^1 


F[z(|-p)]  +  2a1(-^-G[z(|-p)])^,  0< R< A , 

(1  -  P)  F[  z(ap)  ]  +  p  F[  z(  -|p  )  ]  +  ax  (4.7.12) 


+  2(a2p  -  ax)  G[  z(ap)  ]  -  2a2P  G[  z(  |-p  )  ]|» ,  A <  R  <  B  , 


and 


Tee(R) 

lit 


T<^(R) 

lit 


z2(fp)-z  4(-^-p)  —  2ait z  8^P)_Z  4( R~P) ]  +  pi  )  ’  °-R<A» 
Pt  z2(f  P) "  z_4(f  P)  1 "  2a2P[  z"8(|-p)  -  z_4(f  P)  ]  (4.7.13) 


Trr(R) 

+ - 

Pi 


A<  R<  B . 


These  stresses  are  shown  in  Figures  4.9  and  4.10  for  parameter  values  corresponding  to 
the  bifurcation  diagrams  Figs.  4.4  and  4.6,  respectively.  In  Fig.  4.9,  these  stresses  are 
plotted  for  a  case  with  P  <  1  ( |12  <  Pi  )•  When  P  >  1  ( p2  >  px  ),  the  jump  in  Tee  = 
at  the  material  interface  R  =  A  is  in  the  opposite  direction  to  that  in  Fig.  4.9,  while  the 
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slope  discontinuity  in  has  opposite  character.  In  Fig.  4.10,  we  have  P  =  1.0,  and  the 
jump  discontinuity  in  Tee  =  T^  nearly  vanishes.  This  dependence  of  the  hoop  stress 
discontinuity  on  the  relative  sizes  of  |ii,  p2,  i.e.,  on  the  size  of  P  compared  to  unity,  was 
seen  in  the  hoop  stresses  corresponding  to  the  composite  neo-Hookean  sphere  discussed 
in  [21],  Thus,  the  anisotropy  considered  here  has  a  negligible  effect  on  the  jump 
discontinuity  in  the  hoop  stresses  at  the  material  interface.  This  can  be  seen  from 
(4.7.13)  since  it  is  the  Pz2  term  in  the  second  of  (4.7.13)  which  contributes  most 
significantly  to  the  discontinuity. 

Finally,  it  is  easily  verified  that  the  condition  (4.7.8)  holds  for  the  material  (4.6.1), 
so  that  we  do  have  normal  stress  relaxation  at  the  material  interface  subsequent  to 
cavitation  for  all  materials  of  the  type  (4.6.1).  Thus,  as  seen  earlier,  regardless  of  the 
stability  criterion  adopted,  there  are  essentially  three  possible  ways  in  which  the 
composite  sphere  may  deform.  First,  a  smooth  cavitation  may  take  place.  Second,  an 
immediate  snap  cavitation  at  some  load  less  than  the  critical  load  may  occur.  Third,  a 
smooth  cavitation  may  initially  take  place,  followed  by  a  "snap"  from  one  finite  cavity 
radius  to  another  at  some  load  greater  than  the  critical  load.  The  interfacial  normal 
stress  for  these  three  cases  are  shown  in  Figures  4.12,  4.13,  and  4.14  for  spheres  with 
parameter  values  corresponding  to  the  bifurcation  diagrams  Figs.  4.3,  4.4,  and  4.5, 
respectively. 
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Figure  4.1.  Semi-infinite  strip  0  <  /  = — j-  <1,  co  = 

2> 


d2w\l)/dv2  _  |i2(3+4a2) 
d2w\l)/dv2  Hi(3+4fli) 


>0. 


For  parameter  pairs  (co ,/)  in  region  I,  bifurcation  occurs  to  the  right,  while  in  region  II, 
bifurcation  occurs  to  the  left.  For  the  curve  shown, 
PcA*- 1  =  4(3+4ai) pcr ,  a.\  =  0.157  ,  pcr  =  0.18 .  Figure  4.1  corresponds  to  the  case  in 
which  bifurcation  for  the  homogeneous  sphere  composed  of  the  core  material  would  be 
to  the  right  (pcr  >  1/6).  The  point  shown  in  region  I  is  the  parameter  pair  associated  with 
Fig.  4.3. 
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Figure  4.2.  Semi-infinite  strip  0  <  /  = — j-  <1,  CD  = 


d2w\l)/dv2 

d2w\l)/dv2 


M-2(3+4a2) 

-  >  0. 

ji1(3+4a1) 


For  parameter  pairs  (to,/)  in  region  I,  bifurcation  occurs  to  the  right,  while  in  region  II, 
bifurcation  occurs  to  the  left.  Figure  4.2  corresponds  to  the  case  in  which  bifurcation  for 
the  homogeneous  sphere  composed  of  the  core  material  would  be  to  the  left  (pcr  <  1/6). 
The  point  shown  in  region  n  is  the  parameter  pair  (co ,/)  associated  with  Fig.  4.4,  while 
the  points  in  region  I,  ordered  with  increasing  co,  are  the  parameter  pairs  associated  with 
Figs.  4.6  and  4.5,  respectively. 
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Figure  4.3.  Variation  of  the  deformed  cavity  radius  p  =  c/B  with  applied  dead-load 
traction  p0/\X\  for  the  anisotropic  composite  sphere  described  by  (4.6.1).  Here 

U-2  4  ^  U.2 (3+4fl2) 

ai  =0.157,  a2  =0.2477,  B  =  —  =  1.0,  /=  =  cT3  =0.1,  ©=  ~~ - — -  =  1.1. 

H  lit  B3  lli(3+4ai) 

Alsop^/jj.!  =2.61. 


Figure  4.4.  Variation  of  the  deformed  cavity  radius  p  =  c/B  with  applied  dead-load 
traction  pQ/\i\  for  the  anisotropic  composite  sphere  described  by  (4.6.1). 

ll2  A3  n  U.2(3+4a  o) 

Here#!  =  1.475,  #2  =  3.7,  B=  —  =0.1,  /=  =  of3  =  0.1,  (0=  „  -  \  =0.2. 

1  2  H  111  J  B3  Hi  (3+4a ! ) 

Also  pcr/[Li  =  3.56  ,  (pt/[L  1  ,  pf)  =  (1.68  ,  1.04) ,  (p4k  ,  p)  =  (1.82  ,  1.45) . 


Figure  4.5(a).  Variation  of  the  deformed  cavity  radius  p  =  c/B  with  applied  dead-load 
traction  p0/\i\  for  the  anisotropic  composite  sphere  described  by  (4.6.1). 

l t2  4 3  ,  \i2(3+4a2) 

Here  ai=0.5,  a2  =  3.7,  p  =  —  =  1.0,  f=^T  =  a~3  =0.1,  (Q=  /<t  -  \  =3.56. 

1  2  H  Pi  J  B3  Hi(3+4fl!) 

Also pcr/\li  =  2.86  ,  (pn/Ui  »  Pn)  =  (3.49 , 1.21) ,  (pt2/\i\  ,  pf2)  =  (4.35 , 0.52) , 

Pri  =  0.25  ,  Pt2  =  2.01 . 


Figure  4.5(b).  Stability  zones  corresponding  to  bifurcation  diagrams  of  the  type 
displayed  in  Figure  4.5(a). 

u2  A3  -x  |i2(3+4a2) 

Here  a i  =  0.5,  a2  =  3.7,  P  =  —  =  1.0,  /=  ~  =  o'3  =  0.1,  co  =  — -  -  =3.56. 

H  111  B3  )liQ+4ai) 

Also  PtrA^l  =  3.71  ,  Ptr  =  0-29  ,  ptr  =  1.56  . 


(F/l/Pl  >  Prl  ) 


>  Pf2  ) 


Figure  4.6(a).  Variation  of  the  deformed  cavity  radius  p  =  c/B  with  applied  dead-load 
traction  p0/\i\  for  the  anisotropic  composite  sphere  described  by  (4.6.1). 


r-i.  A ^  q  P2(3+4fl2) 

Here  ai=2.5,  a2  =  3.7,  B  =  —  =  1.0,  /=  Ar  =  a-3  =0.1,  0)=  =  1.37. 

H  Pi  J  B3  Pi  (3+4*3 1 ) 

Also PcrA^i  =  4.30 ,  (pt i4n  ,  Pfi)  =  (3.50 ,  1.21) ,  (p,2> 111  ,  ptj)  =  (4.63 , 0.41) . 


Figure  4.6(b).  Variation  of  the  deformed  cavity  radius  p  =  c/B  with  applied  dead-load 
traction  p0fy\  for  the  anisotropic  composite  sphere  described  by  (4.6.1). 

Ll 7  A3  o  (J-2  (3+4(7  2) 

Here  a \  =  2.5,  a2  =  3.1,  B  =  —  =  1.0,  f=^r=  a-3  =0.1,  CO  =  ~~~  -  ■■  =  1.37. 

P  111  B3  lii(3+4ai) 


Also  (pi4ii  ,  Pi)  =  (3.76  ,  1.60) ,  (p2 4ii  >  P2)  =  (4.58 , 0.49) ,  p,2  =  2.18  ,  p2  =2.15  , 
pci  =  1.98  ,  p!  =  1.60 ,  pc2  =  0.64 . 


Figure  4.6(c).  Energy  P  =  — + - 'F(p),  on  the  curve  (4.6.8)  for  the 

Pi  Pi  Pi 
parameter  values  of  Figure  4.6(a). 


Figure  4.6(e).  Stability  zones  corresponding  to  bifurcation  diagrams  of  the  type 
displayed  in  Figure  4.6(a). 

Ll 2  A3  LI2  (3+4fl  2) 

Here fli  =2.5,  a2  =  3.7,  p  =  — =  1.0,  /=  ~  =  cT3  =0.1,  co=  =  1.37. 

lit  B3  lii(3+4ai) 

Also  (pi4ti  ,  Pi)  =  (3.76  ,  1.60),  (p24*i  »  P2)  =  (4.58 , 0.49),  p,2  =  2.18,  p2  =2.15, 
Pci  =  1.98  ,  px  =  1.60 ,  pc2  =  0.64 ,  pt2  /[Li  =  4.63  . 
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Figure  4.7.  Variation  of  the  deformed  cavity  radius  p  =  c/B  with  applied  dead-load 
traction  p0M\  for  the  anisotropic  composite  sphere  described  by  (4.6.1). 

Llo  A3  i  1^2  (3+4^2) 

Here  a,  =  0.5,  a2  =  1.475,  P  =  —  =  2.0,  /=  - 5-  =CT  =  0.1,  (0  =  =  3-56- 

Also pcr/\ii  =  2.86  ,  (pt\/\y\  ,  pit)  =  (4.81 , 0.94) ,  (pr2/U  1  »  Pti)  =  (4.86 , 0.68) . 


Figure  4.8.  Variation  of  the  deformed  cavity  radius  p  =  c/B  with  applied  dead-load 
traction  p0/\i\  for  the  anisotropic  composite  sphere  described  by  (4.6.1). 

LI  2  4  ^  U.2  (3+4(3  2) 

Here  a  j  =0.5,  a2  =0.14,  B  =  —  =  5.0, /=  ^V  =  oT3  =  0.1,  m=  „  -  N  =3.56. 

F  Hi  J  B3  lli(3+4ai) 

Also  Pcrfai  =  2.86 . 
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R/B 


Figure  4.9.  Variation  of  the  stresses  Trr(R),  Tqq(R)  =  T^(R)  with  undeformed  radius  R 
subsequent  to  cavitation  for  the  anisotropic  composite  sphere  described  by  (4.6.1).  The 
parameters  are  those  corresponding  to  Fig.  4.4,  i.e.,  a i  =  1.475  ,  a.2  =  3.7  ,  (3  =  0.1 , 
co  =  0.2  ,  /=  0.1 .  Also  po  /|ii  =  1.85,  pcr /|ii  =  3.56  .  At  the  material  interface,  the 
hoop  stresses  are  discontinuous,  and  the  radial  stress  is  discontinuous  in  slope. 


Figure  4.10.  Variation  of  the  stresses  Trr(R),  Tqq(R)  =  T^(R)  with  undeformed  radius 
R  subsequent  to  cavitation  for  the  anisotropic  composite  sphere  described  by  (4.6.1). 
The  parameters  are  those  corresponding  to  Fig.  4.6,  i.e.,  a\  =  2.5  ,  «2  =  3.7  ,  P  =  1.0 , 
(0=1.37,  /=0.1.  Also  Pq/\Jl\  =  3.95,  pcr/ Pi  =4.30 . 
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Trr(A)/[Li 


Figure  4.12.  Variation  of  the  radial  stress  at  the  interface  Trr(A)  with  applied  dead  load 
traction  p0  /\i\  for  the  anisotropic  composite  sphere  described  by  (4.6.1)  with  parameters 
corresponding  to  Fig.  4.3,  i.e.,  a\  =  0.157  ,  a2  =  0.2477  ,  P  =  =  1-0 ,  €9=1.1, 

/  =  0.1.  Here/?cr/|ii  =2.61 . 
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TAA)/  ui 


Figure  4.13.  Variation  of  the  radial  stress  at  the  interface  Trr(A)  with  applied  dead  load 
traction  pQ  /\i\  for  the  anisotropic  composite  sphere  described  by  (4.6.1)  with  parameters 
corresponding  to  Fig.4.4,  i.e.,  a\  =  1.475,  a2  =  3.7,  P  =  M-2Aii  =0.1,  CD  =  0.2, 

T)  Pt  Per 

f  =  0.1 .  A  jump  discontinuity  exists  at  =  1.82.  Here  —  =  1.68  ,  =  3.560 . 
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Trr(A)/\Lx 


Figure  4.14.  Variation  of  the  radial  stress  at  the  interface  T„.(A)  with  applied  dead  load 
traction  p0  /(ij  for  the  anisotropic  composite  sphere  described  by  (4.6.1)  with  parameters 
corresponding  to  Fig.  4.5,  i.e.,  a  \  =  0.5  ,  a  2  =  3.7  ,  |5  =  1.0 ,  (0  =  3.56  ,  /  =  0.1 .  A 

Ptr  Pt  1  Pt  9  Per 

jump  discontinuity  exists  at  =  3.71.  Here - =  3.49  ,  - =  4.35  ,  - =  2.86 . 

Hi  Hi  Hi  Hi 
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Chapter  5 

CONCLUSIONS 


We  now  wish  to  summarize  some  of  the  important  and  interesting  results  of  this 
dissertation.  The  type  of  anisotropy  considered  here  not  only  gives  rise  to  novel  types  of 
behavior  for  the  cavitation  problem,  but  has  also  proved  a  worthwhile  topic  for 
investigation  in  its  own  right.  As  mentioned  earlier,  a  comprehensive  study  of 
constitutive  inequalities  is  carried  out  in  [41]  for  compressible  isotropic  materials  as  well 
as  for  compressible  materials  of  arbitrary  symmetry  (see  pp.  153-171  of  [41]).  As  no 
comparable  fundamental  studies  of  inequalities  have  been  carried  out  for  incompressible 
materials,  the  results  of  Section  2.5  are  especially  useful  in  several  ways.  First,  we  have 
determined  the  relationships  between  the  infinitesimal  elastic  moduli  and  the  nonlinear 
stored-energy  function,  and  thereby  determined  restrictions  on  the  nonlinear 
incompressible  transversely  isotropic  strain  energy  density,  W(tj,  I2, I3, I4, 15 ), 
consistent  with  the  classical  inequalities  of  the  associated  linear  theory.  In  addition,  the 
methods  of  Section  2.5  employed  to  obtain  these  results  seem  to  be  readily  extendable  to 
materials  with  various  symmetries.  That  is,  the  techniques  of  this  section  ought  to  be 
applicable  to  any  material  with  assumed  symmetries  for  which  there  correspond  linearly 
elastic  moduli  with  known  restrictions,  and  thus  lead  to  useful  physically  motivated 
inequalities  on  the  nonlinear  strain  energy  density.  Finally  in  Section  2.5.3,  we  propose  a 
nonlinearly  elastic  transversely  isotropic  material  model  and  subsequently  give  all  of  the 
restrictions  for  physically  reasonable  response  which  follow  from  the  known  constitutive 
inequalities  of  the  linear  theory.  As  analytic  models  based  on  experimental  results  are 
limited  even  in  isotropic  finite  elasticity,  the  need  for  counterparts  based  instead  on 
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mathematical  considerations  is  unavoidable,  especially  for  illustrative  purposes.  In  view 
of  this,  the  class  of  materials  proposed  in  Section  2.5.3  at  least  goes  beyond  a  purely 
mathematical  construction  to  include  conditions  for  physically  reasonable  response,  and 
thus  should  be  a  useful  material  model  in  other  problems  of  incompressible  transversely 
isotropic  finite  hyperelasticity. 

Chapters  3  and  4  concentrate  on  the  cavitation  problem  for  incompressible 
nonlinearly  elastic  transversely  isotropic  spheres.  As  stated  in  Chapter  1,  this  work  was 
motivated  by  two  previous  investigations  ([4],  [20])  which  found  the  dramatic  material 
instability  of  snap  cavitation  to  occur  in  compressible  homogeneous  transversely 
isotropic  spheres  and  incompressible  composite  isotropic  spheres  respectively.  In 
Chapter  3  we  were  successful  in  obtaining  more  explicit  results  on  the  effects  of 
anisotropy  on  cavitation,  and  in  particular,  on  the  snap  cavitation  phenomena  which  we 
found  (as  expected)  to  occur  for  an  incompressible  homogeneous  transversely  isotropic 
sphere.  It  would  be  interesting  to  investigate  further  the  connection  between  the  problem 
of  Chapter  3  and  that  of  [20,  21]  as  suggested  by  the  concept  of  homogenization  in 
composite  material  mechanics. 


In  Chapter  4,  we  obtained  many  interesting  results  on  cavitational  instabilities  for  an 
incompressible  composite  transversely  isotropic  sphere.  We  summarize  some  of  these  in 
what  follows: 


( i )  If  the  core  material  alone  would  undergo  a  bifurcation  to  the 


right  I 
left 


,  then  the 


composite  sphere  will  also  bifurcate  to  the 


right  I 
left 


provided  the  surrounding  shell  is 


stiffer 

more  compliant 


in  tension  in  the  radial  direction  than  the  core  material. 
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(ii)  If  the  core  material  alone  would  undergo  a  bifurcation  to  the 


"  right] 

•<  left  k  then  the 

l  left  J 


composite  sphere  will  bifurcate  to  the 


'  left  ^ 

■s  right  ► 

L  left 


provided  there  is  a  sufficiently 


r  small! 

-  small  r 
>  large  J 


core  surrounded  by  a  shell  which  is 


more  compliant 
stiffer 
stiffer 


in  tension  in  the  radial  direction 


than  the  core  material. 

(Hi)  If  the  core  material  alone  would  undergo  a  bifurcation  to  the  right,  the  composite 
sphere  will  also  bifurcate  to  the  right  provided  the  volume  fraction  of  the  core  to  the  total 
volume  exceeds  1  /  6pcr . 


In  Chapter  4,  we  also  derived  an  alternative  formula  for  the  critical  load  at  which 
bifurcation  from  the  trivial  solution  occurs  (equation  (4.5.21))  involving  the  stored- 
energy  function  itself  rather  than  its  derivative.  This  formula  aided  in  obtaining  the 
energy  and  stability  results  for  cavitation  at  the  end  of  Section  4.5  for  general 
transversely  isotropic  composite  spheres,  and  has  the  additional  benefit  of  providing  a 
way  to  calculate  the  critical  load  if,  for  example,  the  derivative  of  the  strain  energy 
involved  in  the  original  formula  is  not  readily  obtainable.  While  equation  (4.5.21)  results 
from  a  simple  application  of  the  integration  by  parts  formula  and  the  growth  assumption 
(4.5.20),  the  author  has  been  unable  to  find  such  a  relationship  documented  in  any  of  the 
previous  work  on  cavitation  given  in  the  references.  We  note  that  (4.5.21)  is  clearly 
valid  for  all  combinations  of  homogeneous  or  composite,  isotropic  or  transversely 
isotropic  incompressible  materials  with  the  appropriate  interpretation  for  the  strain 
energy  function  appearing  in  it 
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The  illustrative  example  of  Section  4.6  provided  several  new  descriptions  of  how 
cavitation  in  a  composite  sphere  may  progress  in  a  quasi-static  loading  process.  In 
particular,  we  found  three  possible  way  for  the  sphere  to  deform  when  a  bifurcation  to 
the  right  occurs,  two  of  which  involve  material  instabilities  associated  with  a  jumping 
phenomenon  similar  to  the  snap  cavitation  phenomenon  originally  found  in  [4]  and  [20] 
(recall  Figures  4.5  and  4.6  and  the  discussions  regarding  them).  We  were  also  successful 
in  determining  more  explicit  information  on  how  the  anisotropy  considered  in  this 
dissertation  results  in  the  four  possible  types  of  bifurcation  found  here  (recall  the 
discussion  at  the  end  of  Section  4.6).  Finally  in  Section  4.7,  we  were  able  to  obtain  a 
(general)  condition  on  the  strain  energy  function  guaranteeing  stress  relaxation  at  the 
material  interface  after  cavitation,  and  thus  relate  this  to  the  possibility  for  cavitation  to 
prevent  debonding  of  the  composite  at  the  interface. 

We  wish  now  to  comment  on  the  effects  of  surface  energy  on  cavitation,  since  void 
formation  is  obviously  accompanied  by  creation  of  additional  surface.  Fracture  initiation 
due  to  a  hydrostatic  tensile  stress  field  was  analyzed  in  [45]  using  an  energy 
minimization  approach  involving  three  energy  quantities:  the  potential  energy  of  the 
external  tensile  stress,  the  strain  energy,  and  the  surface  energy.  The  problem  addressed 
in  [45]  was  to  consider  the  large  strain  and  instability  analysis  for  a  small  spherical  cavity 
of  initial  radius  ^  located  in  an  isotropic  incompressible  elastic  medium  (a  neo-Hookean 
sphere)  and  subjected  to  a  hydrostatic  tensile  stress,  p,  at  some  large  initial  outer  radius 
bo,  such  that  b0  » ac.  (Gent  and  Lindley  [14]  first  studied  the  above  problem  by 
calculating  the  stress  at  infinity  needed  to  expand  a  finite  cavity  in  an  infinite  expanse  of 
material  indefinitely.  As  remarked  by  Ball  [5],  this  problem  is  equivalent  to  the 
cavitation  problem  because  "an  infinitesimal  hole  in  a  finite  piece  of  material  behaves 


141 


like  a  finite  hole  in  an  infinite  expanse  of  material."  The  underlying  reason  for  this  is  that 
the  equilibrium  equation  and  boundary  condition  are  invariant  under  the  transformation 

(X,x)  — >  (cX.cx)1).  It  is  then  shown  in  [45]  (see  eqn.  (31))  that 


where  pc  is  the  critical  load  needed  to  expand  the  cavity  indefinitely  ,  E  is  Young’s 
modulus,  and  ae  is  the  hoop  stress  corresponding  to  the  neo-Hookean  sphere.  As  is  well 
known  in  cavitation  problems  however,  the  hoop  stress  at  the  cavity  wall  becomes 
infinite  (see  e.g.  [21],  or  see  Sections  3.7,  4.7  and  Figures  3.6,  4.10  of  this  dissertation). 
(Note  that  infinite  stresses  do  not  present  a  difficulty  within  the  framework  of  nonlinear 
elasticity.)  We  therefore  recover  the  formula  for  the  critical  load  for  the  neo-Hookean 
material  from  (5.1)  on  letting  (Te  -»  «>,  i.e., 

Pc  =  |E.  (5.2) 

Thus  the  critical  load  (5.2)  is  exactly  that  found  from  an  analysis  which  does  not 
incorporate  surface  energy  effects.  That  is,  since  E  =  3(1  for  incompressible  materials, 

the  critical  load  (5.2)  for  the  neo-Hookean  material  can  also  be  represented  as  pc  =  —  (i, 

which  agrees  with  equation  (3.5.21)  when  a  =  0.  As  reported  by  Gent  in  [13],  it  is  also 
interesting  to  note  that  experiments  on  rubbery  solids  seem  to  indicate  good  correlation 
between  the  predicted  (5.2)  and  the  observed  critical  load  for  the  types  of  materials 
tested. 


'  [5,  p.607] 


142 


In  conclusion,  we  indicate  some  additional  areas  of  future  research  related  to  the 
considerations  of  this  thesis.  Most  studies  of  cavitation  within  elasticity  theory, 
including  this  one,  have  examined  this  problem  in  the  context  of  spherically  symmetric 
solutions.  Recently  however,  several  authors  have  obtained  interesting  results 
concerning  cavitation  in  isotropic  elasticity  when  the  type  of  material  instability  is  not 
restricted  to  radially  symmetric  void  formation  (see  e.g.  [3],  [26]).  It  would  therefore  be 
of  great  interest  to  further  pursue  the  goals  of  Chapters  3  and  4  of  this  dissertation  to 
determine  the  effects  of  anisotropy  on  cavitation  by  relaxing  the  assumption  of  radially 
symmetric  solutions.  Another  worthwhile  generalization  would  involve  the  further  study 
of  cavitation  in  compressible  anisotropic  nonlinear  elasticity  following  the  early  work  of 
[4].  Finally,  the  considerations  of  Chapter  2  and  the  consequences  found  for  the 
cavitation  problem  arising  from  anisotropy,  suggest  that  investigations  of  the  effects  of 
anisotropy  for  other  problems  of  nonlinear  elasticity  may  yield  many  novel  phenomena 
and  exciting  new  results  in  the  field  of  finite  hyperelasticity. 
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Appendix  A:  Verification  of  (3.3.15) 


To  establish  (3.3.15),  we  rewrite  (3.3.12)  as 
a  3%  B 


Po  = 


1+ 


BJ 


f  -4  [  (v_ 7  (s)  -  v-1  (s) )  Wi  I  +  (v-5  (s)  -  v(s))  W2  I 
o  L 


+  v-7(s)W5  ij 


ds 

s 


(Al) 


where  we  recall  the  notation  (3.3.8).  From  (3.3.3),  we  have 


so  that 


Then 


and  on  using  (A2)  we  obtain 


v3(s)-l  =  -^- 


2  dv  _  -3c 


3v-t-  = 


ds  s4 


2  .  c3  ds 

V  dv  = - r - , 

s3  S 


(A2) 


ds  v2 


s  1  -  v3 


dv . 


(A3) 


Also,  recalling  the  notation  (3.3.13),  (3.3.14),  the  chain  rule  gives 


W!  (v)  =  (-4v-5  +  4v)  W!  I  +  (4V3  -  4v-3 )  W2 1  -  4v"5  W5  I  .  (A4) 

R  R  R 


Thus,  substituting  (A3)  into  (Al)  and  comparing  with  (A4),  we  see  that  the  integrand  in 
(Al)  is  Wi  (v)/(l  -v3)  which,  on  changing  the  limits  of  integration  according  to  (A2), 
establishes  the  result  (3.3.15)  as  desired. 
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Appendix  B:  Verification  of  (3.5.30) 


Let 


p3  =  e  1  ,  0  <  e  <  1 , 

so  that  e  ->  0  implies  p3  ->  and  rewrite  (3.5.23)  as 


P  =  2(l  +  e-1)2/N 


(1  +  e-1)-1"*  +  -7  (1  +  £"lr4/3 
4 


(Bl) 


+  2a  i 


-  (1  +  e'1)-1^  +  4-  (1  +  e_1r%  -  T  (!  +  e_1)~4/3  +  T  (1  +  e_1)~5/3  +  T  (1  + 

6V3  2  4  5  8 


1 

V3 


2(1  +e-1)h  + 1 

V3  (1+e-1)* 

arctan 

V3 

+  arctan 

2  +  (l  +  £-1)^ 

lb 

► 

(B2) 


Substituting 


(1  +  e-1)m  =  (l±i )“  =  £_m(l  +  e)m  =  £— m(l  +  m£  +  0(E2) ) 
e 


(B3) 


into  (B2)  where  appropriate  and  taking  the  limit  as  e  — >  0  gives  to  leading  order, 

1  An  .  1  s n  .  1 


lim  P  =  lim  2  E 
e~>o  e^o 


e>/s  +  1  e4/3  +  2a 

4 


_  evs  +  i  £«  _  JL  £«  +  ±  e»  +  -V* 

6^3  2  4  5  8 


1 

V3 


+  1  £ 
arctan - — —  +  arctan 


16 


€ 


2  +  £ 


,-16 


(B4) 


=  lim 
£  — >  0 


+  2a 


lA 


-  E~~  -  £-"  +  %  -  -p- 

6V3  V3 


arctan - 


2e“’/3  +  1  .  .  V3  r~* 


+  arctan - 


2  +  £ 
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=  lim  i  2e  'A  +  4a 
e  o 


Now  let 


5k  _  t* _ 1_ 

6^  V3 


2e_w  +  1  ,  «.  43  e 

arctan - —  +  arctan 


■V3 


4 3 


2  +  e 


-'/3 


y+2a.(B5) 


8  =  eV4  so  that  8  = 


(B6) 


Then  8  -+  0  implies  p  and  (B5)  becomes 


lim  P  =  lim  28  1 

5->0  6-»0 


+  lim  4a 

8->o 


57t  8-  1 


6 43  43 


28_1  +  1  .  ^  43  b~l 

arctan - 7= - h  arctan 


43 


2  +  8 


1 


82 


+  2a.  (B7) 


Using  L’Hopital’s  rule,  the  second  term  on  the  right  in  (B7),  after  some  manipulation, 
becomes 


4a  lim  -7 
6-»o2 


5-2 


8~2  +  8"1  +  1 


-1 


8-1 , 


or 


2a  lim 
8  ->  0 


8"2  +  8"1 


8“2  +  8"1  +  1 


(B8) 


Now  using  (B6)2,  the  above  expression  becomes 


-  2a  lim 

p  — >  <* 


P2  +  P 


p2  +  p  +  1 


=  -  2a . 


(B9) 


Therefore,  on  using  (B6)2  and  substituting  (B9)  into  the  second  term  on  the  right  in  (B7), 


we  obtain 
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lim  P(p)  =  lim  2p  -  2a  +  2a  =  lim  2p  ,  (BIO) 

p  — )  oo  p  — >  eo  p  — ^  oo 

and  thus  the  desired  result  (3.5.30),  i.e., 


P  -  2p  as  p  — >  oo . 


(Bll) 
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Appendix  C:  Verification  of  (3.5.16),  (3.5.19),  and  (3.6.6)-(3.6.8) 

Here  we  present  the  details  of  the  derivation  of  equations  (3.5.16),  (3.5.19),  and 
(3.6.6)  -  (3.6.8).  We  consider  the  indefinite  integrals  which  are  needed  to  evaluate  the 
definite  integrals  (3.5.14),  (3.5.18),  and  (3.6.4).  (Constants  of  integration  will  be 
ignored.) 

(i)  Evaluation  of  (3.5.14): 

As  in  [6],  we  write 
r  v  -  v-5  j  r  v6  -  1  j 

3  v3- 1  3  vV-l) 

=  J  —  —  dv  =  J  (v-2  +  v-5)  dv 

=  -  v-1  -  lA  v-4  .  *  (Cl) 

On  using  (Cl),  the  definite  integral  in  (3.5.14)  is  immediately  evaluated  to  yield  the 
desired  expression  (3.5.16). 

(ii)  Evaluation  of  (3.5.18): 

First,  we  record  here  the  indefinite  integrals  (2.128,  1.)  from  the  integral  tables  of 

Gradshteyn  and  Ryzhik,  1990: 

j.  dv  _  _ _ 1 _ b(3m  +  n  -  4)  r  dv 

3  ~~  (n  - 1)  a  vn— 1  zf 1  a(n  “  1)  v^zf  ’  (C2) 

where  Z3  =  a  +  bv3,  a  #0  ,  n^l,  b,  m  >  0  are  constants.  Now  the  integral  of  concern  in 
(3.5.18)  can  be  written  as 
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1=  f  v  5  —  9  dv  =  J-=-f - j  9 

J  V3_J  J  y5  (v3  _  J  v9 


dv 


dv 


(v3-!) 


=  Ii-I2. 


(C3) 


To  evaluate  I2,  we  use  (C2)  with  n  =  9,  m  =  1,  b  =  1,  a  =  -1  to  get 


I2  = 


1 


dv 


8v8  J  v6  (v3  -  1) 


(C4) 


The  integral  in  (C4)  is  evaluated  as  in  [5]  by  using  (C2)  with  n  =  6,  m  =  1,  b  =  1,  a  =  -1 
to  get 


dv 


1 


J 


dv 


J  v6  (v3  -  1)  5V5  J  v^v3  —  1) 

Reapplying  (C2)  with  n  =  3,  m  =  1,  b  =  1,  a  =  -1  to  the  integral  in  (C5),  we  obtain 


(C5) 


I 


dv 


1  +i dv 


v3  (v3  -  1)  2v2  ^  v3  —  1 


(C6) 


and  this  last  integral  can  be  evaluated  using  standard  integral  tables.  For  example,  from 
(2.126, 1.)  of  Gradshteyn  and  Ryzhik  we  obtain 


r  dv  1 ,  (1  +  v  +  v2)1/4  1  4  V3"  v 

J  — - =  -  —In  - ; — - - arctan 


v3-l 


v-  1 


V3 


2  +  v 


(C7) 


Then  on  using  (C4),  (C5),  (C6),  (C7),  we  have 

1111,  (1  +  v  +  v2)  1  4  ^v 

1 2  =  — r  +  — r  +  — T  ~  Tln  - o  — F7  arctan  - 

8v8  5v5  2v2  6  (v  - 1)2  VJ  2  +  v 


(C8) 


To  evaluate  Ii ,  we  use  (C2)  with  n  =  5,  m  =  1,  b  =  1,  a  =  -1  to  get 


li  =  +  J  -5-? - , 

4v4  J  v^v3-!) 


(C9) 


and  reapply  (C2)  with  n  =  2,  m  =  1,  b  =  1,  a  =  -1  to  the  integral  in  (C9)  to  obtain 
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dv 


=  1+J- 

V  J  ' 


vdv 


J  v^v3  -  1)  v  ■*  v3  -  1 

This  last  integral  is  given  by  (2.126, 2.)  of  Gradshteyn  and  Ryzhik  to  be 

r  v  dv  1  ,  (v  - 1)2  ,  1  „  2v  + 1 

J  — - =  —  In  — \ arctan 


v3-!  6  v2  +  v  + 1  ^3 


€ 


(CIO) 


(Cll) 


Then  on  using  (C9),  (CIO),  (Cll),  we  have 


Il=^  +  I  +  ilnJ-Di 


4v* 


6  \  +  v  +  1 


1  2v+ 1 

—  arctan  — =— 

S3  S3 


(Cl  2) 


Substituting  (Cl 2),  (C8)  into  (C3)  and  noting  that  In 


obtain 


T  1  ,  1  ^  1  ,  2v+ 1  1 

I  =  — — -  h - h  — 7=-  arctan 


(v-D2 

V2  +  V  +  1 


,  (v2+v+l) 

-  In  - - we 

(v-1)2 


4v4  v  ^3 


V3  8v8  5V5  2v2  € 


1.1  „  V3"v 

+  ~  arctan— - .  (C13) 


2  +  v 


On  using  (C13)  and  l’Hopital’s  rule  where  appropriate,  the  definite  integral  in  (3.5.18)  is 
immediately  evaluated  to  yield  the  desired  expression  (3.5.19). 

(Hi)  Evaluation  of  (3.6.4): 


We  consider  the  integral  in  (3.6.4)  as 

T  f  v-2  +  2v4  -  3v2  r  v-6  -  2v~2  +  v2 

W  •dV  +  aJ 


(V3  -  l)3 


(V3  -  l)3 


dv  =  I4  +  a  I5 


(C14) 


Then 


T  r  v  2  +  2v4  -  3v2  ,  2V3  -v-1 

I4  =  J - ^ ; - dv  = 


(v3  -  l)2 


(l-v^v 


(Cl  5) 


as  given  by  equation  (3.3)  of  [20].  To  evaluate  I5,  we  consider 
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I5=J 


dv 


-2J 


dv 


vw  (v3  -  l)2  J  v2(v3  -  l)2 


v2  dv 
(V3  -  l)2 


=  I6-2I7+I8.  (C16) 


Then  I8  can  be  directly  integrated  to  yield 

Ig  = 


-1 


3(v3  - 1)  ’ 

and  I7  is  given  by  (2.131, 2.)  of  Gradshteyn  and  Ryzhik  as 


(C17) 


It  =J 


dv 


v^v3  - 1)2 


1  4  2 

--•7  v 

v  3 


1 


1 _ t  v  dv 

-  1  3  v3  -  1 


v3 


(Cl  8) 


where  this  last  integral  in  (C18)  was  evaluated  earlier  in  (Cll).  Finally,  (2.128,  2.)  of 
Gradshteyn  and  Ryzhik  gives  1$  as 


I6=| 


dv 


-1 


8 


J 


dv 


(C19) 


v6  (v3  -  l)2  Sv^v3  -  1)  3  J  v6(v3  -  1) 

where  this  last  integral  in  (C19)  was  evaluated  earlier  by  (C5)  -  (C7).  Combining  (C14) 
(C19),  (Cll),  and  (C5)  -  (C7),  we  obtain 


T  2v3  -  v  -  1  ,  , 

I3  = - ^ - H  a  i 

( 1  —  v3  )v 


-1 


8 


3V5  (v3  -  1)  3 


1  +  1 


1  ,  (l+v  +  v2)14  1 

—  In  - — - =•  arc  tan  _ 

3  v-1  yg  2  +  v 


5V5  2v 


_ 1 

1 

K> 

1  4  2 

j  v3  ”  1 

v  3  V 

(C20) 


8 


1  1n  (v-1)2  ,  1 

6  v^  +  v  +  1  V3 


arctan 


2v+  1 

V3 


3(v3  - 1) 


On  observing  that  the  logarithm  terms  in  (C20)  cancel  identically,  and  on  using  (C20) 
and  l’Hopital’s  rule  where  appropriate,  the  definite  integral  in  (3.6.4)  is  immediately 
evaluated.  The  desired  expressions  (3.6.6)  -  (3.6.8)  can  then  be  obtained  by  utilizing  the 
definitions  (3.5.22)  and  (3.6.5). 
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Appendix  D:  Verification  of  (4.5.9) 


Here  we  present  the  details  of  the  derivation  of  (4.5.9).  Recall  from  (4.5.6)  and 
(4.4.11)4  that 

gp)=  , -t-(p)  ^p^p)  .  pan,  (CD 

d2W  (1)  /  dv2 

where  <b(p) ,  'F(p)  are  given  by  (4.5.7),  (4.5.8),  respectively.  From  (4.5.8),  it  is  easily 
seen  on  taking  derivatives  that 

'F(Of)  =  #(Of)  =  0  ,  ¥  (Of)  =  -6 .  (D2) 

•  M  ••• 

We  thus  confine  attention  to  the  evaluation  of  d> ,  <t>  ,  d>  as  p  — >  0+  to  establish  the 
desired  results. 


(i)  Verification  ofL( 0+)  =  0: 

By  the  chain-rule,  we  have  d>(p)  =  3p2dd>/dp3,  where  dd>/dp3  is  given  by  (4.5.17), 
so  that 


<J>(Qf)  =  lim  3p2dd>/dp3  =  0 , 

p->Of 


(D3) 


where  we  have  used  (4.5.22)  in  the  final  step  in  (D3).  From  (Dl)  we  have 


£<p) — ,f(p)  +?*»), 

d2W  (l)/dv2 


(D4) 


so  that  upon  evaluation  of  (D4)  as  p  — »  0+  and  substitution  from  (D2)j,  (D3),  we  have 
the  desired  result  (4.5.9)  i.e., 


E(0+)  =  0 . 


(D5) 
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( ii)  Verification  oft  (0+)  =  0: 

By  the  chain-rule,  we  have 

••  dO  a  d^O 

t(p)=6pd7+9P  W'  (D6) 

where  again  the  first  derivative  on  the  right  is  given  by  (4.5.17),  and  the  second  by 
(4.5.23).  Substituting  (4.5.23)  into  (D6)  and  taking  the  limit  as  p  -»  0+  gives 


••  d<T> 

0(0+)  =  lim  6p — t- 

P  — >  o+  K  dp3 


+  lim  3p 

p  — >  o+ 


W?(  (1  +  a3p3)^  )  -  w|(  (1  +  a3p3)*  )  W?(  (1  +  p3)*  ) 


3^3\2/3 


(1  +  crp3) 


a+p3)2/ 


where  we  note  from  an  argument  analogous  to  that  used  in  (D3)  that  the  first  term  on  the 
right  in  (D7)  is  zero,  while  (3.4.2)!,  i.e.,  Wt(l)  =  0,  (i  =  1,2),  gives  that  the  second  term 
in  (D7)  is  also  zero.  Thus, 


0(0+)  =  0 


From  (D4)  we  have 


S(p)  =  2-f(P)  ,+pfop),  (D9) 

d2W  (1)  /  dv2 

so  that  upon  evaluation  of  (D9)  as  p  ->  Of  and  substitution  from  (D2)2,  (D8),  we  have 
the  desired  result  (4.5.9)2,  i.e., 


2(0+)  =  0 


(DIO) 


(in)  Verification  oft  (0+)  =  6(pcr  -  p): 


By  the  chain-rule,  it  can  be  shown  that 
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•••  dO  'x 

4>(p)  =  6^-  +  54p3-?-%- 

dp3  d(p3r 


f  27p6 


d(p3)3 


(Dll) 


From  (4.5.22),  we  know  that  the  limit  as  p  ->  0+  of  the  first  term  on  the  right  in  (D1 1)  is 
6pcr,  so  that  the  desired  result  will  follow  from  showing  that  the  last  two  terms  in  (Dll) 
are  zero  as  p  -»  Of.  We  first  address  the  second  term  on  the  right  in  (Dll).  It  is  easily 
seen  that,  on  substitution  from  (4.5.23)  and  on  using  (3.4.2)!,  we  have 


lira  54p3 

p  ->Of 


d2<l> 

d(p3)2 


=  0. 


(D12) 


To  address  the  last  term  in  (Dll),  we  first  differentiate  (4.5.23)  to  obtain 


where,  in  (D13)  we  have  introduced  the  notation 

[f-g]  (x)  =  f(x)  -  g(x) .  (D14) 

Since  the  quantities  d2W*(l)  /  dv2,  (i  =  1,2),  are  finite  from  (3.4.2)2  and  W*i(l)  =  0  from 
(3.4.2)i,  substitution  from  (D13)  into  the  last  term  in  (Dll)  and  taking  the  limit  as 
p  — >  Of  results  in 
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lim  27p6 

p  — >  Of 


d3$ 

d(p3)3 


=  0. 


(D15) 


Thus,  we  have  shown  that  the  last  two  terms  on  the  right  in  (D1 1)  vanish  as  p  — » 0+  and 
so 


3>  (Of)  =  6Pcr .  (D1Q 

From  (D9)  we  have 


2(P)  = 


$(p) 


*  i 


d  W  (1)  /  dv 


+  P*<P). 


(D17) 


so  that  upon  evaluation  of  (D17)  as  p  — >  0+  and  substitution  from  (D2)3,  (D16),  we  have 
(recalling  (4.4.11)3)  the  desired  result  (4.5.9)3,  i.e.. 


*f  (Of)  =  6(Pcr-p). 


(D18) 


155 


REFERENCES 


[1]  Abeyaratne,  R.  and  Hou,  H.-s.  (1989).  Growth  of  an  infinitesimal  cavity  in  a  rate- 
dependent  solid.  J.  Appl.  Mech.,  56 , 40-46. 

[2]  Abeyaratne,  R.  and  Hou,  H.-s.  (1991a).  Void  collapse  in  an  elastic  solid.  J. 
Elasticity,  26 , 23-42. 

[3]  Abeyaratne,  R.  and  Hou,  H.-s.  (1991b).  On  the  occurrence  of  the  cavitation 
instability  relative  to  the  asymmetric  instability  under  symmetric  dead-load  conditions. 
Q.  J.  Mech.  Appl.  Math.,  44 , 429-449. 

[4]  Antman,  S.  S.  and  Negron-Marrero,  P.  V.  (1987).  The  remarkable  nature  of  radially 
symmetric  equilibrium  states  of  aeolotropic  nonlinearly  elastic  bodies.  J.  Elasticity,  18  , 
131-164. 

[5]  Ball,  J.  M.  (1982).  Discontinuous  equilibrium  solutions  and  cavitation  in  nonlinear 
elasticity.  Phil.  Trans.  R.  Soc.  Lond.,  A  306 , 557-610. 

[6]  Chou-Wang,  M.-S.  and  Horgan,  C.  O.  (1989a).  Void  nucleation  and  growth  for  a 
class  of  incompressible  nonlinearly  elastic  materials.  Int.  J.  Solids  and  Structures,  25  , 
1239-1254. 

[7]  Chou-Wang,  M.-S.  and  Horgan,  C.  O.  (1989b).  Cavitation  in  nonlinear 
elastodynamics  for  neo-Hookean  materials.  Int.  J.  Engineering  Science,  27 , 967-973. 

[8]  Chung,  D.-T.,  Horgan,  C.  O.  and  Abeyaratne,  R.  (1986).  The  finite  deformation  of 
internally  pressurized  hollow  cylinders  and  spheres  for  a  class  of  compressible  elastic 
materials.  Int.  J.  Solids  and  Structures,  22 , 1557-1570. 

[9]  Chung,  D.  -T.,  Horgan,  C.  O.  and  Abeyaratne,  R.  (1987).  A  note  on  a  bifurcation 
problem  in  finite  plasticity  related  to  void  nucleation.  Int.  J.  Solids  and  Structures,  23  , 
983-988. 

[10]  Ertan,  N.  (1988).  Influence  of  compressibility  and  hardening  on  cavitation.  ASCE  J. 
of  Engineering  Mechanics,  114, 1231-1244. 

[11]  Ericksen,  J.  L.,  and  Rivlin,  R.  S.  (1957)  Large  elastic  deformations  of  homogeneous 
anisotropic  materials.  J.  Rational  Mech.  Anal.,  3 , 281-301. 

[12]  Eubanks,  R.  A.  and  Sternberg,  E.  (1954)  On  the  axisymmetric  problem  of  elasticity 
theory  for  a  medium  with  transverse  isotropy.  J.  Rational  Mech.  Anal.,  3 , 89-101. 


156 


[13]  Gent,  A.  N.  (1990).  Cavitation  in  robber:  a  cautionary  tale.  Rubber  Chem.  Technol. 
63 ,  G49-G53. 

[14]  Gent,  A.  N.  and  Lindley,  P.  B.  (1958).  Internal  rupture  of  bonded  rubber  cylinders 
in  tension.  Proc.  R.  Soc.  Lond.,  A  249 , 195-205. 

[15]  Gurtin,  M.  E.  (1981).  An  Introduction  to  Continuum  Mechanics.  Academic  Press, 
San  Diego. 

[16]  Haughton,  D.  M.  (1986).  On  non-existence  of  cavitation  in  incompressible  elastic 
membranes.  Q.  J.  Mech.  Appl.  Math.,  39 , 289-296. 

[17]  Haughton,  D.  M.  (1990).  Cavitation  in  compressible  elastic  membranes.  Int.  J.  of 
Engineering  Science,  28 , 163-168. 

[18]  Horgan,  C.  O.  (1992).  Void  nucleation  and  growth  for  compressible  nonlinearly 
elastic  materials:  an  example.  Int.  J.  Solids  and  Structures,  29 , 279-291. 

[19]  Horgan,  C.  O.  and  Abeyaratne,  R.  (1986).  A  bifurcation  problem  for  a  compressible 
nonlinearly  elastic  medium:  growth  of  a  micro-void.  J.  Elasticity,  16 , 189-200. 

[20]  Horgan,  C.  O.  and  Pence,  T.  J.  (1989a).  Void  nucleation  in  tensile  dead-loading  of  a 
composite  incompressible  nonlinearly  elastic  sphere.  J.  Elasticity,  21 , 61-82. 

[21]  Horgan,  C.  O.  and  Pence,  T.  J.  (1989b).  Cavity  formation  at  the  center  of  a 
composite  incompressible  nonlinearly  elastic  sphere.  J.  Appl.  Mech.,  56 , 302-308. 

[22]  Horgan,  C.  O.  and  Pence,  T.  J.  (1989c).  Void  nucleation  due  to  large  deformations 
in  nonlinearly  elastic  composites.  Proceedings  of  the  4th  Japan-US  Conference  on 
Composite  Materials,  Washington,  D.C.,  June  1988  (ed.  by  Vinson,  J.  R.)  pp.  232-241. 
Technomic  Publishing  Co.,  Lancaster,  PA. 

[23]  Hou,  H.-s  and  Abeyaratne,  R.  (1992).  Cavitation  in  elastic  and  elastic-plastic  solids. 
J.  Mech.  Phys.  Solids,  40 , 571-592. 

[24]  Hou,  H.-s  and  Zhang,  Y.  (1990).  The  effect  of  axial  stretch  on  cavitation  in  an 
elastic  cylinder.  Int.  J.  Nonlinear  Mechanics  25 , 715-722. 

[25]  Huang,  Y.,  Hutchinson,  J.  W.  and  Tvergaard,  V.  (1991).  Cavitation  instabilities  in 
elastic-plastic  solids.  J.  Mech.  Phys.  Solids  39 , 223-242. 

[26]  James,  R.  D.  and  Spector,  S.  J.  (1989).  The  formation  of  filamentary  voids  in  solids. 
J.  Mech.  Phys.  Solids  39 , 783-814. 

[27]  Jaunzemis,  W.  (1967).  Continuum  Mechanics.  McMillan,  New  York. 

[28]  Meynard,  F.  (1992).  Existence  and  nonexistence  results  on  the  radially  symmetric 


157 


cavitation  problem.  Quart.  Appl.  Math.  50 , 201-226. 

[29]  Ogden,  R.  W.  (1984).  Non-Linear  Elastic  Deformations.  Ellis  Horwood  Limited, 
England. 

[30]  Pericak-Spector,  K.  A.  and  Spector,  S.  J.  (1988).  Nonuniqueness  for  a  hyperbolic 
system:  cavitation  in  nonlinear  elastodynamics.  Arch.  Rat.  Mech.  Anal,  101 , 293-317. 

[31]  Podio-Guidugli,  P.,  Vergara  Caffarelli,  G.  and  Virga,  E.  G.  (1986).  Discontinuous 
energy  minimizers  in  nonlinear  elastostatics:  an  example  of  J.  Ball  revisited.  J. 
Elasticity,  16 , 75-96. 

[32]  Sivaloganathan,  J.  (1986a).  Uniqueness  of  regular  and  singular  equilibria  for 
spherically  symmetric  problems  of  nonlinear  elasticity.  Arch.  Ration.  Mech.  Analysis, 
96, 97-136. 


[33]  Sivaloganathan  J.  (1986b).  A  field  theory  approach  to  stability  of  radial  equilibria  in 
nonlinear  elasticity.  Math.  Proc.  Camb.  Phil.  Soc.,  99 , 589-604. 

[34]  Sivaloganathan,  J.  (1991).  Cavitation,  the  incompressible  limit,  and  material 
inhomogeneity.  Quart.  Appl.  Math.  49,521-541. 

[35]  Spencer,  A.  J.  M.  (1972).  Deformations  of  Fibre-reinforced  Materials.  Oxford 
University  Press,  Oxford. 

[36]  Spencer,  A.  J.  M.  (1982).  The  formulation  of  constitutive  equations  for  anisotropic 
solids.  In  Mechanical  Behavior  of  Anisotropic  Solids  (Proceedings  of  Euromech 
Colloquium  115,  ed.  by  J.-P.  Boehler)  pp.  3-26.  Martinus  Nijhoff,  The  Hague. 

[37]  Steigmann,  D.  J.  (1992).  Cavitation  in  elastic  membranes  —  an  example.  J. 
Elasticity  28 , 277-287. 

[38]  Stuart,  C.  A.  (1985).  Radially  symmetric  cavitation  for  hyperelastic  materials.  Ann. 
Inst.  Henri  Poincare- Analyse  non  lineaire,  2 , 33-66. 

[39]  Stuart,  C.  A.  (1993).  Estimating  the  critical  radius  for  radially  symmetric  cavitation. 
Quart.  Appl.  Math.  51  (in  press). 

[40]  Tian-hu,  H.  (1990).  A  theory  of  the  appearance  and  growth  of  the  micro-spherical 
void.  Int.  J.  of  Fracture,  43 ,  R51-R55. 

[41]  Truesdell,  C.  and  Noll,  W.  (1965).  The  Non-Linear  Field  Theories  of  Mechanics. 
Handbuch  Der  Physik,  m/3.  Springer- Verlag,  New  York. 

[42]  Tvergaard,  V.  (1990).  Material  failure  by  void  growth  to  coalescence.  In  Advances 
in  Applied  Mechanics  (ed  by  J.  Hutchinson  and  T.  Wu),  vol.  27,  pp.  83-151.  Academic 
Press,  San  Diego. 


158 


[43]  Tvergaard,  V.,  Huang,  Y.  and  Hutchinson,  J.W.  (1992).  Cavitation  instabilities  in  a 
power  hardening  elastic-plastic  solid.  European  J.  of  Mechanics  A  (Solids)  11 , 215-232. 

[44]  Walker,  J.  L.  (1956).  Structure  of  ingots  and  castings.  In  Liquid  Metals  and 
Solidification,  pp.  319-336.  American  Society  for  Metals,  Cleveland. 

[45]  Williams,  M.  L.  and  Schapery,  R.  A.  (1965).  Spherical  flaw  instability  in 
hydrostatic  tension.  Int.  J.  Fracture  Mechanics,  1 , 64-71. 


